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§.1 Introduction 

The development of processes to produce fullerenes and carbon nanotubes has 
largely been empirical. Fullerenes were first discovered in the soot produced by laser 
ablation of graphite [1] and then in the soot of electric arc evaporated carbon [2]. 
Techniques and conditions for producing larger and larger quantities of fullerenes 
depended mainly on trial and error empirical variations of these processes, with attempts 
to scale them up by using larger electrodes and targets and higher power. Various 
concepts of how fullerenes and carbon nanotubes were formed were put forth, but very 
little was done based on chemical kinetics of the reactions. This was mainly due to the 
complex mixture of species and complex nature of conditions in the reactors. 
Temperatures in the reactors varied from several thousand degrees Kelvin down to near 
room temperature. There are hundreds of species possible, ranging from atomic carbon to 
large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to 
complexes of metals and carbon. Most of the chemical kinetics of the reactions and the 
thermodynamic properties of clusters and complexes have only been approximated. In 
addition, flow conditions in the reactors are transient or unsteady, and three dimensional, 
with steep spatial gradients of temperature and species concentrations. All these factors 
make computational simulations of reactors very con p'ex and challenging. 

This article addresses the development of the chemical reactions involved in 
fullerene production and extends this to production of carbon nanotubes by the laser 
ablation/oven process and by the electric arc evaporation process. In addition, the high- 
pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. 
The first one addresses the thermochemical aspects of modeling; and considers the 
development of chemical rate equations, estimates of reaction rates, and thermodynamic 
properties where they are available. The second part addresses modeling of the arc 
process for fullerene and carbon nanotube production using 0-D, 1-D and 2-D fluid flow 
models. The third part addresses simulations of the pulsed laser ablation process using 
time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser 
ablation process. The fourth part addresses steady state modeling in 0-D and 2-D of the 
HiPco process. In each of the simulations, there is a variety of simplifications that are 
made that enable one to concentrate on one aspect or another of the process. There are 
simplifications that can be made to the chemical reaction models, e.g. reduction in 
number of species by lumping some of them together in a representative species. Other 
simulations are canned out by eliminating the chemistry altogether in order to concentrate 
on the fluid dynamics. When solving problems with a large number of species in more 
than one spatial dimension, it is almost imperative that the problem be decoupled by 
solving for the fluid dynamics to find the fluid motion and temperature history of 
“particles” of fluid moving through a reactor. Then one can solve the chemical rate 
equations with complex chemistry following the temperature and pressure history. One 
difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to 
take dilution and mixing into account. This changes the ratio of carbon species to 
background gas. Commercially available codes may have no provision for including 
dilution as part of the input. One must the write special solvers for including dilution in 
decoupled problems. 
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The article addresses both fullerene production and single-walled carbon 
nanotube (SWNT) production. There are at least two schemes or concepts of SWNT 
growth. This article will only address growth in the gas phase by carbon and catalyst 
cluster growth and SWNT formation by the addition of carbon. There are other models 
that conceive of SWNT growth as a phase separation process from clusters made up 
carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We 
will not deal with that concept in this article. Further research is needed to determine the 
rates at which these composite clusters form, evaporate, and segregate. 

§.2 Reaction Schemes and Thermodynamic Properties 

This section deals with the formation of fullerenes and SWNTs via precursor 
species, clusters of carbon, and metal catalysts. Two basic processes are considered; one 
being processes that involve carbon vaporized at high temperatures such as the arc 
process, and the laser ablation process. The other is a gas phase process in which carbon 
is supplied by carbon monoxide at high pressure, the HiPco process. 

§.2. 1 Carbon and Fullerene Kinetics 

Various workers have studied the kinetics of the formation of carbon clusters. 
Bemholc and Phillip . [3] solved the Smolnchowski 'quitions for carl *n cIum-Cis C„ up to' 
n=25 using a reaction probability, or kernel, based on a scaled derivative of the Gibbs 
free energy. Their work included both negative and positive ions, in addition to neutral 
clusters. Using this technique, they predicted the distribution clusters, including the 
presence of magic numbers. Creasy and Brenna [4,5] used a simple model in which 
clusters grew by attachment of only the small clusters C, C2, and C3. The reaction rates 
were estimated from gas kinetic rates. Subsequently, Creasy [6] developed more complex 
reaction models that allowed for multiple steps, in an attempt to account for “magic 
number” clusters. His analysis produced cluster and fullerene distributions up to about 
n=450 that depended on the initial concentration of small carbon molecules. The higher 
the initial density, the greater was the average size cluster, and the more the distribution 
shifted toward larger clusters at steady state. The reaction rates in their model did not 
depend on temperature. 

§.2.1.1 Model of Krestinin and Moravsky 

Krestinin, et al. [7] developed a reaction scheme for fullerene formation. Krestinin 
and Moravsky [8] applied the model to an arc process, where they achieved reasonable 
agreement with measured fullerene production, considering the approximations made in 
the arc flow field. They were able to explain the relatively constant ratio of C60 to C70 
seen in arc experiments. The objective of their model was to depict carbon vapor 
condensation in an arc reactor. Because of the lack of thermodynamic data for the various 
clusters, the model was formulated, to the extent necessary, not to rely on that data. Many 
reactions are written as separate forward and backward reactions, or only in terms of 
condensation and not decomposition of clusters. They rejected formulating a scheme that 
includes all carbon clusters, as well as their ions, due to the complexity of experimental 
verification. As experimental data becomes available and simulations are made, then the 
model can be validated and rate coefficients and reactions can be refined. Their rates are 
based on measurements using ionic gas chromatography published in the literature 
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[9,10,1 1,12]. Carbon clusters up to n=31 are assumed to be highly reactive chains, cycles, 
and polycycles. Clusters with n=32 to 79 are closed shells, where the main growth or 
decomposition is from C2 attachment or detachment, respectively. That odd-n clusters are 
less stable than even clusters is attributed to a lack of resonance stabilization. Their 
model reflects a 251-kJ/mole lower enthalpy of formation of even-n numbered clusters as 
compared with odd-n numbered clusters. Upon decomposition, odd clusters tend to eject 
C, rather than C2 [13]. Clusters with n=80 or greater are assumed to be soot Z. In the 
calculations of Krestinin, et al. [7] they assumed that soot takes on a size distribution, and 
that soot reactions are heterogeneous. The reaction rate then depends on the surface area 
of the soot particle. (In later applications [14] of this model to nanotube production 
discussed in this article, the model is modified so that soot is treated as a single gas phase 
species.) 

Reaction equations and rate coefficients from [7] are given in Arrhenius form in 
Table §.l for the fullerene model. 


k = A exp(-E/RT) (§.l) 

The thermodynamic properties for these clusters were determined from existing 
data for C, [15], €2,' [16], C 3 - C r n, [17], and-C^p and C70F, [18]. The enthalpy of- • 
formation for clusters with n— L 1 -59, 61 69, and 71 79 were estimated by interpolating 
the values between n=10 and 60. The entropy of formation was estimated by interpolation 
between n=5 and 60. The specific heats were set at values corresponding to T>2000 K 
where fully excited vibration and rotation is achieved. Then Cp/R=3n-2 for n>6. The 
complete set of thermodynamic coefficients for the fullerene model is listed in Table §.2. 
These coefficients are in the old NASA format where the various thermodynamic 
functions are calculated by the formulas: 


Cpi/R = A„ + A 2i T + A 3i T 2 + A* T 3 + A Si T 4 (§.2) 

hj/RT = A h + A 2i T/2 + A 3i T 2 /3 + A 4i T 3 /4 + A 5i T 4 /5 + A 6l /T (§.3) 

Sj/R = A,jln(T) + A 2i T + A 3i T 2 /2 + A 4i T 3 /3 + A 5i T 4 /4 + A 7i (§.4) 

where R is the universal gas constant, hj is the specific enthalpy of the i lh species, and Sj is 
the specific entropy. For the clusters for which thermodynamic properties were 
interpolated, it can be seen from the values in Table §.2, that C P j and Sj are constant with 
respect to temperature. Thus, they are rather approximate. In most simulations, this does 
not pose a problem because of dilution. The fully excited approximation for C p is not 
very significant due to dilution and the high temperatures used in production of fullerenes 
and nanotubes. In such case the energy equation for the flow does not depend 
significantly on the energetics of the carbon reactions. 

§.2.1.2 Reduced Carbon Cluster Scheme 

Because of the significant computer time required for models containing large 
numbers of species used in computational fluid dynamic (CFD) simulations, large models 
are impractical for simulating complex time-dependent flows, such as, pulsed laser 


3 


ablation. To overcome this limitation of CFD codes it is almost imperative that the 
number of species in the chemical reaction models be reduced to a manageable number. 
An effort to do this was accomplished for simulating the HiPco process and likewise for 
the carbon fullerene and carbon/metal-vapor SWNT processes. Gokgen, et al. [19] 
showed that in simulations of the flow in the HiPco reactor the number of catalyst 
particles Fe n for n>2 could be lumped into a single representative iron cluster Fe c , and the 
200 SWNTs in their model could be reduced to a single representative nanotube, CNT. 
They compared the production of CO2 using both models and showed that they produced 
similar results for a constant temperature case. This enabled Gokgen, et al. to use the 
model in their simulations of the HiPco reactor, thus saving much computer time. Their 
full and reduced models will be discussed in more detail later. 

In a similar way we have reduced the fullerene model of Krestinin, et al. [7], In 
the reduced fullerene model we have lumped all carbon cluster C n with n>3 into a single 
representative species CC. For carbon balance in the reaction equations, we specify the 
number of carbons in CC to be 40. The feedstock for growth of CC and the fullerene C60 
are the small carbon molecules, C, C2, and C3. These small molecules coalesce to form 
CC, which then forms C6o- A soot Z is also formed from CC. This soot is considered a 
single species of constant n=80. The reduced fullerene model is given in Table §.3, and 
the corresponding thermodynamic coefficients are given in Table §.4. The .... ... ,, ... . ... 
thermudynainic coefficient ; for CC and’Z conesp'Oi d ‘o C40 and Cg) 'in the full fuherene 
model. These coefficients were derived from interpolations of properties as indicated 
above. 

§.2.2 Metal Catalyst Schemes 

When pure carbon is used as the feedstock for reactions in high temperature 
processes, the product is usually amorphous carbon chains and rings, fullerenes, graphitic 
particles, and soot. At some conditions multiwall carbon nanotubes (MWNT) are formed. 
However, if some metal catalysts are added, there is a preference to form single-wall 
carbon nanotubes. SWNTs have many interesting and unique properties. They may be 
metallic or semi conducting, depending on their chirality and diameter. The present 
understanding of their formation does not allow us to predict from chemical kinetics their 
relative abundance in a given production process. To these authors’ knowledge, no one 
has developed a chemical kinetics model of particular types of nanotubes. Therefore, in 
the following models, carbon nanotubes will be treated as a single species, with no 
consideration of types. At present, there are several concepts of how SWNTs form; but 
we will consider here only those that depend on gas phase reactions in which metal 
clusters are involved. It is assumed in this article that gas phase reactions govern the rate 
of the formation of SWNTs and not diffusion or segregation of carbon from metal 
catalyst particles. 

Metal catalyst clusters are formed from metals vaporized from electrodes in the 
arc process, from composite targets in laser ablation processes, or from decomposition of 
metal-containing precursor gases, such as iron pentacarbonyl. 

§.2.2.1 Iron Cluster Formation and Evaporation 

Several authors have developed the kinetics of iron cluster formation and 
evaporation. Girshick and his coworkers published a series of articles on iron cluster 
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nucleation and growth [20,21,22,23,24,25,26]. They used modified classical aerosol 
cluster theory to develop nucleation, growth and evaporation rates. They compared 
plasma generated iron atom nucleation and cluster evolution with measuremenfs. Their 
basic equations for the rates of cluster growth are given by the collision frequency. 
Evaporation is considered as a function of cluster size by accounting for its volume and 
surface tension. They expressed the net rate of change of iron monomers as 


Dn ] 

Dt 


Pg ;=i 




(§.5) 


while for dimers and all larger clusters they write 

P 


~~ = E Pa n k n j ~ PPh YjPjk n i + R Li n k + i ~ K n k 

Dt Z /+ j=k 


(§- 6 ) 


i+ j-k 


In these equations R is the rate of monomer generation/depletion by chemical processes, 
nj is the number of j-mers per unit mass of gas, (3jj is the collision frequency for i-mers 
with j-mers, Ej is the evaporation coefficient. Adding the Kronecker delta § 2 j accounts for 
two monomers formed from each dimer. The collision frequency function is taken from 
ideal gas kinetic theory for the free molecule regime 
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The evaporation coefficient is written 

R j = fitj- i n < ex p(®L/ 2 ~ O' - O'" J ) 

where n s is the equilibrium saturation number density, © is a dimensionless surface 
energy 

0 = ^L 

kT 

where a is the surface tension, and Si is the surface area of the monomer. 

The saturation number density n s can be calculated using the perfect gas law and 
the vapor pressure 

f-E 

Ps = P°s ex P 


(§•7) 

(§• 8 ) 

(§•9) 


kT 


(§•10) 


where p s °is the coefficient of the vapor pressure expression. 

Rate coefficients for gas phase kinetics can be derived in Arrhenius form from 
these expressions for (3;j and Ej. The evaporation rate coefficient can be expressed as 


kj = AjT ~ 112 exp 


J_ 

kT 


where 

and 
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Aj = r x 
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The rate coefficients developed from this analysis were published for iron clusters in [" ]. 


§.2.3 High Pressure Carbon Monoxide Disproportionation (HiPco) Models 


The high-pressure carbon monoxide (HiPco) process begins with iron 
pentacarbonyl diluted in carbon monoxide. This mixture is injected at near room 
temperature into a high temperature reactor, where it mixes rapidly with hot carbon 
monoxide. As the temperature increases above about 600 °C, Fe(CO )5 decomposes into 
FeCO and Fe, which then react to form dimers of iron. These iron dimers coalesce with 
iron atoms and other dimers to form Fe 3 and Fe 4 . Coalescence continues and larger and 
larger clusters of iron are formed. 

For modeling the HiPco reactor, Dateo, et al. [ ] adopted the rate coefficients of 
Krestinin, et al. [29] for iron cluster condensation and evaporation. These coefficients are 
based on assumptions, such as, clusters up to n=4 are considered to be in the gas phase, 
above that they are considered to be condensed phase with rate coefficients for reactions 
with smaller clusters having reached their asymptotic limit. Evaporation is assumed to be 
negligible for clusters larger than n=9; and condensation and growth of iron clusters 
occurs mainly through reactions with Fe 4 and smaller clusters. The rate coefficients that 
they used were developed through considerations of the dependence of equilibrium 


constants on the formation energy and standard reduced thermodynamic potential. They- 
accounted for the’ collision frequency and its dependence on size when computing the rate 
coefficients for cluster growth. These reactions and rates from [29] are given in Table 5. 

After Krestinin, et al. [29] published their results for iron cluster chemistry, 
Vlasov, et al. [30] developed rate coefficients for iron cluster kinetics that takes into 
account nonequilibrium effects associated with the difference between the gas 
temperature and the internal temperature of the iron clusters. They made their 
calculations of decomposition of clusters within the framework of the statistical theory of 
unimolecular reactions, in which they treated only the attachment and detachment of 
monomers. The nonequilibrium effects are manifested by depletion of internal energy by 
detachment and gain of internal energy by attachment of monomers. Nonequilibrium 
effects are seen most vividly for clusters up to n=20. After that, there seems to be a 
smooth variation in the equilibrium constant and rate coefficients. They plotted their 
equilibrium constants, decomposition rate coefficient, and recombination coefficients as 
functions of cluster size for temperatures from 1000 K to 3000 K in Figs. 3, 4 and 5, 
respectively, of their article [30], One present author (CDS) curve fit these figures and 
obtained Arrhenius expressions for the rate coefficients. These coefficients for n<51 are 
given in Table §.6 5. For n>50 the constants in these expressions can easily be 
extrapolated. A comparison of evaporation rate coefficients at T=1380 K is given in Fig. 
§.l; and recombination coefficients at T=1500 K are given in Fig. §.2. It is apparent there 
is a strong influence nonequilibrium effects associated with atom attachment and 
evaporation (Vlasov), particularly for small clusters. The recombination rates of Girshick, 
where classical gas kinetics governs attachment, are probably overestimated for small 
clusters. 


§.2.3.1 Decomposition of Fe(CO) 5 

There are several sources of reaction rates for Fe(CO)s decomposition. Krestinin, 
Smirnov, and Zaslonko [29] included a two-step decomposition scheme when simulating 
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the decomposition of Fe(CO )5 and the condensation of Fe. Their decomposition 
reactions are: 

Fe(CO) 5 -4 FeCO + 4CO (§.14) 

FeCO + M — > Fe + CO + M (§.15) 

The rate coefficient for (§.14) was taken from [31]. They inferred the rate 
coefficient for (§.15) from the shock tube data. Their complete reaction set is given in 
Table §.7 6. Note that the formation of dimers in their scheme is allowed by the reactions 

FeCO + FeCO -> Fe 2 + 2CO (§.16) 

FeCO + Fe -> Fe 2 + CO (§.17) 

Since FeCO can react to form Fe 2 , the direct nucleation of Fe 2 from two Fe atoms is not 
required. Due to the nature of the closed shell electronic structure of Fe it is unlikely that 
two ground state atoms are likely to form Fe 2 . It probably requires one or both atoms to 
be in excited states to achieve Fe recombination [27]. The rate coefficient for the reaction 

Fe + Fe + M -> Fe 2 + M (§.18) 

given in Table §.5 6, was inferred from shock tube data. Since formation of iron clusters 
is very rapid in the shock tube experiments, it is possible that reactions (§.16) and (§.17) 
dominate. This was discussed in [27]. 

Smirnov inferred iron pentacarbonyl decomposition rates from shock tube 
measurements X32] of ab&dfptioh's^ \ T Fe atoms and parent n-oiefciiles 

Vlasov, et al. [33] included them in a review oi metal compound decomposition. 
However, there was a difference in the Fe(CO )5 decomposition rate of one order of 
magnitude. It is not clear which rate is correct. Rumminger, et al. [34] compiled a set of 
reactions for the inhibition of flames by iron pentacarbonyl. They included recombination 
reactions from Seder, Ouderkirk, and Weitz [35]. Using more recent bond energy 
measurements from Sunderlin, et al. [36] and accurate ab initio calculations of Ricca 
[37], Dateo, et al. [28] calculated the equilibrium constants for 

decomposition/recombination reactions of Fe(CO) n . With the association rates of Seder, 
et al. [37] they calculated the dissociation rates that are denoted “Ames” in Table §.7. 
Since, the HiPco reactor that they modeled is at many atmospheres pressure, they did not 
give rate coefficients for the low-pressure fall off regime. Rate coefficients for the iron 
carbonyl dissociation reactions of Smimov [32] are also given in Table §.7 6. 

Iron clusters are formed from FeCO and Fe combining and agglomerating. 
Krestinin, et al. [29] developed a set of iron cluster reactions and thermodynamic data. 
Their reactions and rate coefficients are given in Table §.5. They contend that due to low 
pressure, FeCO decomposition is controlled by the bimolecular activation step. Thus, the 
rate of decomposition of FeCO is commensurate with the reactions 

FeCO + FeCO => Fe 2 + CO + CO (§.19) 

FeCO + Fe => Fe 2 +CO (§.20) 

Because of this, Fe 2 and larger clusters are produced in large quantities during FeCO 
decomposition. Therefore, it is not necessary that direct formation of Fe 2 nucleation occur 
via Fe-atom recombination. Further growth of iron clusters is dominated by coagulation 
of these clusters. Evaporation of small clusters was estimated based on classical theory. 
However, These rates did not lead to good agreement with measurements. Therefore, 
they were adjusted as given in Table §.5. The original rate coefficients are given in 
parentheses. 


7 


§.2.3.2 CO Attachment to Iron Clusters 

In the development of the model for the high-pressure carbon monoxide reactor 
Dateo, et al. [28] assumed that CO attaches to iron clusters beginning at Feio- There is 
some data [38,39] that suggests that CO will attach to smaller clusters and may attach to 
larger clusters in greater numbers. However, multiple bonding of CO probably does not 
alter significantly the rate of production of carbon nanotubes (Boudouard reaction). To 
keep their model simpler, only Fe n CO are considered for n>9 and no multiple 
attachments of CO to the clusters. This reaction forms the first of a three-step model for 
the production of carbon nanotubes. 

§.2.3.3 Formation of Carbon Nanotube Reactions 

The second of the three steps for the formation of carbon nanotubes in the Dateo 
model is detachment of CO from Fe n CO. The final step is the conversion of Fe n CO to a 
carbon nanotube. 

Fe n CO + CO-> J3. Fc n C m + ^ C0 2 + (l - £)Fe n CO (§-21) 

where (3=l/(2nNT-l) and nwT is the average number of carbon atoms in a carbon 
nanotube formed in Fe p -C^T. In the calculations of Dateo, et al. [28], they let ont= 3000. 

In the calculations - of-Scott, et al. [27 j, they let hnt 999 due to a limitation, in. the input. «... ... 

format of CHEMKIN® [40], 

§.2.3.4 Inerting or “Death” of Clusters 

Daniel Colbert of Rice University suggested by that carbon nanotubes stop 
growing due to lack of iron clusters or due to carbonizing of the iron cluster to which the 
nanotube is attached during growth. Iron clusters may become over coated with carbon 
that blocks the formation of nanotubes. This results in their “death” as a useful catalyst; 
and nanotube growth stops. To account for this possibility, the following two types of 
reactions are included in the Dateo model 

Fe n — Fe n * (§.22) 

Fe n C N T-+Fe n *C NT (§.23) 

Fe n and Fe n Cnt are clusters and nanotubes that have been poisoned with carbon to the 
extent that they cannot catalyze the reactions and cease to grow. The HiPco reaction 
model of Dateo, et al. [28] is summarized in Table §.8. 

§.2.3.5 Reduced HiPco Model 

The full model of Table §.8 is very large, having 971 species and 1948 reactions. 

To simplify the model and make it useable in a 2D CFD code, Dateo, et al. [28] collapsed 
or lumped all the iron clusters having n>2 into a single representative cluster Fee- They 
also collapsed the other clusters as well, leaving the following fourteen species: iron 
carbonyls (Fe(CO) 5 , Fe(CO) 4 ,Fe(CO) 3 ,Fe(CO) 2 , and FeCO), CO, C0 2 ,Fe, Fe 2 ,Fe c , 

FecCO, FecCNT. Fee , and Fcc’Cnt- Their reduced model, listed in Table §.9, has only 
twenty-two reactions. 

§.2.4 Carbon Vapor Models of Carbon Nanotube Formation 

Carbon vapor models are required for processes that start from vaporized carbon 
such as in laser ablation processes and the arc vaporization process. These processes, and 
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how these models can be applied to them, are described in more detail in the subchapters 
of this chapter. Single wall carbon nanotubes require a feedstock of carbon vapor plus 
vaporized metal catalysts. From this feedstock clusters of carbon and nickel condense to 
form SWNTs and impurities, such as soot and carbon-coated metal clusters. The general 
concepts of SWNT formation fall into two or three categories. In order to form carbon 
nanotubes it is necessary for the carbon to be catalyzed by nickel and other metals. 
Nucleation and growth of carbon nanotubes are not well understood; and there are many 
possible scenarios for carbon and nickel to condense in a reactor and then to grow 
nanotubes. An early scheme, called the “scooter mechanism” was envisioned by 
Smalley’s group at Rice University. In it, carbon was envisioned to form fragments of 
graphene sheets that would tend to close into fullerenes. However, nickel or other metal 
catalyst atoms would attach to an open edge of the carbon cluster and “scoot” or move 
around, holding open the end of the tube, while small carbon clusters, e.g. C 2 would 
attach to the nascent nanotube that would continue to grow. In this scheme, carbon would 
condense first, then nickel would tend to agglomerate until a cluster of nickel at the end 
of the nanotube are so large as to soot up and cease to allow the nanotube to grow; or the 
metal cluster would detach from the nanotube, thus ending its growth. This scenario is 
supported by the observation that laser ablation nanotubes are almost never seen attached 
to metal clusters. It also conforms to the tendency for carbon to agglomerate first, 
follow ed by nickel, ‘ due to the stfongei to the higher 
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concentration of carbon than catalysts. 

Another scheme that has much support involving vaporized carbon in SWNT 
production is one in which metal catalysts are envisioned to form along with carbon with 
the clusters at high temperature. As the mixture cools, carbon precipitates out of the 
cluster and the energetics and geometry favor the formation of nanotubes. This is 
essentially a phase change phenomenon due to cooling of a mixture, especially eutectic 
mixtures of nickel, cobalt, and carbon. This scheme is supported by the observation that 
the ends of nanotubes, particularly bundles, are sometimes seen attached to metal 
clusters. Some clusters are seen to have a number of bundles attached in sort of a “sea 
urchin” pattern. See Fig. §.3, [courtesy of Alexander Moravsky, MER Corp.] and Fig. 4 
from [41]. 


a « 


§.2.4.1 Nickel Cluster Nucleation and Growth 

Nickel is used as a catalyst in both the laser ablation process and the arc 
vaporization process. Other metal catalysts that are commonly used are cobalt and 
yttrium. Since cobalt is a neighboring transition metal it has similar properties to nickel. 
Yttrium properties are less well known. For these reasons, only nickel will be discussed 
here. To assess why nickel did not produce nanotubes in the HiPco reactor, Scott and 
Smalley [42] tested the effect of the diatomic bond energy and the bond energy of nickel 
atoms with CO. The bond energy of Ni to Ni is larger than that of Fe to Fe and also, the 
bond energy between CO and Ni is larger than between CO and Fe. To assess the effects 
of bond energies, simulations of the HiPco reactor were made using each metal as 
catalyst. To do the nickel catalyst simulation, it was necessary to create agglomeration 
and evaporation rates for nickel. The Girshick technique described in equations (§.7-§.13) 
was used. The surface tension was adjusted until the activation energy for Ei equaled the 
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Ni2 bond energy. The activation energy for NiCO was set equal to its corresponding bond 
energy. 

In principle, more accurate evaporation rate coefficients for small n can be 
calculated from the cross section data of Amentrout’s group [43,44], They measured the 
dissociation of Ni atoms and small clusters from nickel clusters from n=2 to 18 in a 
crossed beam apparatus in which they obtained collisions of the nickel clusters ions with 
xenon. They corrected their data for internal energy (pressure) effects; and they also 
suggested the following simple equation to fit their data. 

ct(e)= — — — — ~ ^ (§.24) 

E 


where N and CTo are adjustable parameters, E is the relative kinetic energy, E; is the 
internal energy of the cluster, and Eo is the collision induced dissociation threshold 
energy. This provided an accurate and convenient formula for calculating rate 
coefficients from the expression 


k.. = 



exp 


2kT 


c{\)d\ 


(§•25) 


where k e is. the evaporatiomrate.coeffiejent,. ji is the reduced mass, v = ^J2E isihe, - * 
relative velocity, and cr(v jTis the collision cros's-sectiori as a function of relative velocity. 


The main difficulty with this method is the need for the internal energy, which is not 
known for many clusters. 

Because of the large number of possible nickel clusters, there seemed to be a need 
for reducing the model in some way. The option chosen [42] was to lump clusters greater 
than n>8 into sizes 16, 32, ... 2048. We call this the “binary” model. Evaporation 
equations for these clusters are written 

Ni n — > Nin/2 + n/2 Ni (§.26) 

which, in an approximate fashion, accounts for the nickel clusters not included in the 
model. Agglomeration rates for all nickel clusters and evaporation rate coefficients for 
n>16 were determined from Girshick’s method, equations (§.7-§.13). These evaporation 
coefficients are given in Table §.10a; and the agglomeration coefficients are given in 
Table §. 10b. The rate coefficients for n>16 have been multiplied by n to account for the 
missing clusters between each n and 2n. 


§ 2.4.2 Combined Carbon and Nickel Clusters to Form Carbon Nanotubes 

One kinetics model that takes advantage of the fullerene model of Krestinin, et al. 
was proposed by Scott [14], The nickel model of Scott and Smalley [42] is added to the 
set of reactions for carbon to produce fullerenes and other large carbon clusters is. 

Carbon and nickel clusters combine to form the nuclei for the growth of carbon 
nanotubes. The rates of nickel cluster formation and growth are as discussed previously, 
the growth of the combination carbon/nickel clusters can be obtained from the aerosol 
theory of Girshick as given in equation (§.7). Carbon clusters in the form of soot Z, and 
fullerenes C6o and C70 are assumed to be the precursors that combine to form 
nickel/carbon clusters, ZNic and CFNic- These clusters are the species that nucleate and 
grow nanotubes CNT. These reactions are given in Table §.11. For specificity, soot is 
assumed to have 80 carbon atoms, while ZNic is assumed to have 80 carbon atoms and 
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1024 nickel atoms. CNT is assumed to have 999 carbon atoms. Whereas, this is much 
fewer than in actual nanotubes, it can represent them; and the number of them will be 
correspondingly larger. The number 999 is purely arbitrary and can be selected to be 
something greater. However, in [42] it is chosen because the computer code used to solve 
the kinetic rate equations has a three-digit limit in the input format for number of atoms 
in the species definition. 

§.2.4.3 Reduced Combined Carbon and Nickel Cluster Nanotube Scheme 

As we saw in the reduced HiPco model, we can attempt to reduce the number of species 
and reactions by lumping most of the metal clusters as well as carbon clusters into 
representative clusters. Both active and inert (dead) clusters that combine the catalyst and 
carbon can also be lumped together into a single representative clusters. The species in 
this reduced model are Ar/He, C, Cz, C3, Cc, Z (soot,) C 6 of> Ni, N 12 , Nh, Nic, ZNic, DZ, 
and DZNic- The subscript C on these symbols denote “cluster.” The “D” on the symbol 
indicates a dead particle that no longer reacts. As in the reduced HiPco model and the 
reduced fullerene model, we the representative lumped clusters have specified numbers 
of atoms each, just for the specificity required to satisfy element balance in the reaction. 

•'Reaction rate coefficient? may-tmve'tobe adjusted slightly to make the results^ ~- 
agree with those of the full model. This can be attempted by solving the problem 
parametrically, using a simple temperature and dilution profile that is representative of a 
given reactor or process. The reduced model for SWNT production based on carbon 
vapor and metal catalysts is given in Table §.12. Nickel is the metal chosen because it is 
most often used. 

§.2.4.4 Results Comparing Full and Reduced Fullerene Models 

To assure that a reduced model actually is capable of predicting the production of 
nanotubes one can calculate both models for some representative simple cases such as a 
temperature profile linear in time. A comparison of the complete Krestinin, et al. [7] 
model with the reduced fullerene model was presented in [14]. The time evolution of a 
carbon/nickel vapor from 3500 K to 1500 K was calculated using a zero-dimensional 
code (AURORA code of CHEMKIN® package [40]) for both models. The calculation of 
the production of carbon nanotubes and other clusters agree very well. Figs. §.5 and §.6 
show comparisons of the time evolution of various species and the production of SWNTs. 
Whereas, the comparison is very good (within no more than 50%) the production of 
carbon nanotubes and other clusters must be compared with experiment to calibrate the 
rate coefficients used in the models. 

The models discussed in section §.2 have been applied to various problems 
associated with the production of carbon nanotubes. The application of these models has 
taken various levels of fidelity and sophistication. Assessments of the models, 
particularly comparisons between full and reduced models are given in this section. 
Section §.3 contains results applied to the arc process, section §.4 addresses laser ablation 
process results, and §.5 deals with the HiPco process. Two approaches are considered for 
solving the coupled fluid dynamics equations and the chemical kinetics of the models. 
Since complete fluid dynamics in three dimensions plus large models having hundreds of 
species presents a formidable computational task, we approach the solution by 
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simplifications of two types. One is to reduce the size of the model by lumping many 
clusters into a single representative cluster, as in the “reduced” models. One can then use 
these reduced models in highly complex geometrical situations using multidimensional 
fluid dynamics equations. On the other hand, if high fidelity of chemical species and 
reactions is desired, then the chemistry and fluid dynamics solutions can be decoupled. 
The fluid dynamics equations are solved without chemical reactions; and then the 
complex chemical kinetics is developed along streak lines. 


§.3 Analysis and Modeling of Carbon Arc Reactors 

The arc discharge, developed initially for fullerene synthesis [2], can also produce 
single walled carbon nanotubes (SWNT) [ 45 , 46 , 47 , 48 ], by simply adding catalysts to the 
graphite and by changing process conditions. Compared with other synthesis techniques, 
the arc leads to temperatures higher than 5000 K, as measured by optical emission 
spectroscopy [49,50,51], thus insuring total and fast vaporization of graphite anode 
containing bimetallic catalysts in a background gas of helium, [46] argon [46, 32 ], or 
mixtures of them [52, x ’]. Catalysts that are commonly used are nickel and yttrium at 
about 'one f6T6ufW61'af pbfceht ['.6], Opffrnail cdtiil Tons ne 6bt.-l.ncd with this Catalyst *' 
composition [46], for helium at a pressure of 660 mbar, and a current of 100 A. 

Nanotubes are collected in a soot, called “collaret,” that is close to the cooled cathode. 
They seem to grow either in bundles or as individual tubes that coalesce readily into 
bundles. Laboratory scale reactors can produce up to one gram of this collaret per 10- 
minute batch run. The as-produced material contains nanotubes that are 1 to 1.4 nm in 
diameter and a few micrometers in length. Due to its relatively low cost, the arc discharge 
process is used to produce nanotubes industrially. Although prices are still high, they are 
in constant decrease. For example, the MER Corporation in the United States offers a 
soot synthesized by the electric arc method (Fig. §.7) that contains 10 to 40% by mass of 
nanotubes at $50 per gram, compared with $2000 per gram a few years ago! Contrary to 
the competing process by chemical vapor deposition CVD in fixed [53, 55 , 5 VW 9 ] or 
fluidised beds [ 60 , 61 ], nanotubes produced by arc are reputed to have excellent structural 
properties. They are quite straight and present fewer topological defects, probably 
because of the fast growth conditions. However, the major inconvenience of the arc 
process remains the presence of impurities such as catalysts and amorphous carbon mixed 
with the nanotubes. To take advantage of the good structural quality of nanotubes 
produced by arc, two solutions are possible. The first concerns the development of post- 
synthesis purification methods to remove these impurities. However, due to the large 
variety of carbon impurities (Fig. §.8), such as amorphous carbon, graphitic 
nanoparticles, polyhedral carbon particles, onions, single wall nanohoms (SWNHs), 
fullerenes, and some other forms of carbon with metal catalyst enclosed, purification is a 
very tedious task and no purification technique is universally admitted today. In addition, 
this large variety of impurities renders the quantitative techniques developed to measure 
the amount of carbon nanotubes in the sample such as thermogravimetric analysis (TGA) 
or near-infrared (NIR) spectroscopy [ 62 ] not very reproducible. For these reasons, it is 
better to improve the process to increase nanotube yield during their formation. Hence, 
many experimental and theoretical studies have attempted to improve nanotube yield 
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during arc synthesis, and to provide a comprehensive understanding of the growth 
mechanism. For the future, there is an important challenge for scaling up, optimizing, and 
controlling the arc process to make it commercially viable. To perform these objectives, 
several experimental approaches were developed to study the effect of the different arc 
parameters on nanotube yield and structure. These studies were also complemented by 
several models, either phenomenological, to understand how and where nanotubes are 
formed, or macroscopic to understand such processes as mass and heat transfer in the arc 
and their interconnection. When implemented into computational fluid dynamics CFD 
numerical codes, these models can predict temperature, species distribution, and fullerene 
yield or nanotube growth rate. In this section, we will summarize first the experimental 
approaches used to enhance fullerene or nanotube yield in the arc. Then, by a 
comparative survey of all synthesis processes, we will point out the specificity of the arc. 
Phenomenological models, as well as macroscopic multi-dimensional models at OD, ID, 
and 2D will be presented. These models were developed to predict the temperature and 
chemical species distributions, fullerene yield or nanotube growth rate. 

§.3. 1 Review of the Arc Process 

Several laboratory-scale methods have been proposed to produce SWNTs. 
Examples include condensation of laser l. 63 , 64 " 1 vapori zetPcarbon in the presence-oF •- ~ 
catalysts, decomposition of hydrocarbons on supported [54-56], entrained [ D ", 6L ', 67 , o8 j or 
fluidized bed catalysts, [59,60] decomposition of carbon monoxide on supported catalysts 
[58] and a high-pressure technique called HiPco. This technique involving decomposition 
of carbon monoxide on gas phase metallo-organic carbon pentacarbonyl Fe(CO )5 was 
developed more recently by Smalley’s group at Rice University [ 69 , 70 ]. Compared with 
these processes, the unique aspects of arc synthesis are the presence of the ions in the 
discharge and the formation of catalysts atoms and small carbon clusters C and C 2 at 
relatively high temperature. After giving a brief history of the arc, we will present in this 
section the uses of arc technique to produce carbon materials as fullerenes, multi-walled 
or single-walled carbon nanotubes. We will discuss in more detail the specific conditions 
of single-wall carbon nanotube formation. 

§.3.1.1 History 

Sir Humphry Davy of London, England discovered the electric arc between two 
carbon electrodes in 1813. In 1844, the French physicist Leon Foucault studied physical 
phenomenon in the arc lamp. The first industrial use of the arc was achieved in 1862 to 
produce acetylene by a process developed by Marcellin Berthelot. In 1880, Louis Clerc 
replaced the oxygen flame for welding applications by the plasma of an electric arc. At 
the end of the 19 lh century, carbon arc welding (CAW) was intensively used for 
locomotive maintenance because of the weld joints were hard and brittle due to the 
carbon flaking off into the weld puddle [ 71 ]. Since 1944, these discoveries opened the 
way to the thermal plasma synthesis of carbon materials. 

The German physicist, Otto Hahn, obtained carbon chains when trying to create 
some heavy atoms by fixing neutrons during the evaporation of metals in a carbon arc. 
Harold Kroto and Dave Walton also observed these chains in 1970 during experiments 
simulating conditions of combustion in the formation of red giant stars. In 1990, fullerene 
synthesis was performed by Kratschmer et al. [2] (Fig. §.9) of the Max Planck Institute in 
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Heidelberg, Germany under a helium atmosphere. One year later, in 1991, nanotubes 
were observed for the first time by S. Iijima [ 72 , 73 , 74 , 7;> ] of the Japanese company, NEC, 
and by Bethune et al. [ 76 ] Due to their unique combination of properties, nanotubes 
generated a lot of interest in the scientific community, and opened new fields of studies in 
science and technology. Since 1991, nanotube synthesis in the arc has been extensively 
developed and studied by several research teams around the world. 

§.3.1.2 Experimental approaches 

Before presenting the different modelling approaches, we will briefly discuss the 
experimental approaches used to improve fullerene and nanotube arc processes. 

§.3.1.3 Fullerenes vs. multi-walled carbon nanotubes 

When the electric current is changed from direct (DC) to alternating (AC) current, 
nanotubes instead of fullerenes are formed in a deposit on the cathode. Hence, the arc 
discharge technique, first used to produce fullerenes [2] was then extended to producing 
multi-walled carbon nanotube [2, 77 ] synthesis by changing the conditions of the discharge 
in the same apparatus. Fig. §.8 illustrates ideal structures obtained in the soot. Multi- 
walled carbon nanotubes are formed inside a hard deposit that grows on the cathode at a 
speed of about 1 mffl/mrn Thty ^re'Orgatrizedm^rmdles ' vith diameter ’ between^ and • 
30 nm [72,77, ]. Ebbesen ct al. [ ] and Ajayan and Iijima [77] proposed improvements 
in the classic process of Kratschmer and Huffmen. These studies showed that optimal 
growth conditions are different from of those of fullerenes. For example, the total 
pressure of helium is 660 mbars for nanotubes and only 200 mbars for fullerenes. The 
physical reason of this difference remains not very elucidated. Among the limited number 
of publications on this topic, we notice Zhang et al. [ 79 ] work on the effect of the helium 
pressure on the MWCNT yield. They indicated that the increase in pressure enhances 
nanotube yield and decreases fullerene formation. Cadek et al. [ 80 ] reported a pressure of 
660 mbar and a current density of 195 A/cm 2 as optimal conditions for MWCNT 
synthesis. Nanotube synthesis under an argon atmosphere was studied by Borisenko et al. 
[ 8I ] who showed that argon can improve MWCNT yield by a factor of 5-10 compared to 
helium in similar arc discharge conditions. A hydrogen discharge was also tested by 
Ando et al. [ 82 ] and Zhao et al. [ 8j ] to grow nanotubes with smaller diameters. Double 
walled carbon nanotubes were produced using a mixture of argon and hydrogen by 
Hutchison et al. [ 84 ] Finally, Jung et al. [ 8:> ] Tang, et al. [old 42 not Tang] ???studied the 
morphology and the structure of MWCNT produced in a mixed atmosphere of helium 
and acetylene in the arc. They noticed that the quantity of carbon nanoparticles, formed 
inside the cathodic deposit, increases when acetylene is added to the ambient gas. An 
important modification of the arc method was achieved using liquid nitrogen [85, 86 ] with 
the advantage of an easier recuperation of the products and a possible continuous 
synthesis. Another less efficient, but also less expensive method, uses water instead 
liquid nitrogen [ ]. 

§.3.1.4 Single-walled carbon nanotubes 

Single-walled carbon nanotubes were discovered by Iijima [71] in 1991 when adding iron 
to a graphite anode in the arc. Later, all the process parameters such as catalysts and inert 
gas composition, current density, pressure, gravity, etc. were varied to enhance nanotube 
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yield. Several metallic particles were used as catalysts including cobalt [76], nickel [ ], 
yttrium [ 90 ], manganese [ 91 ], scandium [ 92 ], lanthanum [92, 9j> ] vanadium [ 94 ], cerium [ 9; ], 
gadolinium [ 9<5 ], 52 and zirconium [90 49 ]. The presence of these catalysts leads to a large 
variety of metallic particles and carbides in the plasma and in the collected material. 
Currently, the most efficient catalyst is obtained by mixing two elements such as nickel- 
yttrium [46, 971 , nickel-cobalt [ 98 . 99 ], iron-nickel, [ 100 , 101 ] or rhubidium-platinium [ 102 ]. 
According to Saito et al. [102], the efficiency of the Rh-Pt couple is comparable to Ni-Y. 

It is also noticed that these catalysts can change the diameter distributions in the interval 
0.7-3 nm. On the other hand, Sugai et al. [ Kb ] used a pulsed discharge to synthesize 
fullerenes and SWNTs using nickel and cobalt catalysts under an inert atmosphere of 
helium, argon, or krypton. Argon was found to be the best gas for fullerene production 
and krypton the most efficient for nanotubes. According to the authors, increasing the 
pulse duration increases nanotube and decreases fullerene yields, making a competition 
between these two kinetic pathways. The conventional electric arc is an unsteady process 
because of instabilities of cathodic spot. This leads to a non-homogeneous distribution of 
the electric field in the plasma. The cathode spot instability results in random and erratic 
motion of a luminous spot along the cathode surface. Hence, Lee et al. [ I04 ] used a 
technique of rotation of the anode to create uniform plasma. This rotation generated a 
t«FbulcFiGe tbat accek v rates carbon-snee^ s.-p&rpen4i6u!arly to the anode fivihis situation; - 
carbdn is not condensed on the surface of the cathode but collected on a graphite 
collector placed in the periphery of the plasma. According to the authors, increasing the 
speed of rotation increases nanotube yield and decreases nanotube diameter. This result 
was confirmed by Bae and al. [ 1(b ] Effects of gravity have been studied by Kanai and al. 
[ 106 ] who found that the absence of gravity forces can reduce convective fluxes of the 
inert gas increasing plasma volume. Consequently, the authors found an increase of the 
yield of nanotubes and an increase of their diameter. A key parameter for the arc is the 
erosion rate of the anode. This parameter depends on the input power, the chemical 
composition of the anode, the length of the arc, and the nature of the inert gas, the 
pressure, the cooling of electrodes as well as the geometry of the reactor. Not all these 
parameters are independent, increasing the complexity of the problem. Therefore, Zhang 
et al. [ 107 ] studied the dependence of the speed of erosion of the anode with the pressure 
and the nature of the inert gas (argon or helium). The authors found that for low pressures 
between 100 and 300 mbars, the erosion rate of a pure graphite anode is lower in argon 
than in helium. This tendency is reversed for pressures between 300 and 900 mbars. 
Takizawa et al. [ I08 ] measured the quantity of soot formed in the reactor while varying the 
percentage of nickel and yttrium in the composition of the anode. They demonstrated that 
the addition of yttrium traces (0,1% at.). with the nickel increases soot yield by a factor of 
2-3. This yield was not important when only nickel or yttrium has been used as catalyst. 

§.3.2 Analysis of the Arc Discharge 

§.3.2.1 Description of the experiment 

The carbon arc apparatus consists of a static water-cooled reaction chamber with 
two cooled graphite rods. For the anode, a graphite rod filled with catalysts (Co, Fe, Ni, 

Y) is generally used. The cathode is made of pure graphite or copper. It was reported that 
copper does not improve the results and leads to a less stable arc [ 109 ] 63 . The distance 
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between the anode and the cathode can be adjusted by moving the anode toward the 
cathode manually [46,52] or automatically [49] using an optoelectronic system. A direct 
current (DC) for nanotube synthesis or alternating current (AC) for fullerene production 
passes through the electrodes; and plasma is created in the inter-electrode region. 
Efficient operation is assumed to exist when the discharge is stable and the anode erosion 
rate is constant. This can be achieved by maintaining a constant voltage between the 
electrodes, which is closely related to the stability of the electrode spacing. The plasma is 
first ignited by contact between the anode and the cathode, which elevates the 
temperature of the contact point until evaporation of the anode material. Then, the anode 
is moved back to maintain a desired gap between the growing deposit on the cathode and 
the burning anode. An active plasma zone bounded by the deposit and the anode is 
created. This hot plasma zone produces carbon and catalyst vapors, which then diffuse to 
the cooled reactor regions. Carbon species and catalysts build up on the end of the 
cathode (called cathode deposit) and deposit forms a collaret around the cathode deposit. 
Soot is deposited on the reactor walls by free convection. The high temperature near the 
anode and the high energy density in the plasma insure vaporization of most of the anode 
material. The water-cooled cathode leads to high quench rates and high levels of super 
cooled or supersaturated vapor with nanotube formation. The quench process is 
unceMF&lled, but-arc- products usually -aie’.-sooi- ©a -th®- reactor walls, web-like •structures 
between the cathode and the chamber walls, a hard gray deposit at the cathode’s end, and 
a rubber-like collaret rich in nanotubes around this deposit. Fig. §.10 Figure 4 shows 
SEM and HRTEM micrographs of as-produced nanotubes. 

§.3.2.2 Comparison with the high temperature processes 

The nanotubes produced by arc discharge are similar to those obtained at high 
temperature process as laser or solar vaporization; and it is interesting to compare all 
these processes to point out the specificity of the arc. In high temperature processes, the 
temperature is higher than 3500 K, permitting the total vaporization of the target or anode 
formed by graphite and catalyst. To achieve such high temperatures, high energy 
densities are needed to totally atomize the solid and to form atoms of carbon C and 
catalyst M. When the temperature decreases in the cold region of the reactor carbon 
reacts to form molecules of superior size and clusters composed only of carbon, such as 
C 2 , C 3 ... C6o-.- CNT..., soot or composed of only metal as M 2 , M 3 , M C | US ter or, finally, 
a mixture of them C n M m . To understand these processes, it is necessary to quantify mass 
and energy fluxes at different conditions. In Table §.131, a comparison of fluxes of mass 
and energy for the most important high temperature processes is given. In this table, I and 
V are respectively the current and potential drop between electrodes, r A is anode radius, 
d A c is the distance between anode and cathode, P the pressure, X is the laser wave length, 
v the frequency of the laser, E pu i S e is the pulse energy, t pu i S e is the duration between two 
pulses and <t> is the laser or solar spot diameter. We notice that the arc process is by far 
the most efficient process in terms of erosion rate and produces 4.2 g/h of soot containing 
nanotubes. Nevertheless, for industrial applications, arc processes have at least two 
disadvantages. (1) This process is discontinuous and requires cycles of 
production/cleaning. (2) The nanotube concentration in the soot is relatively poor. 
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§.3.3 Specificity of the arc 


§.3.3.1 Specific conditions for SWNT growth 

Typical experimental conditions [ no ] for SWNT synthesis in the arc are the 
following: 

(i) The inert gas is helium. 

(ii) The potential drop between the electrodes is V= 40 V. 

(iii) The anode diameter r A =3 mm and the cathode diameter rc= 16 mm. 

(iv) The current density j = 100-350 A/cm 2 . 

(v) The interelectrode distance d A c = 1-3 mm 

(vi) The plasma temperature ~ 6500 K. 

(vii) The erosion rate 70 g/h. 

(viii) The pressure 660 mbar. 

(ix) The velocity of material leaving the anode as calculated by Farhat, et al. 
[107 661 u an0 de=78 m/s and estimated by Krestinin and Moravskii [7, 1 1 1 ], 
u a node=2-80 m/s for fullerene conditions. 

(x) The dilution factor x describing how much time the carbon material leaving 
the interelectrode space is mixed with helium. x= 20-30. [7,1 1 1] 

fx/H .^lxmgdime f n’a^ i s a'p. h ftnel or jr«U,ouucc^ By Krestinin and Moravski L' 
[7,1 1 1] to measure the time that carbon material needs to be totally mixed 
with helium when leaving the interelectrode space: x m j x =2-8 ms [7,1 1 1], 


Other parameters such as water flow-rate used to cool the cathode and loss by radiation 
from the plasma are reported in the literature [ 112 , 113 ]. Gamaly and Ebbesen [ 114 ] proposed 
to estimate gas radiation simply from Stephan-Maxwell equation Qrad = ctT 4 =5.67xl0" 12 
T 4 J/cm 2 s K 4 with T~ 4000 K. However, arcs are not usually optically thick. Therefore, 
this is not a good assumption. 

§.3.3.2 Space charge, potential and electric field distribution 

Space charge, potential, and electric field distributions have been described by 
Gamaly and Ebbesen [114] for multi-walled carbon nanotube conditions. According to 
them, the space distribution of the potential has a steep drop near the cathode in a region 
of positive space charge where practically the entire potential drop occurs as showed in 
Fig. §.11. They calculated this sheath length at A=12 pm for typical multi-walled carbon 
nanotube conditions with V=20 V and j= 1 50 A/cm 2 . This is also the distance of the main 
potential drop in the interelectrode region. Gamaly and Ebbesen [114] also estimated the 
average electric field in this region at E=2xl0 4 V/cm and consider that the electric field 
in the outer region is several order of magnitude lower. 

§.3.3.3 Ionisation state and density of the plasma 

The interelectrode region is composed of a mixture of inert gas, carbon and 
catalyst atoms and molecules. Fig. §.12 shows qualitatively 15 pictures of the 
development of the arc in the argon as the distance between the electrodes is changed. 
Light is emitted from atom and ion lines, as well as from C 2 and, possibly, CN molecules. 
The first ionisation potential of some elements present in the plasma are summarized in 
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Table §.14. For a helium/carbon discharge without catalysts used for MWCNT synthesis, 
Gamaly and Ebbesen [1 14] estimated that the carbon ion C + is the major ion, with the 
density in the plasma of nc+ = SxlO 13 cm'' and near the cathode of no = 6.9xl0'~ cm°. 
The ion velocity was estimated at Vj on = 1.8X10 5 cm/s. At these conditions, the flux of 
carbon to the cathode is 6.25xlO iy carbon/s and the number density of the plasma is 
dominated by helium with nn e = 6.4xlO ls cm ’. Nevertheless, Scott et al. found that 
yttrium and nickel ions are the major ions in the discharge when catalyst is added to the 
anode for SWNT optimal conditions. Figs. §.13 and §.14 show measured optical 
emission spectra and that calculated with the atomic emission spectroscopy Fortran 
program named AES and developed by J. Homkohl of the Tennessee Space Institute. 
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§.3.4 Modeling Arc Process 


After presenting the main results obtained by the analysis of the experimental arc 
process, we need to develop a systemic approach in order to establish a relation between 
nanotubes or fullerene yield and arc parameters. Before establishing such relations, it is 
legitimate to ask some questions. 

• Where and how nanotubes are formed in the reactor and how can this 

perfect linear and symmetric structure be formed in such an anisotropic 




plasma7 


• Why in the same reactor, but with different conditions we obtain 
fullerenes C 60 and C 70 , MWCNT or SWNT? 

• What are the growth precursors, and by which kinetic mechanisms are 
they formed and transported in the plasma? What is the role of the ions 
in the growth of nanotubes? 

• What process parameters control this growth and how do they influence 
nanotube diameter and chirality? 


Due to the complexity of time and space phenomena in the arc, and the 
interconnection of several individual processes such as electric, mass, and energy 
transfer, answering all these questions is a difficult task without the help of modeling. 
Very few works are reported in the literature to determine arc characteristics in the 
specific conditions of fullerene or nanotube growth. We selected some approaches, 
developed mainly by Gamaly and Ebbesen [1 14 77 ] for multi-walled carbon nanotubes, 
Bilodeau et al. [ u5 ] and Kristinin and Moravskii [7,108] for fullerenes, Hinkov, [ u6 ] 
Scott, [ ] and Farhat et al. [110] for single wall carbon nanotubes. 

§.3.4.1 Phenomenological growth models 

A quantitative description of nanotubes was first described by Gamaly and 
Ebbesen [114] based on the approximation of the distribution of density, velocity and 
temperature of carbon vapors, electric charge, potential, and electric field in the arc. On 
this basis, the authors presented a scenario for the nanotube formation, growth, and 
termination in time and space scales. It is doubtful that some of the details of this model 
are correct. However, it points out some useful facts. Typical multi-walled nanotubes 
conditions reported are obtained with the vaporization of pure graphite anode in a helium 
discharge chamber with P=500 Torr, anode to cathode distance d A c<lmm, a potential 
drop of V=20 V and a current density of j=l 50 A/cm . At these conditions, average 
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plasma temperature was assumed to be constant and ~ 4000 K. The deposit of area 0.5 
cm 2 containing MWCNT grow at the cathode with a measured rate of 1 mm/min =16 
[xm/s. Assuming the average material density of 1.5 g/cm 3 , this corresponds to a flux of 
carbon of 1.3xl0 20 atoms/cm 2 /s or a flow rate of 6.25xl0 19 atoms/s. At plasma 
temperature of 4000 K, calculated C + density in the plasma is nc+ = 5x10 15 cm 0 and near 
the cathode nc+ = 6.9xl0 13 cm°. The ion velocity was estimated at Vj on = l.SxlO 3 cm/s 

1 O O 

and the plasma is dominated by helium with nHe = 6.4x10 cm"' . Since the temperature 
of the plasma is close to the melting (4100 K) and boiling (4470 K) temperatures of 
graphite, carbon begin to evaporate producing a thin layer of saturated carbon vapor near 
the cathode surface. By a kinetic theory calculation of the evaporation, the authors 

i o n 

determine that the carbon vapor sheath density is n vap ~ 1.35x10 cm" which is very 
close to the density of surrounding gas. The carbon vapor initially expands from the gas 
and forms a thin shield over the cathode surface. There the gas is cooled in a region very 
close to the cathode of only 2 to 3 pm in thickness having a temperature gradient of 26 
K/cm. The saturated carbon is in fact diluted with helium in such way that the density 
after mixing is n vap ~ 4.8xl0 17 cm"' 1 . According to the authors, the growth is due to the 
competitive input of two groups of carbon having different velocity distributions. The 
first group comes from the anode and has a Maxwellian velocity distribution. The second 
group is composed of ions accelerated in the gap between the positive' space. charge-anri - 
the cathode. Due to the large difference in the first ionization potential between carbon 
(1 1.26 eV) and helium (24.87 eV), the major ion in the plasma is C + . At a potential drop 
of V=20V, the ion velocity, calculated from the kinetic energy, equals 2xl0 6 cm/s, which 
is ten times higher than thermal velocity. The presented scenario of deposit growth 
consists of cycles of four events. 

1) Seed structure formation during the establishment of a steady ion current by collisions 
between C* + C* leading to C 2 which is the feedstock for the seed structure growth. The 
small characteristic time of C 2 formation (1C)" 7 - 10" 8 s) permits forming about 10 7 seeds 
in a micrometer-thick layer. 

2) TlMtTHbaUoubefgDdTO tirdvvTtiiiKttiQestariTdaHm^iDftdinlaliiadTQaigecdljTceitbt ranKlariiffree or 
ppdn tAeflrisointearbpti arm liraJawsttjiSi ngrpfi sHaidjsc tagahei tlkmthitincRnaBi (Bah gsra in 
eartgrnaii by the authors is due 
to the splitting of the current into thinner filaments threads that moves randomly to the 
cathode. 

4) The tube ends are capped by a process involving the rearrangement of carbon having 
Maxwellian distribution in the absence of current. 

Based on this model, and the assumption that the characteristic reaction time for 
carbon-carbon attachment is proportional to the carbon-carbon collision time, the authors 
estimated the time of creating a 1 pm length tube with a diameter of 5 nm at -0.02 s. 
Their conclusion is based on the qualitative agreement between the predicted time of 
these two groups with the experimental time -0.07 s. Gamaly and Ebbesen pointed out 
the importance of cooling the cathode in reducing the velocity of carbons with 
Maxwellian distribution, thus reducing the perturbation amplitude and improving yield 
and quality of nanotubes. They finally concluded that the extremely high local field of 

o lift 

-10 V/cm as suggested by Smalley [ ] to explain open end nanotube growth is four 

orders of magnitude higher than their calculated average electric field of 20 V / 10 mm 
~2xl 0 4 V/cm. 
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§.3.4.2 Plasma models 

The majority of the work of modeling electric arc discharge in presence of carbon 
and inert gas assumed local thermodynamic equilibrium (LTE). [110, 7, 111, 114, 115, ll9 ] 
This assumption consists of considering a unique temperature to represent the plasma. 

Thus the temperature of electrons T e , ions Tj and atoms and neutral molecules of the gas 
T g are equal. This assumption has been justified by the work of Bilodeau et al. [115] 
concerning fullerene modeling in the arc. According to them, pure helium plasmas show 
departures from LTE at pressures of 10 4 Pa. However, their calculations demonstrate that 
carbon species dominate in the arc region, resulting in an increase of the electric 
conductivity in comparison to pure helium. As a result, the electric field is limited to 
values of 10 V/cm, thus limiting deviations from LTE. In order to verify this assumption, 
Hinkov [116] calculated features of a plasma used at specific conditions of SWNT 

A . . 1 90 

growth using the AURORA plasma software in the CHEMKIN collection [ ]. This 

software predicts the steady state or time averaged properties of a well mixed or perfectly 
stirred plasma reactor. By modeling non-thermal plasmas, AURORA determines the ion, 
electron, and the neutral radical species concentrations, and the electron temperature. 

Well stirred reactors are characterized by a reactor volume, residence time or mass flow 
rate, heat loss or gas temperature, surface area, surface temperature, the incoming 
temperature ami mixture composition . as well a s the newer deposited into the plasma for - 
non-thermal systems. Fig. §.15 Figure 9 is a schematic representation of the well- mixed 
plasma reactor [120] and its adaptation to predict plasma compositions [116]. The reactor 
of volume V=0.6 cm'’ and containing a reactive cathode surface area of A m =2 cm 2 is 
continuously fed by a mass flow of 1800 seem coming from the anode erosion. The 
composition of the inlet flow is calculated at 4000 K using Ivtanthermo ['"] equilibrium 
computer code with the initial composition of the anode C:Ni:Y 94.8:4.2:1 at.%, diluted 
by a factor of 20 in the helium, as suggested by Krestinin and Moravskii [11 1]. The 
model takes into account 92 species including carbon radicals C, C2, C3, .... C79, 
fullerenes C6 of and C 70 f, inert gas He and ions Ni + , Y + , He + , C + , C 4 + , C 5 + and C6o + . All 
these species involves 562 electron and neutral reactions listed respectively in Tables 
§.15 and §.1. In addition to the homogeneous gas chemistry, the surface reactions listed 
in Table §.16 5, were added to account for the growth and etching of carbon nanotubes 
assumed to grow at the cathode surface from small carbon atoms. The sticking coefficient 
of unity for the deposition and 0.1 for etching was arbitrarily chosen. Finally, ions 
recombine in the sheath by the application of the Bohm condition for ion fluxes to 
surfaces. Indeed, as suggested by Meeks et al. [120], it is reasonable to constrain the ion 
flux to a surface according to the Bohm criterion. This condition results in the maximum 
net flux of a particular ion to a surface equal to the product of the ion density and the 
Bohm velocity. No data are available for recombination of nickel and yttrium ions in the 
sheath, hence an arbitrarily Bohm criterion of BOHM=0.4 was chosen by Hinkov [116]. 
The calculations were performed with the following conditions, inert gas=helium, 
measured erosion rate G er osion= 17,2 mg/s at 100 A, assumed dilution factor Xdiiution= 20 
and a total pressure P= 660 mbar, with and without the LTE assumption. In the first 
calculation, the LTE assumption with T g =Tj on =T e =6500 K was assumed. The output 
results from AURORA are summarized in Table §.17 and reveal a residence time of 
about 28 ps and a weak electronic density of 6.4xl0 14 cra°. The predicted total ion 
current density of 10.4 A/cm 2 corresponds to a total current of the order of 20 A, based 
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on a cathode surface area of 2.01 cm 2 . This is much lower than 100 A, making the LTE 
assumption questionable. The second calculation was performed in the same conditions 
by maintaining gas temperature at T=6500 K and varying electron temperature from 0.56 
eV (6500 K) to 5.17 eV (60000 K) corresponding to varying the current from 20 to 200 
A. The estimated ion mass fractions, plotted in Fig. §.16, reveal that the major ion is Y + 
for the current intensity I <50 A and Ni + for the current intensity I >50. This result is 
explained by the low quantity of yttrium (1% at.) added to the anode material, and to the 
weak energy of ionization of the yttrium. Increasing the current intensity, one forms more 
and more Ni + ions from nickel atoms. Fig. §.17 gives the evolution of calculated current 
density and electron temperature with electron density. In the optimal domain of current 
(80-120 A) for nanotube synthesis, the calculated electron density is between 1 .25x10' 3 
and 1.5xl0 13 cm' 3 , and the electron temperature is 1.7 - 2.4 eV. Table §.18 summarizes 
the plasma composition calculated for 1=100 A where the electron density of n e - 
1.4xlO l3 cm° is in good agreement with the measurements of Akita, et al. [51] in the 
interval 10 13 - 10 16 cm' 3 for T= 6000 K as measured by optical emission spectroscopy in 
the following nanotube conditions: inert gas = helium, pressure P = 300 Torr, G er osion = 50 
mg/s, and the flow rate of carbon without catalysts is 5500 seem. The calculations 
indicate that when the catalyst is present the major ion is Ni + with a molar fraction of 

Howed by Y A withamo-r fraction o' 7 ; 7 ||V ' mdC* with a molar fraction of 
1.3x10 . From these results, even if the jlTE condition is not satisfied, the weak rale of 
ionization of the carbon (<0.3%) indicate clearly that C + ions can not explain the fast 
nanotube growth rates, as was implied by the model of Gamaly and Ebbesen, described in 
Section §.3.4.1. 

§.3.4.3 Zero-dimensional models 

A zero-dimensional model for arc nanotubes synthesis was developed by Farhat et al. 
[110] and Hinkov [116] to calculate the equilibrium chemical composition in one point in 
the plasma, assuming LTE and given the gas temperature and pressure. The objective is 
to study the sensitivity of the kinetic model developed by Krestinin and Moravskii [119] 
in the specific conditions of nanotube growth. 

§.3.4.3.a Model Equations 

The multidimensional conservation model is reduced to zero spatial dimensions 
(0D) by assuming steady state, and a given temperature and zero velocity. The 
conservation equations can be written as a balance between an accumulation term on the 
left hand side equation (§.27) and chemical source term. 

dY, M f 

— ~ 0=2,..., ns) (§-27) 

dt p 

In this equation, n s the number of gas species, Yj is the mass fraction, 0)j is the chemical 
molar production of species i (mol. cm' 3 . s' ! ), Mj is the molecular weight (g. mol' 1 ), 
respectively of species I, p is the total mass density (g. cm' 3 ), and t is the time (s). 

The source term 0)j is calculated by considering R elementary reversible or irreversible 
reactions involving n s reacting chemical species A;. 

ns j ns 

< r=1 *> «•»> 

/=i 2 ,=i 
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where v’; r and v”j r are the stoichiometric mole numbers of the reactants and products 
respectively and Aj the chemical symbol for the i th species. The chemical production rate 
(Oj of i lh species can be written as a summation of the rate-of-progress variables for all the 
reactions involving the i th species : 

n r 

^•=S Al/ ir^r (r = l,...,n r ) (§.29) 

r=l 

where AVj r = Vj r " - Vj r '. The rate-of-progress variables q r for the r th reaction is given by the 
difference of the forward rates 1 and the reverse rates 2 as: 

ns ns 

^Kn^-KYlcr* (r = 1. •">*) C§-30) 

i=l 1=1 

Where ki r and k 2 r are the forward and reverse rate coefficients of the r th reactions and C; 
is the molar concentration of i lh species. Reverse reactions are written explicitly in the 
forward sense. The forward rate constants are calculated for each reaction r by assuming 
Arrhenius temperature dependence: 

k lr =A r T fi 'cx pP^-) (r = 1, ...,R) (§.31) 

Kl 

The' rate AtFprogress variables^ are ctrmputed'using CHEMKIN database and software-' 

package. 

The condition of total mass fraction given by. 

ns 

(§32) 

(=i 


§.3.4.3.b Kinetic and thermodynamic data 

The gas phase chemistry used in this approach involves neutral carbon species 
and includes small clusters (C) - Cio), cycles and polycycles (Q ] - C 31 ), fullerene shells 
(C 32 -C 46 ), fullerene clusters (C 47 -C 79 ) and two fullerene molecules C 6 of and C 7 of- In 
addition, argon and/or helium, as well as nickel and yttrium, are considered as inert 
chemical species. There are at least 554 chemical reactions, listed in Table §.l 4, that 
describe the chemical kinetics of carbon vapor condensation and clusters of different 
sizes in the arc. The source of these data and the formation enthalpies are discussed by 
Krestinin and Moravsky [119], The symbols <-> and — > in each reaction indicates that this 
reaction is reversible or irreversible respectively. Finally, this model includes 81 carbon 
species (Ci, C 2 ..., C 79 , C 6 of, C 70 F) involving all the reactions of Table §.l. The 
thermodynamic properties of species are fitted in the temperature range of 300 to 20 00 0 
K from Ivtanthermo, [122] JANAV [ I23 ], NASA [ l24 ], and CHEMKIN [40] databases and 
given in CHEMKIN (old NASA) format. This format presumes that the standard-state 
thermodynamic properties as Standard state molar heat capacity at constant pressure of 
the ith species c P i, standard molar enthalpy H°, and standard molar enthalpy 5° are 
thermally “perfect”, in that they are only functions of temperature and are given in terms 
of polynomial fits (§.2)-(§.4). For small carbon species, argon, helium, nickel and yttrium 
the polynomial coefficients are listed in Table §.4. 
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§.3.4.3.c Solution procedure 

The solution of differential-algebraic-equations (DAE) representing mass 
conservation equations (§.27) and (§.32) employs a backward differential formula 
method implemented in LSODI package [ 125 ] . This software was designed to solve stiff, 
DAE and the solution procedure attempts time-integration from an initial guess of the 
variables. The initial composition of the gas was estimated at a given temperature and 
pressure using IVTANTHERMO [121] equilibrium computer code, with the initial 
composition of the anode C:Ni:Y 94.8:4.2: 1 atom. %, diluted by a factor of 20 in helium. 
Time-stepping improves the starting point by relaxing it closer to the steady state 
solution. This solution method is much faster and much more stable in reaching steady- 
state convergence than the classical Newton iteration method. The model takes into 
account 84 equations corresponding to all the species including: radicals C, C 2 , C 3 , ..., 

C79, fullerenes C 6 of and C 70 F, nickel Ni, yttrium Y and inert gas He. These species involve 
the 554 reactions of Table §.l. The typical time evolution of the solution for carbon atom 
mole fractions calculated at T=4000 K and P=660 mbar, is given in Fig. §.18. The 
characteristic reaction time to obtain steady-state C mole fraction is about lp.s. At a 
cooler temperature, T=2000K, we start to form fullerenes from carbon clusters. Figs. §.19 
and §.20 show the time evolution of carbon clusters Csg,C59, C(,q and the fullerenes C6of 
and C70F. respective! v. Fullerenes have a characteristic reaction time of about' 10 us — 
Finally, in Figs. §.2 la and b are plotted mole fractions of small carbon species and total 
fullerene yield versus temperature calculated in helium at 660 mbar. We can see from 
these figures that fullerenes start to be formed at 3000 K. 

§.3.4.4 Turbulent Fan Jet Model 

When the electrode spacing is very small, carbon is ejected from between the 
electrodes as a radially expanding jet. Krestinin and Moravsky [7,1 1 1,1 19] developed a 
model to describe the chemical processes taking place in a turbulent fan jet leaving the 
interelectrode space. 

§.3.4.4.a Model formulation 

The turbulent jet model accounts for the main processes controlling fullerene 
formation in electric arc, namely, 1 ) cooling and mixing of carbon vapor in a buffer gas, 

2) reactions of cluster growth and decomposition under non-isothermal conditions, 3) 
formation of soot particles and heterogeneous reactions on their surface. In this model, 
the flow from the arc zone is idealized as an axially symmetric fan jet represented in Fig. 
§,324:&fc> the reactor, the temperature gradient in the gap, as well as 

the rateTfidi^ft^'titi0Hsiqaaf1iM8'4it<h^iU©'ne>l^hi(wk bdUflda^a^leclddtuFhulrmtdel 
sxjisgstffafefgM d wafctcueqosti 1 ttasjgo iftm ftsrrnt 0 f an initial volume AV 0 at the 
temperature To with an initial high velocity Uo of a jet surrounded by a pure and cold 
inert gas at T at . The velocity profit U in tl^plane (z,r) representing a jet section is shown 
schematically in Fig. = ^ k ^ P ^ +M ^ ( i=l >-’ n *) (§' 33 ) 

pC ' u fr = ^ rPC 'lr +/r (§34) 
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where e T is the coefficient of turbulent momentum transfer, k c and kj are empirical 
coefficients relating momentum transfer to mass and heat transfer, 7, and T respectively 
the mass fractions of the species i and the temperature at the position longitudinal, r, and 
transverse, z, coordinates of the jet. The terms, /• and/ 7 - are the rates of mixture 
component formation and the rate of heat release due to chemical reactions, Mj are the 
molecular weights of the mixture components, p is the gas density, C p is the mixture 
specific heat, and U is the mass-mean gas velocity. To solve these partial differential 
equations, the semi empirical theory of free turbulent jet developed by Abramovich 
1 126 1 is applied. In this theory, the ratio between the half-width of a free turbulent jet b 
and the mixing length L is assumed constant at each longitudinal position r of the jet. L/b 
= const ~ 0.2. In addition, due to the jet symmetry, derivatives of all the variables with 
respect to z are equal zero at the jet center plane. Hence temperature and concentration 
vary little within the turbulent mixing length L (r) and can be assumed constant and equal 
to their values at z=0. Thus, integrating the equations (§.33) and (§.34) with respect to z 
over interval [-L,L] yields to an approximate description of evolution of an elementary 
gas volume moving near the symmetry plane of a fan jet. By replacing the derivatives 
with respect to z by finite differences: 


(§■35) 


c> T I Yy-Y, 

— b 

partial differential equations (§33) and (§34) are transformed via a procedure described 
in detail in [7] to a set of ordinary differential equations. These equations are written in 
CHEMKIN notation as following, 


at p 




(§36 13) 


£ K c r , ~ = £ r, C PJ k T (r. -T)+o>,h,^- 

CM ,=1 P 


(§37 14) 

1=1 Ul 1=1 P 

where Y, and Yo,i are the mass fractions of mixture components including carbon black 
particles in the reaction volume and buffer gas respectively, T is the temperature, 0 )j is the 
chemical molar production of species i (mol. cm' . s' ), Mi is the molecular weight (g. 
mol 1 ), p is the total mass density (g. cm'’), hj is the specific enthalpy of the i lh species 
(erg/g), C p> i is the specific heat at constant pressure of the i lh species (ergs/g/K)) and t is 
the time (s). The mass and heat transfer coefficients k c and k T are calculated by the 
empirical formulas: 

*.=- 032 " (§38) 


K = 


T ■ + 2 1 

mix 

0,64 


T +2t 

mix 


(§39) 


were x^x is the time of turbulent mixing calculated from the flow velocity in the jet 
within the arc zone Uq using equation (§.40). 
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(§•40) 
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The flow velocity Uois calculated by equation (§.41) assuming that the carbon vapor in 
the arc zone is equilibrated with electrode graphite and taking the experimental value of 
the formation rate of fullerene carbon black, V soot 


= 


2tt r A d AC P 


(§•41) 


where dAC is the anode to cathode distance, Ta is the anode radius, R is the universal gas 
constant, P is the helium pressure in the reactor, and To the arc temperature calculated 
from the equilibrium pressure of the carbon cluster vapor over graphite. V soo t is 
determined from the measured rate of ablation of the anode 


§. 3.4.4. c Numerical results 

The representative experimental conditions of Krestinin and Moravskii are 
P=100-760 Torr, To=3600-3900 K and Uo=2-80 m/s. Under these experimental 
conditions, the model equations (§.36) and (§.37) were solved using the chemical gas 
model in Table §.l and Table §.19. Fig. §.23 compare calculated and experimental 
fullerene yields as obtained by Krestinin and Moravsky [8] at different pressures. 


§.3.4.5 One-dimensional models 

In order to analyze the effect of arc process parameters on temperature and" • «-• 
species profiles in (ID), Farhat et al. [Ill] used the highly structured computer packages 
SPIN and SURFACE CFEMKIN developed by Kee et al. [40] at SANDIA for rotating 
disc CVD processes. The complex gas-phase chemical reaction mechanism representing 
carbon condensation from the arc was added into numerical simulations in a one- 
dimensional (ID) model. This model was formulated under specific conditions of 
nanotube growth in the inter-electrode region. It solves for species, temperature, and 
velocity profiles in a steady-state, one dimensional stagnation-point flow, with 
temperature dependant fluid properties. 

§.3. 4. 5. a Model formulation 

This model is a boundary value problem consisting of a set of ordinary 
differential equations, solved by a finite difference procedure. It assumes local thermal 
equilibrium (LTE) and solves the steady-state axial and radial momentum, species, and 
energy equations in one spatial dimension between the anode and the cathode. The model 
accounts for carbon deposition at the cathode by a set of surface reactions that simulates 
nanotube growth. The steady state assumption is justified by the continuous adjustment 
of the inter-electrode gap leading to a constant erosion rate of the anode, hence a constant 
condensation of carbon vapor close to the cathode. Local thermodynamic equilibrium 
was assumed based on the Bilodeau et al. [115] model for fullerene synthesis by arc 
discharge in the same range of pressure as nanotube synthesis. 

§.3.4.5.b Model equations 

The governing equations that are solved in the SPIN code include continuity, 
radial momentum, species conservation, and thermal energy, and are given below. 


Continuity 


1 dp du u dp_ =Q 


p dt dx 


p dx 


(§•42) 


Radial momentum 


dV d dV dV 2 1 dp 

p — = - ( M )-pu- pV 1 — -f = 0 

dt dx dx dx r dr 


(§•43) 


Species continuity 


dY i d{pY'V,) dY, l£ 
y dt dx ^ dx 
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Thermal energy 

dT d 

P° p 


dt d x 


A 


dT 


d: 


- pc r u - jr % < ' Cr , pY , v , jp‘ Sfh ‘ ) * s ‘' {x) ~ Q ' M=0 
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Although these equations are stated in their transient form, the solution is obtained when 
aH'time derrvattres'are'/.eterfrTv ne mdep< •-•««« Wes~aw*the> — 

distance normal to the cathode x (Fig. §.24), and the time t. The dependent variables are 
axial u and radial V velocities, gas temperature T, gas-phase species mass fractions Yj. 

The mass density is given by p and the specific heats at constant pressure by c p . In the 
radial momentum equation, p is the spatially varying component of the pressure [40]. The 
molecular weight and specific enthalpy of species i are given by Mj and hi respectively. 

The viscosity and thermal conductivity are p and A. The net chemical production rate of 
species i by gas-phase reaction is 0)i. The species diffusion velocity Vj is calculated from 
mixture diffusion coefficient and species gradient [40]. The source term S q (x) in the 
energy equation accounts for the electrical energy dissipated in the arc. It was assumed 
that it is distributed in the form of Gaussian centered at x s with a (half) half-width of w s 
by the equation: 





-3(jc-xj) 2 

[T 
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j— e 
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Here, q accounts for the total power integrated over its full spatial extent 

r +w s 

q = S (x)dx and includes 100% of the net power added to the arc. The center of the 

J-w. ^ 


interelectrode gap was chosen as the peak in the distribution x s , and w s chosen to adjust 
the flatness of the distribution. In the present calculations w s was 0.15 cm, resulting in a 
very uniform distribution. The total energy q was calculated from measured electric 
power dissipated in the arc. It was corrected by the loss from arc plasma due to the gas 
convection. The radiative term Q ra d in the energy equation accounts for the net loss of 
energy by gas radiation. It was estimated from a curve fit of Owano [ 
argon plasma by: 


1 27 ] 93 in atmospheric 


^ mu .-141170 W 

Qrad =1.065x10 crexp( ) (— ) 

T m 
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with T the gas temperature in K and a a non -equilibrium factor accounting for the 
deviation from the equilibrium; a=l for local thermal equilibrium (LTE). Gas radiation is 
negligible when helium is used. Gas kinetics, thermodynamic and transport properties of 
species are estimated in the same way as in section §.3.4.3 

§.3.4.5.c Surface chemistry 

Carbon nanotube growth is considered as a boundary condition at the cathode by 
a set of surface reactions simulating nanotube growth. As a first approach, the open-end 
tube growth (scooter) mechanism developed by Smalley [118] and illustrated in Fig. §.25 
by an African calabash structure decorated by hexagons and open on the top, was 
implemented as a boundary condition. Open-ended nanotubes were first explained by the 
high local electric field in the region of nanotubes growth preventing their closure ends. 
More recently, catalysts such as Ni were thought to be attached to hold open the ends of 
nanotubes and promoted growth. It was assumed that nanotube growth occurs at the 
cathode surface (Fig. §.26) and that surface chemistry is controlled by the local 
terminated bond and not by the bulk nanotube bonds. Nanotube growth is based on the 
adsorption and desorption of three small carbon clusters (Ci to C 3 ) to simulate nanotube 
growth. The surface reactions considered here are li v d in Table §.20, where CP 
represents the nanotube radical site at the open edge of the tube and CNT is assumed to 
be the unique “bulk” species. In the reaction (SI) for example, carbon atoms C react with 
a radical site CR to generate a new radical site CR and one incorporated carbon into the 
bulk nanotube CNT. The pre-exponential factor of these reactions was fixed arbitrarily at 
2.50x10 ! 1 and was varied as a parameter from 3.0xl0 4 for carbon addition of C, C2 and 
C3 on soot [119] to 1.0x10 forcarbon addition on diamond [ , ]. Since no other 

carbon phases are incorporated in this model, the total number of sites T is calculated by 
geometrical considerations from Hamada’s indices n and m. We consider a triangular 
arrangement of carbon nanotubes within a bundle [ 13 °] (Fig. §.25) and divide the number 
of atoms per hexagon by the area of the hexagon using the 
3 N 1 

equation T = — 1 = — x with N at the number of atoms per unit cell 

3 "V 3 ( d CNT + d NT _ NT ) N Av 

(nanotube), dcNT the nanotube diameter, dNT-NT=0.340 nm the distance between two 
adjacent nanotubes [129] and N Av the Avogadro’s number. For example for a (10,10) 
nanotube, cIcnt=1-357 nm and N at =20, hence the site density is T = 6.66xlO" 10 moles/cm". 
For comparison, this is lower than the diamond site density [128] of 2.61x10" moles/cm . 
The surface site density is often assumed to be conserved; and the surface-species 

dZ s 

conservation equation is given by — — =— ^ = 0 where Zcr is surface species site 

dt T 

fraction and T is the surface site density [40]. The chemical production Scr of surface 
species by surface reactions is given by the sum over the rate-of-pro^ress variables for all 
surface reactions. The rate of production scr expressed in moles/cm 2 /s is converted to 
linear nanotube growth rate G in gm/s by using nanotube bulk mass density Pcnt = 2.20 
g/cm J and molecular weight Mcnt= 12.01 g/mol using the equation 
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§.3.4.5.d Boundary conditions 

The axial gas velocity at the anode was estimated from the measured erosion rate 

o 

O and the total gas density p by 11 anode ~ ^ where A is the anode surface area. For a 

stagnation-point flow problem, axial velocity at the cathode is zero. The anode 
temperature was always fixed at the average vaporisation temperature of graphite, 
Tanode=4000 K. The temperature of the cathode Tc is predicted as a part of the solution 
by adding a radiative energy balance as a boundary condition. Close to cathode surface, 
the diffusive heat flux in the gas-phase is balanced by the thermal radiative heat loss to 

d T 

the anode given by: A—— =<7£(T c i ~T^ N0DE )F CA where a is the Stefan- 

OX 

CATHODE 

Boltzmann constant, £ is the surface emissivity and Tanode is the anode temperature to 
which the cathode radiates. The radiative exchange form factor is Fca- At the anode, mass 
, / of (000 K. gi >cn a dih lion r a- 
x. This factor accounts for ihe mixing of anode material with inert atmosphere and is 

moles (C + Ni + Y) 




defined by the ratio r =- 


-xlOO . In the case of evaporation of the 


moles{C + Ni + Y + He) 
mixture of 94at.%C/5at.%Ni/lat.%Y, Krestinin and Moravsky [ 7,1 11] reported that 
carbon vapor is diluted as much as x=20-30 times by mixing with inert atmosphere. For 
the reduced model considered here, mass fractions at the anode were calculated using 
IVTANTHERMO [121] computer code. For the full model, mass fractions Yj were 
calculated from a zero dimension (0D) model involving the following set of ordinary 

BY. 

differential equations: p — =M , (O l : = 0 (i = l,...,w ) where the net chemical 

dt 

production rate tOj is calculated, as discussed above, from all the reactions listed in Table 
§• 1 . 

§.3.4.5.e Numerical results 

Numerical simulation of nanotube growth in the arc reactor has been carried out 
by solving the set of model equations and the boundary conditions using the SPIN 
computer code for the full chemical model. The chemical model of Krestinin and 
Moravsky, developed originally for fullerene synthesis, was tested in typical single 
walled nanotube growth conditions in the interelectrode region. For these calculations, 
helium used was as buffer gas, with a total pressure of P=660 mbar and interelectrode gap 
of 3 mm. The measured electric power dissipated in the arc for 100 A electric current was 
q=1.24xl0 0/ Watt/m 2 ; and the dilution factor at the anode was fixed at x=20. In these 
conditions, the measured erosion rate was 0=20.4x1 O' 03 g/s and the calculated mass 
density of the gas at the equilibrium temperature of 4000 K is p=9.24 xlO' 06 g/cnr . For a 
surface anode area of A=0.29 cm , the axial velocity u An ode=7818 cm/s was estimated at 
the anode from mass conservation. When small gaps are considered, the form factor Fca 
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should be accounted for. However, in these calculations it was set to unity. Calculated 
temperature profiles in the 3 mm interelectrode gap for pure helium are shown in Fig. 
§.27. A maximum temperature of ~ 6500 K is obtained at ~ 1 mm from the cathode and is 
in fairly good agreement with measured rotational temperature from optical emission 
spectroscopy performed by the Huczko group [49] (5300 K ± 500 K) in the center of the 
plasma with 1=54 A and P=880 mbar). Nanotubes are created between 1200 K and 1800 
K. [ I31 , lj2 ] Such temperatures are reached very close to the cathode, justifying the 
assumption of nanotube growth at the cathode surface. The calculated mole fractions of 
major carbon species are shown in Fig. §.28. Due to the high temperature that exists 
throughout most of the gap, the discharge is dominated by atomic carbon. The regions 
close to the cathode and anode are cooler than the center of the plasma. This enhances 
carbon atom recombination and explains the two peaks of C2mole fraction shown in Fig. 
§.28. The first peak is due to the recombination of C atoms close to the anode and the 
second peak occurs at 0.3 mm from the cathode, and is due to the competition between 
gas and surface chemistry. This second peak constitutes a feedstock for nanotube growth. 
Radiation intensity contours measured at 515 nm by the Huczko group [49] indicate 
qualitatively the presence of two zones rich in C2, in agreement with our calculations. It 
can be noted from Fig. §.16 that decreases of C, C2 and C3 carbon species close to the 
cathode are due to the surface chemistry, andthat full* rer.es start to be formed at 0.3 mnr 
close to the cathode. However, the (ID) model as developed here, is not able to calculate 
the distribution of fullerenes in all of the reactor. The region where fullerenes start to be 
formed is important from growth mechanism point of view and could be investigated by 
optical spectroscopy measurements. The calculated normalized species abundance at the 
cathode is plotted versus cluster size in Fig. §.29. It shows a roughly Gaussian 
distribution of even-numbered clusters with 28-58 carbon atoms and two peaks for C60 
and C70. Inserted in Fig. §.29 are the distributions of carbon clusters as measured under 
various experimental conditions for laser ablation of a graphite target by Kroto et al. [1] 
that show a similar behavior. The calculated number densities of major carbon species are 
listed in Table §.21. Number densities of the buffer gas and the catalysts calculated from 
equation of state are fairly constant nHe=l-4xl0 18 cm' 3 , nNi=2.0xl0 14 cm' 3 and 
n Y =3.2xl0 14 cm' 3 . 

Finally, the calculated nanotube growth rate using (10,10) is plotted in Fig. §.30 
for dilution factors 5, 10 and 20, as the pre-exponential factor of surface reactions varied 
from 3.0xl0 4 carbon addition of C, C2 and C3 on soot [1 16 s8 ] to l.OxlO 13 corresponding 
to carbon addition on diamond [127,128], The results show that the predicted growth rate 
ranges from a few (im/min to 1000 /xm/min. 

§.3.4.6 Two-dimensional models 

A two dimensional (2D) approach was developed by Bilodeau, et al. [115] for 
simulating carbon arc reactor for fullerene synthesis in helium or argon. The model 
solves velocities, temperature, and total concentration of carbon species in the two- 
dimensional axisymmetric space between the electrodes. The model accounts for carbon 
evaporation from the anode, deposition on the cathode, and condensation in the regions 
surrounding the arc. 
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§.3.4. 6. a Model formulation 

Bilodeau, et al. [1 15] assumed an axisymmetric laminar flow, steady state and 
local thermodynamic equilibrium (LTE). The anode erosion rate was assumed uniform 
and was obtained experimentally. The deposition on the surface of the cathode is 
governed by diffusion, with velocities of the order of 400 m/s. The drift of ionic species 
having a velocity 10 m/s due to the electric field was neglected. The input of energy in 
the arc is due to ohmic heating and to the enthalpy flux of the electrons. Enthalpy 
diffusion due to species transport is considered and radiation losses are considered using 
the net emission coefficient method with a plasma thickness of 0.5 mm. As shown in Fig. 
§.31, two calculation domains were chosen to predict temperature and velocities 
distribution. The first domain, noted (1) in Fig. §.31 represents the interelectrode gap of a 
thickness chosen equal to 1 or 4 mm. The second domain noted (2) in Fig. §.31 
represents the entire reactor. 


§.3.4. 6. b Model equations 

Continuity 
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Axial and radial . momentum c observation 
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Energy conservation: 

f (pfo)=f(— W—+|— £&“£(— -pD c ){h c -h g )§(O c -47ie n +S h (§.51) 

c <J 2 ec c 

p p p 

Carbon species conservation 


V . (pv(ti c ) =' V . (pD c Vw c ) + S m 


(§•52) 


In these equations, is the gas velocity, P is the local pressure, p. the viscosity, p the 
mass density, g the gravitational acceleration, h is the specific enthalpy, he and h g the 
enthalpy of pure carbon and pure buffer gas, respectively, k is the thermal conductivity, 
c p the heat capacity at constant pressure, e n is the net emission coefficient, e the charge of 
the electron, ks is the Boltzmann constant, a the electrical conductivity, o>c is the mass 
fraction of carbon species and Dc is the diffusion coefficient of the carbon species. The 
axial current intensity 

j x = aE is calculated from the linear electric field E, the conductance G, the total current 
intensity I, and the radius of the electrodes. 


£ = 
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B is the magnetic field intensity: 



(§•54) 


with po the magnetic permeability of vacuum. Sh and S m are respectively heat and mass 
source terms. The mass source term S m represents the carbon evaporation near the anode 
(positive term) or the condensation of the condensation of carbon species near the 
cathode (negative term). 
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§.3.4.6.c Numerical results 


The model equations presented above were solved using the SIMPLE method 
described by Patankar [ lj3 ]. Boundary conditions are the following: 

On the reactor wall, the temperature T=350 K, the axial and radial velocities 
V z =V r =0 and the carbon composition Xc=0. 


d(p 

On the reactor axis, the radial velocity Vf=o and radial gradients — — =0 ; with 

or 


<N(T, Xc, V z , J z ) 

At the anode and cathode tip T=3300 K on the side and T=2800 K on the axis. 


Typical temperature profile maps are shown on Fig. §.32 for helium and argon in the two 
simulations domains (1) and (2) discussed above. For a 1-mm gap the maximum 
temperature is -12,000 K. When the gap distance is increased to 4 mm, this maximum 
temperature increases to over 17,000 K near the cathode. Bilodeau, et al. attribute this 
difference to the lower carbon concentration and higher electrical resistivity. When argon 
is used instead of helium with a gap width of 4 mm, the temperature range (1500 K-5000 
K) is wider in argon than in helium. In the temperature range 2000-3000 K assumed by 
the authors to be favorable for the formation of fullerene precursors, the carbon species 
.coiicentratiomis-higb&r miielium. than in argon. . . 

§. 3.4.6. d Flow and heat transfer modeling 


More recently, Hinkov [1 16] compared the temperature profiles obtained in 
nanotube conditions for helium at 660 mbar and argon at 100 mbar. The model consists 
of Navier Stokes equations implemented in FLUENT computer code. For these 
calculations, only two species, namely carbon atoms and helium or argon without any 
chemical reactions were considered. A constant heat flux of 6.5x10 ergs.cm'-.s' for 
argon and 1.2xl0 10 for helium was considered from the measured power. The wall 
reactor temperature was 300 K. The anode and cathode temperatures are 4000 K and 
2000 K, respectively. From Fig. §.33, the calculated maximum temperature is 14700 K 
for argon and 9080 K for helium. The lower thermal conductivity of argon leads to a 
much larger hot zone close to the cathode. This result could explain the higher observed 
nanotube yield when helium is used rather than argon. 


§.3.5 Arc Modeling Concluding Remarks 

Mathematical modeling of the carbon arc has been the subject of numerous 
studies where the objective was to better understand fullerene and nanotube growth, as 
well as to predict the chemical composition of species and the temperature distribution in 
the discharge. Several approaches were developed and permitted isolating the effects of 
thermal, kinetic, mixing and electric processes on vaporization and condensation of 
carbon and metallic clusters. Compared with other high temperature techniques, the 
uniqueness of the arc is the presence of ions. When a metallic catalyst is mixed to the 
graphite in the anode, the major ions as measured and calculated are found to be the 
catalyst ion Y + or Ni + and not C + . Hence, in addition to its catalytic role, these metals 
play the role of current carriers in the discharge. In addition, the energy of ionization of 
the catalyst is weak; and the process of the anode vaporization is more efficient. 
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Otherwise, the experimental results indicate that the limiting process in nanotube growth 
is the rate of anode erosion. It will be interesting to explore the effect of other types of 
catalysts with lower energies of ionization. In addition, we can conclude that the low 
density of carbon ions in the discharge does not explain the fast growth of nanotubes. 
Hence, carbon nanotube growth occurs from neutral precursors and could be modeled by 
a set of surface reactions simulating open nanotube growth exposed to a mass flux from 
the anode. Since in this approach nanotube surface chemistry is controlled by the local 
terminated bond and not by the bulk nanotube bond, a mechanistic approach based on the 
formal resemblance between the bonding and the structure between open nanotube and 
other carbon surfaces was proposed to explain nanotube growth. Predicted growth rates 
are in the range of 100 to 1000 microns/min. Nevertheless, a debate still exists on the 
reason why the nanotube is maintained open during the growth. Some authors attributed 
this fact to the high local electric field in the growth region close to the cathode; but the 
predicted values of this local field varies within four orders of magnitude from 10 4 tol0 h 
V/cm. This suggest that further modeling efforts are needed to estimate with accuracy the 
electric field distribution in the connection with thermal, kinetic and turbulent mixing in 
multi-dimensional configuration. Further development of a modeling approach is in the 
direction of selecting more representative gas phase and surface reactions and using ab- 
initio calculations of mechanisms and kinetic data. There is a real need for accurate - *• ■ ■ 


kinetic' dafaTor the adsorption of small carbon clusters on open rianotubesbf different 
diameters and chiralities or on pentagon-heptagon defects of closed nanotubes. Since we 
have evidence that SWNTs may grow in the gas phase, further work in modeling the arc 
process will include gas phase production as proposed in Ref. [14]. There, gas phase 
reactions for the production of SWNTs are proposed to grow from carbon/nickel clusters 
using a formalism for the reactions like the one developed for the Boudouard reaction 
used in the high-pressure carbon monoxide process (HiPco) by Dateo et al. [28]. The 
calculations will be extended to two dimensions, and possibly 3D in order to take into 
account both collaret and core deposit growth in the cathode region. Even though, growth 
chemistry and surface site densities are very different for SWNT growth in the collaret as 
compared with MWNT formation in the hard deposit, the formalism discussed in this 
paper remains valid, and could take into account simultaneous growth of these species in 
2D or 3D configuration. 


§.4 Computational Fluid Dynamics Analysis of Transient Carbon 
Plumes in Laser-Ablation SWNT Production 

The use of computational fluid dynamics (CFD) simulation methods for the 
analysis of the flowfields resulting from the laser ablation of carbon targets in the 
production of single-walled carbon nanotube production is a relatively new practice, with 
most studies having been conducted only since the turn of the century. CFD techniques 
for the solution of the governing equations of fluid motion have been present since the 
advent of high-speed computers in the 1960’s. The National Aeronautics and Space 
Administration (NASA) and various aircraft corporations originally developed CFD 
methods primarily by for the solution of external flowfields surrounding aircraft and 
high-speed space vehicles leaving or returning to the Earth’s atmosphere. Additionally, 
CFD methods for the solution of internal flows such as those occurring through gas 
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turbine engines were developed at the same time as external solution techniques or 
“codes.” Recently, CFD methods have gained popularity in other fields of endeavor that 
are less “aeronautical” in nature such as the analysis of various forms of materials 
processing, or when the fluid dynamics equations are merged with equations of 
electrodynamic forces such as the Maxwell equations (in a group of techniques referred 
to as "magnetohydrodynamics" or MHD), to such diverse fields of study as the flow of 
charged particles in the magnetosphere of the Earth. Two excellent references for general 
conventional internal and external CFD techniques are the works of Hirsch [ lj4 , lj3 ] and 
Hoffman and Chiang [ l36 ]. The former publication is more theoretical in nature than the 
latter, however the latter provides a better introduction for the casual CFD code user. 

Recent work in the application of CFD methods to laser ablation flowfield studies 
have been accomplished by both the solution of the Euler equations [134] which ignore 
viscous effects in the solution technique and by solution of the Navier-Stokes equations 
[134] which do incorporate viscous effects. Both methodologies will be examined as they 
have been applied to laser-ablation SWNT production in the following paragraphs. All 
numerical techniques are approximations of the full partial differential equations of fluid 
motion, be they Euler or Navier-Stokes equations, and the solution of the governing 
PDE’s calculated by computer codes will be highly dependant upon the boundary 
conditions applied. AIKTansientrsokitton methods are highly dependant oirtwo faetors -*• - • 
uidai \ cuiidniv.iis and me lime step Sizetakeri af eacfTtini< -.cp. Transient metho 
rely on having accurate assessment of the current flowfield and then extrapolate the 
current solution into the future by known mathematical methods. Therefore, no 
discussion of the CFD techniques as applied to laser ablation flowfields would be 
complete without a discussion of the boundary condition determination for the flowfields 
- such discussion will follow the examination of the basic numerical solution 
methodology for the Euler and Navier-Stokes equations. 

§.4.1 Inviscid Solution of Carbon Plumes in Laser Ablation 

Lobao and Povitsky [ L ’ ] have successfully used a unique combined Eulerian and 
Lagrangian approach to the solution of the flowfields resulting from the rapid 
vaporization of carbon targets in the SWNT oven. Note should be made that the phrase 
“Eulerian” in this sense refers to analysis by a fixed control volume, not to the governing 
inviscid fluid dynamics conservation equations - although Lobao and Povitsky do indeed 
use the Euler equations in their methodology. Lagrangian methods differ from Eulerian in 
that individual particles are tracked as they flow through space; and overall 
thermophysical parameters are then calculated by statistical techniques applied to the 
many particles in the flowfield. 

Lobao and Povitsky make the assumption that the inviscid Euler equations suffice 
for an analysis of the laser ablation of carbon due to the rapid time frames in which the 
ablation process occurs. The laser irradiation times for the laser ablation process as 
practiced at NASA JSC are approximately 10 ns in duration, although the material may 
continue to ablate for several nanoseconds after the irradiation terminates. It should be 
noted however that viscous shear stresses, x, are a function of velocity gradients wherein: 


3v,. 

r„ = u — - 
,J ^dx, 


(§•55) 
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for a Newtonian fluid, where vj is the velocity component in the “i” direction, xj 
represents the “j” spatial direction, and p. is the viscosity. A direct physical dependency 
upon time in the viscous components of the flow physics does not exist then. Given the 
extreme velocity gradients that occur in the laser ablation process, with inlet velocities of 
the ablative carbon being on the order of kilometers per second flowing into essentially 
quiescent argon gas It is possible that even given the short time duration of the flow 
viscous effects may be of significant enough magnitude that they cannot be ignored. 


§.4.1.1 Solution Methodology 


Lobao and Povitsky used the two-dimensional compressible Euler equations with 
a generalized curvilinear coordinate system in £,, r\ coordinates. The equations can be 
represented in vector form by: 


dQ dF ^ 3 G _ 

where the vectors of conserved variables (Q, F, G, and S) are given by: 
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In this formulation, p represents the local gas density, u and v are the local velocity 
components, U and V are the contravariant velocities, E is the total energy, and P the 
pressure. C| and C 3 are the mass fractions of the two species. £ x , £, y , r| x , and % are the 
metrics of the coordinate transformation. J represents the Jacobian of the coordinate 
transformation, and 8 is the radial axisymmetric coordinate in the r| direction. The first 
equation represents the conservation of mass while the second and third equations are the 
conservation of momentum in two dimensions. The fourth equation is the conservation of 
energy, while the final two equations are the “conservation of species” equations for the 
two chemical species, C 3 and Ar, used in the Lobao and Povitsky analysis. The gas is 
assumed to be a calorically perfect gas in their analysis. Excellent references for further 
elaboration on the Euler equations and their generalized coordinate versions can be found 
in [135, ]. Lobao and Povitsky use a second-order upwind scheme for the conservation 

equations called “MUSCL” [134] and they use a central relaxing Total Variation 
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Diminishing (TVD) scheme best described in their paper [136]. The Lagrangian scheme 
used by Lobao and Povitsky to solve for particle motion in the flowfield is not well 
described in their paper other than the fact that it is based upon a second-order Runge- 
Kutta scheme, and yields particle streak lines in the resulting carbon plume. 

Lobao and Povitsky considered an axisymmetrical flow domain that was 25 cm in 
the longitudinal direction by 5 cm in the lateral dimension with a structured grid giving a 
grid size of 250 x 50 computational nodes. Initial conditions inside the flow chamber 
were a pressure of 1 atm and a temperature of 1500 K for the argon in simulation of a 
high-pressure ablation process and also for a pressure of 10' 6 atm for a low-pressure 
ablation case in their paper. The boundary conditions at the site of the laser irradiation are 
given as 100 atm of pressure at a temperature of 5000 K. The irradiation boundary 
condition was allowed to propagate for 20 ns with an initial time step of 10' 11 sec. The 
ablation boundary condition is then replaced by the flow chamber conditions and the 
entering plume of C3 is allowed to propagate downstream. Time steps in the Lobao and 
Povitsky method were controlled by adjustment of the Courant-Friedricks-Lewy, or CFL, 
number given the fixed even grid spacing. 

§.4.1.2 Results 

Lobao arid Tovitiky simulated a uumbf r of Variations upon the standard ' 
conditions used in the laser ablation at NASA JSC and qualitative comparisons were 
made between the parametric variations, although no comparisons were made to existing 
experimental data [ Ij9 ] for shock front expansion. As previously mentioned, two 
variations on the chamber pressure were simulated and showed consistent physical results 
in that the C3 plume expanded more rapidly in the lower chamber pressure than in the 1 
atm pressure case. 

Other variations on the standard operating conditions that Lobao and Povitsky 
examined included the influence of carbon injection velocity on the flowfield dynamics. 
Initially, the entrance velocity of the carbon plume was set to zero and the plume allowed 
to propagate by virtue of the pressure differential. Additional injection velocities were 
simulated in their paper although few results were shown. A final variation examined by 
Lobao and Povitsky was the effect of the interaction of multiple laser irradiations upon 
overall plume dynamics. The standard processing of SWNT calls for irradiation of 10 ns 
duration for the first pulse followed 50 ns later by a second pulse of 10 ns duration. This 
process then repeats on a 60 Hz cycle. Lobao and Povitsky allow for a similar double- 
plume injection however using a gap of 16 p.s between irradiation intervals. The main 
effect of the multiple plume interaction noted by Lobao and Povitsky was the increase in 
plume temperatures observed from the reflected shock waves re-compressing the plume 
from wall and plume interactions. 

§.4.2 Navier-Stokes Solutions of Carbon Plumes in Laser Ablation 

Greendyke et al [ 140 ) 14, j 142 I4 -'' ? 144 j ,4; ’] have conducted a number of studies of the 
carbon plume resulting from laser ablation of a carbon surface in SWNT production 
using a full Navier-Stokes solution of the flow equations. In preliminary work [139], it 
was determined that the viscous stress terms could not be ignored given that the 
background flow was approximately 3 cm/s yet the carbon plume entrance velocity was 
on the order of 5 km/s. such large differences in velocity would necessarily lead to large 
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viscous stress terms (as explained in the previous section) given the magnitude of the 
inevitable velocity gradients. 

In his first study [139], Greendyke used the VULCAN code [ 146 ], which was * 
originally designed for the study and analysis of internal supersonic combustion 
processes in hypersonic propulsion applications. The code was selected for its ability to 
handle mixed elliptical, parabolic, and hyperbolic flowfield equations. The code solves 
the generalized curvilinear Navier-Stokes equations described by: 

a +(E-E.\+(F-F,\+(a-a,\=s 

where Q t , E, F, G, and S are vectors describing the conserved variables included in the 
conservation of mass, momentum, energy, and species conservation equations as best 
described by White and Morrison [146], E v , F v , and G v represent the viscous flux vectors 
and are also described in [146]. In addition to the full Navier-Stokes equations, the 
VULCAN code also contains several turbulence model formulations, and Greendyke 
selected the k-co model of Pope (described, along with other turbulence models in [ ,47 ]) 
for his studies. The VULCAN code utilizes several solution methodologies for the 
Navier-Stokes equations, and Greendyke selected the time-accurate 4 th order Runge- 
Kutta methodology option of the VULCAN code that is based upon a method first .„ . , 
developed by Jameson et al [ ,4£ J~for the solution of the Euler equations. IirgeiteraL~ • — «« 

however, any CFD code capable of mixed flow regimes would work well in laser 
ablation plume development studies - a very good synopsis of modem CFD codes can be 
found in Laney [ 149 ] and at the accompanying web site [ 1;>0 ]. 

§.4.2.1 Boundary Condition Determination at Ablative Surfaces 

Before any CFD code can be accurately run for the solution of carbon plumes, it 
is necessary to correctly determine the boundary conditions at the ablative surface - 
without the boundary conditions correctly determined, no CFD code can predict reliable 
results. However, very little experimental information was available for the boundary 
conditions at the surface of the laser ablation carbon target. The only known 
thermophysical quantities had been the carbon mass ablation rate of 1.6x1 O' 6 gm/laser 
pulse in argon, and the approximate energy of 300 mJ over the 5 mm laser spot size. 
Greendyke et al’s original studies [140,141] assumed that the carbon mass was injected 
into the flowfield through the 10 ns duration of the laser pulse. A simultaneous solution 
of both the ideal gas equation 


PV = NRT 
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and the Clausius-Claperon equation: 
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resulted in a density of 10.59 kg/m’, a gas temperature of 521 1 K, a vapor pressure of 377 
atm., and a carbon injection velocity of 1900 m/s at the target surface. In those studies 
however, the resulting propagating carbon plumes exceeded the shock front locations 
data gathered by Puretzky, et al. [139] past the 200 ps post-ablation time. 

In later studies [142,143], the ablation time was assumed to be 15 ns - the 
additional 5 ns worth of ejected ablation material coming after the termination of the laser 
pulse. The resulting solution of ideal gas and Clausius-Claperon equations showed that 
while the temperature of the injected carbon was relatively insensitive to the additional 
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ablation time (T=4950 K), the vapor pressure of the carbon plume dropped considerably 
to approximately 100 atm. The corresponding density was an order of magnitude less, 
1.039 kg/m ', and the plume injection velocity rose to 5228 m/s due to the lower 
molecular weight. Chemical equilibrium was assumed at the ablation surface and the 
CHEMKIN [40] code was used to determine inlet carbon species mass fractions for C 
through C 5 (no data being available in CHEMKIN for C 6 ). C 3 and C 5 were found to be 
the dominant species at the ablation surface inlet with the mass fraction of C 3 being 0.463 
and C 5 being 0.426. The next contributors to the inlet flow were C 2 and C 4 with species 
mass fractions of 0.044 and 0.059 respectively. The atomic carbon mass fraction was 
nearly negligible at 0.008. 

§.4.2.2 Flow Domain Gridding 

After the correct determination of flowfield boundary conditions, the domain of 
interest needs to be resolved into a computational grid for the calculation of flowfield 
characteristics. Greendyke [140] originally used a 2D axisymmetric solution 
methodology for the flowfield volume with a 2D grid. At the time of his first 
investigation into laser ablation carbon flowfields, little was known of the resulting 
flowfield and the degree to which the carbon plume would propagate into the flowfield 
abwftsfreafh of ffie'fc'arb> a fargefcf or «he region n the 'outer quartz tube of the laser 
ablation oven. Fig. §.34 shows tiie original axisymnietric grid that was used by 
Greendyke in the original investigations. Fig. §.34 however displays a symmetrical grid 
about the y=0 axis. Only the portions for y > 0 are used in the actual axisymmetric 
flowfield calculations. One problem encountered was that the grid only extends to 10 cm 
in front of the carbon target and boundary conditions at the upstream portion of the 
background argon flow reflected the pressure wave resulting from ablation boundary 
conditions back into the propagating carbon flow. In later studies, Greendyke et al [142- 
145] used a “primary grid” (Fig. §.35) with 0.5 mm grid node locations for the first 10 cm 
in front of the ablation surface for the inner quartz tube alone as in the first work [140- 
141], with a “secondary grid” of 0.5 mm vertical spacing and 5 cm horizontal spacing out 
to a point 55 cm upstream of the carbon target. The elongated flowfield region under 
consideration alleviated the problem of reflecting pressure waves returning to influence 
the carbon plume’s propagation into the argon background flow. 


Multiple gridding software options are available to the researcher, with varying 
degrees of complexity to use. The most convenient of the available gridding software 
codes is the Mesh Generator software [ l31 ] available as an “add-on” package to Amtec 
Engineering Inc.’s Tecplot & software [ l32 ]. The Mesh Generator ™ is capable of outputting 
grid information in both the ASCII format and the PLOT3D form (PLOT3D is a graphics 
package for CFD use that is available from the Open Channel Foundation - see their web 
site at [ l33 ]) used by many commercially available CFD codes such as the VULCAN code 
[146], 

§.4.2.3 Flowfield Solution Procedures used by Greendyke et al 

In all of the studies by Greendyke, et al. [140-145] the same procedure 
was used throughout. After determination of the boundary conditions, the first priority of 
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the simulations was the solution of the “base” argon flow. The global elliptical viscous 
method of the VULCAN code was used with diagonal approximate-factorization in 
axisymmetric fashion for this part of the simulation. Full k-co turbulence modeling and a 
mixture of thermally perfect gases options were chosen. The solution was allowed to 
converge until a sixth-order reduction in the L2 norm was observed. 

Actual simulation of the carbon plume resulting from the laser-ablation of the 
target could proceed with the base flow established. To simulate the carbon plume a jet of 
carbon was allowed to enter at a density, with the boundary conditions described above, 
from the first five grid point locations in the grid block immediately in front of the carbon 
target corresponding to 2.5 mm in actual length. Flow enters axisymmetrically having a 
radius of the laser beam’s spot size. The time-accurate Runge-Kutta solution procedure 
was applied to the restart files from the base flow simulation for 100 (or 150 in later 
studies) time steps of 0.1 ns each — corresponding to the actual duration of a laser pulse. 
The restart files from this “carbon-injection” simulation were themselves used in 
subsequent simulations with the carbon mass flow injection deleted in the input data file 
for VULCAN. An isothermal wall boundary condition was used in place of the carbon 
injection condition with a temperature set to slightly higher (1773 K) than the initial 
ambient conditions in the nanotube production chamber. This boundary condition was 


held constant until 60 ns ot simulated trmt had passed and then' the ablation boundary - 
condition was re-iniiiated for another i5 ns of simulated time io conespond to the second 
laser ablation of the processing on a 60 Hz cycle as practiced at NASA Johnson Space 
Center. The flow with the developing carbon plume was simply allowed to develop 
naturally to whatever point in time required by manipulation of the time step size and 
maximum number of iterations allowed in the VULCAN input data file. 


§.4.2.4 Results of Navier-Stokes Simulations 

For the original work with full Navier-Stokes solution methods, Greendyke et al 
used only C3 as the injected species at the carbon ablation boundary points to determine 
flowfield characteristics and plume propagation into the background argon gas. No 
chemical reactions were allowed in these studies. Fig. §.36 shows the result of the C3 
injection simulations at 200 ps post-injection time, and indicates the classic “smoke ring” 
profile of the carbon plume observed by multiple experimental investigations. Further 
confirmation of the accuracy of the Navier-Stokes modeling was presented by 
Greendyke, et al. [143] using only C3 as the injected carbon species without chemical 
reactions.. Good comparisons were obtained with experimental observations of carbon 
plume propagation conducted by Puretzky, et al. [139] for the earlier post-injection times. 
As the flowfield solutions were allowed to progress however, computational results for 
the leading edge of the carbon plume began to overtake experimental results at the 1- 
millisecond point of post-injection time. Temperature contours compared favorably to 
experimental results as well in this study as can be seen in Fig. §.37. 

In later studies [142-145], Greendyke et al added chemical reaction mechanisms 
to the flowfield solution methodologies. The first reaction model [141] was a twelve 
species, 14-reaction model for the formation of carbon molecules up to Cc including the 
ion species of C + , C/, and Cs + as well as electron concentrations. The chemical reaction 
model was a reduced form of the reaction mechanism for the formation of C60 and C70 
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fullerenes developed by Krestinin and Moravsky [8], The primary purpose of the 
inclusion of higher species and chemical reactions was two-fold. The first purpose was to 
see what effect condensation of the carbon species would have on flowfield dynamics 
and temperature profiles. The second purpose was to determine if formation of the C6 
molecule, since it is not injected in the boundary conditions, would provide some 
indication to the thermophysical conditions for higher carbon molecule formation. The 
formation of C6 would then serve as an “indicator species” for more complex carbon 
molecule formation. One discrepancy in this initial study with chemical reaction 
mechanisms was immediately observed - temperature profiles with chemical reaction 
mechanisms quickly exceeded the temperature profiles from earlier work - this result was 
first thought to be the result of energy being released into the flowfield from the 
condensation of lighter carbon species. However, in later work, Greendyke, et al. [144] 
discovered that the high temperatures were the result of excessive time steps being used 
in combination with chemical reactions. If the time steps for the initial laser ablation 
boundary conditions were decreased to 0.01 ns, temperature profiles in argon background 
gases returned to similar, although not identical, profiles for the original work with C3. A 
curious result of this work was that even with the higher temperature profiles, leading 
edge propagation of the carbon plume was comparable to the earlier work [141], but 
de dated'TbTri Pnretzky x t af s exjritr^entbl’v 0T’:'[f3 c 'j at post-ablation -oirset times nr- or 
beyond, 1 millisecond. A possible cause foi the higiiei plume propagation speed ol the 
calculated results than experimental could be the curve-fit data incorporated into the input 
files used in the VULCAN code for diffusion of the higher species of carbon. No data 
was available for the higher carbon species and data for C and C2 was simply 
extrapolated for higher species. As is generally known, lighter species will diffuse faster 
than heavier species. 

It is impractical to use the full C60 and C70 model for fullerene formation of 
Krestinin and Moravsky. All flowfield solution methods such as the VULCAN code have 
to solve the species continuity equations at each grid point in the considered flowfield 
domain. Computational intensity therefore increases with the number of grid points and 
species by the relationship: 


CPU °c (v*s ) 2 

where CPU indicates CPU processor time, N is the number of grid point locations in the 
domain, and S is the number of chemical species. As a result, for the flowfield grid 
described in [138], computational times of up to 30 minutes were required for a single 
iteration at a time step of 10 13 s thereby rendering such CFD studies impractical. 

In the same studies [143-145], work began by Greendyke et al on duplicating a 
series of parametric experimental studies conducted at NASA Johnson Space Center [ 154 ]. 
The latest work [145], represented the culmination of the computational simulations to 
duplicate the set of experiments. For this study, the 12 species model originally used was 
replaced by an 11 species model with 14 reactions - which is given in Table §.22. The 
goals of this study were to gain insight into the thermophysical characteristics of the 
flows under differing pressures, background gas temperatures, different background gas 
chemical species, and differing diameters of the quartz inner tube of the laser ablation 
oven. A secondary goal was again to attempt to define “indicators” in the flowfield 
solutions of the carbon plume that would provide insight in the refinement of the laser 


39 


ablation SWNT processing and allow improvement of the method for increased 
production. Unfortunately, in the set of experiments that this study was based upon it was 
not possible to determine the actual quantities of SWNT produced in the experiments. 
Only the total carbon ablated was available, thereby negating the ability to do a direct 
comparison between computational studies and SWNT production. 

However, certain insights were still available from the computational study. One 
of the main results of Greendyke’s study was a direct correlation between plume 
expansion into the background gas and the amount of recovered carbon from the 
experiments. A plume expansion factor, (j), was defined by Greendyke as: 
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oi, the. local suniraatkn^of thedrufividual-carbun -spades’ mass fractions. The plume, 
expansion factor provided a relative measure from one case to the next Tthc propagation 
of carbon throughout the flowfield assuming the carbon plume originates at the origin, as 
it does in all cases in [145], The plume expansion factor provided a direct relative 
measure of the expansion of carbon throughout the flowfield. 

Varying the thermophysical characteristics of the background gas resulted in 
widely different plume development at similar points in time. Figs. §.38-40 indicate the 
difference in both leading edge propagation of the carbon plume for argon, helium, and 
N 2 background gases at 100 ps post-ablation time. The resulting plume expansion factors 
for the variation of background gases is seen in Fig. §.41, which shows that the higher the 
plume expansion factor, the lower the mass of recovered material from the laser ablation 
when compared to the experimental results of Arepalli et al [154], Similar results were 
obtained for the other variations of background gas temperature, pressure, and the inner 
tube diameter of the laser ablation oven with one noticeable discrepancy - the case where 
the inner quartz tube was reduced from its nominal 1 inch diameter used in standard 
processing to a diameter of Vi inch. Experimental results showed that the narrower 
diameter case should result in a reduced quantity of mass being recovered from laser 
ablation, while the inverse relationship seen with the plume expansion factor plotted in 
Fig. §.42 indicates that the mass recovered should have been higher than the nominal one 
inch diameter case. With the exception of this one case, the plume expansion factor did 
provide the best “indicator” yet of the amount of ablative carbon material that could be 
recovered from laser ablation studies. However, the plume expansion factor did not 
provide a direct quantitative relationship for the carbon recovered, but it did provide a 
qualitative relative measure from one case to the next. A possible cause of this 
discrepancy may be the result of the high concentration of carbon in the small diameter 
case resulted in more absorption of laser energy by the ablation products. This reduced 
the energy flux to the target, resulting in less ablation. 

Other results from the latest study [145] did provide thermal and species 
concentration profiles for the parametric variations that can be used in separate de- 
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coupled solutions with more complex chemical models (such as the Krestinin and 
Moravsky model in [8]) using the CHEMKIN software [40], At the present time, 
Greendyke and Scott are continuing these decoupled studies with a refined grid mesh in 
the immediate ablation area of the carbon target for times ranging from ablation onset up 
to the first few microseconds after ablation onset. The flowfield solutions at different 
points in post-ablation onset times that result from this study can be combined into a 
single solution file via use of the Tecplot® [150] software. The resulting solution data file 
can then be input into the CFD Analyzer™ [ 15:> ] add-on feature to the Tecplot package. 
Once the user has specified the location of basic flowfield properties in the file, such as 
pressure and temperature, the CFD Analyzer™ package can automatically integrate 
streaklines or particle trajectories over the time domain of the composite flowfield 
solution file. All flowfield properties and chemical compositions along the streaklines are 
then returned by the software for analysis with the decoupled methodology previously 
described. Results from decoupled calculations along streamlines will be discussed in the 
next section. 


§.4.3 Chemical Kinetics Along Streak lines In Pulsed Laser Ablation 


For CFD problems having more than one dimension, adding a large chemical 
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weeks of computer time to solve. In many cases it is useful to decouple the problem and 
solve the fluid dynamics and energy equations without chemical reactions. The history of 
temperature, pressure, and mixing (dilution) with an external gas can be determined along 
streak lines. For steady state problems such as the arc or HiPco process, as was described 
in section §.5.3.1, one can easily integrate the position divided by velocity to obtain the 
time at corresponding points in the flow. The CFD solution provides all other properties, 
except the detailed chemical species. One can then use a chemical kinetics code, such as 
the AURORA code of CHEMKIN to solve the chemical rate equations along the 
streamline, as a function of time. For time dependent solutions of the flow equations, as is 
needed for dynamic situations such as the laser ablation process, a similar technique can 
be used to follow a sample “particle” of the flow as time evolves. However, the process 
requires having the solution at various times during the development of the flow. With 
these solutions described in §.4.2 one can do a similar integration of the distance/velocity. 
For example, in some preliminary calculations we have the evolution of the solution 
along a number of streak lines that start next to the laser ablation target. These streak 
lines are shown in Fig. §.43. Integrating the distance divided by velocity we obtain the 
time interval between point and therefore the time at each location. Then we have the 
flow field properties as a function of time. The temperature (Fig. §.44) and pressure 
profiles were then input into the AURORA code to solve the chemical rate equations for 
the evolution of the species along the streak line. One such preliminary solution is shown 
in Fig. §.45 for the complete model along one of the streak lines. The rate coefficients for 
this calculation were not accurate; therefore, it only represents the methodology, not 
accurate results. 


§.4.4 Conclusions from CFD Modeling 

The present limitations of reduced carbon chemical reaction mechanisms and the 
resulting computational overhead that would result from higher-order chemistry models 
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severely hinder the analysis of laser ablation SWNT production by CFD methods when 
CFD is used by itself. However, this does not mean that useful information for the study 
of SWNT production is unobtainable - the use of the plume expansion factor formulated 
by Greendyke, et al. [145] and outlined in the previous section does yield insight into the 
possible enhancement of SWNT production by aerodynamic means alone. The 
information provided - chemical profiles of “building block” species, temperature 
profiles, relative localized densities of the chemical species, etc. can provide input data 
for de-coupled studies with chemical reaction codes that do not have the dependency of 
CPU time on grid resolution and number of grid points considered. Indeed, it is possible 
that in the future, unconventional CFD methods such as Direct Simulation Monte Carlo 
methods [ 157 ], or molecular dynamics methods [ l;>8 ], could more completely simulate at 
least the initial onset of ablation in the flowfield since they do not have the same 
relationship between the number of chemical species and computational overhead - 
although such methods have difficulties of their own that make their use questionable at 
the present time. 


§.5 Computational Simulation of the HiPco SWNT Production 



A new method for the production of SWNT that is receiving much attention at the 
current time is the “High Pressure Carbon Monoxide” or HiPco process, first developed 
by Nikolaev, et al. [ 159 ]. Two of the main problems in the production of SWNT are the 
ability to produce them in economically viable quantities and to produce them with 
predetermined structures that have unique properties such as their chirality. While the 
HiPco process does not necessarily address the latter problem, it is a possible answer to 
the former. The HiPco process is a continuous gas production reactor that utilizes a 
relatively cold input stream of CO and iron pentacarbonyl, Fe(CO)s, that is injected into 
the reactor with impinging streams of heated CO in a “showerhead” configuration (Fig. 
§.46). As the iron pentacarbonyl is heated by the impinging streams, it decomposes 
releasing the Fe atoms, which then condense into metallic iron clusters. The clusters 
serve as catalysts for SWNT development. The SWNT begin to develop from the iron 
clusters when CO molecules interact with the iron catalyst and form an iron cluster with 
the developing SWNT forming on it and CO 2 molecules. 

Computational simulations of the HiPco using reacting-gas computational fluid 
dynamics (CFD) have not been common. One of the primary investigations into the 
process has been the work of Dateo, et al. [27] and Gok 5 en, et al. [19], who conducted 
extensive work first on a reduced chemical reaction model for the HiPco process and then 
combined this model with a simplified CFD code to produce a parametric analysis of the 
HiPco process. Work by other researchers such as Povitsky and Salas [ 160 ] have focused 
on the actual aerodynamics of the HiPco process to insure the best mixing of the iron 
pentacarbonyl and carbon monoxide streams for SWNT production. Scott, et al. [31] also 
conducted research on the HiPco process with limited CFD methods that utilized a 
chemical reaction model that is decoupled from the flowfield solution method. 
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§.5.1 Reacting Gas Modeling of the HiPco Process 

Dateo, Gok^en and Meyyappan published a two-part paper [27,19] that modeled 
their chemical reaction mechanism [27] in the first part, and the application of that 
chemical reaction mechanism [19J couples with a CFD solution methodology for the 
analysis of the HiPco reactor chamber. The first paper established their chemical reaction 
mechanism that can be broken down into 3 fundamental parts: the decomposition of the 
iron pentacarbonyl steps, iron catalyst growth and evaporation, and finally, carbon 
nanotube formation. 

Dateo, et al. noted that the iron pentacarbonyl will begin to decompose at 
temperatures above 500 K, which is the situation as the inlet Fe(CO)s stream initially 
encounters the hotter stream of CO. The decomposition mechanism is modeled in five 
steps: 

Fe(CO) 5 Fe(CO) 4 + CO 
Fe(CO) 4 oFe(CO) 3 + CO 
Fe(CO) 3 oFe(CO )2 + CO 
Fe(CO) 2 <-»FeCO + CO 
FeCO + M <->Fe + CO + M 

Chemical reaction rates for the dissociation were estimated by Dateo et al by use of 
experimental. data from Lewis et al ['*'] andJingdldug and Lineberger [ u ':]. , fhe_nivexse_ 
reaction rates were taken from Seder et al [ 163 ]. 

Iron cluster formation in the second part of their model was based upon the rates 
of Krestinin, et al. [29] wherein their model is based upon the reaction mechanisms: 

Fe m + Fe n ^ Fe m+n (n = l-4) 

Fe m ^ Fe m _i + Fe 
FeCO + Fe m — ^ Fe m+ j + CO 
FeCO + FeCO Fe 2 + 2CO 

It is noted that iron evaporation will occur at higher temperatures, while the condensation 
of carbon into SWNT requires higher temperatures for proper formation. 

Condensation of the carbon atoms is accomplished by a modification of the 
Boudouard reaction [ 164 ]: 

2CO(g) C(s) + C0 2 (g) 

into the three-step process suggested by Boudart [ l6;> ]: 

Fe m + CO — ¥ Fe m • CO 
Fe m • uO — ^ Fe m + CO 
Fe m • CO + CO — > Fe ra • Cnt + C0 2 

The final composition of the reaction model in [27] consists of 971 species in 1948 
chemical reactions. Since computational effort in chemically reacting flowfields is 
proportional to the square of the number of species multiplied by the number of grid 
points, the initial model in [27] was reduced for CFD analysis. The reduction was 
accomplished by grouping the iron cluster species that contain three or more atoms of 
iron into a single species. This re-grouping of the iron related species results in a much 
more computationally feasible model that has 14 species and 22 reactions - well within 
the capabilities of most modem CFD codes. 

Both the full and the reduced models were used in a simple “0-D” parametric 
study that examined variations of temperature, pressure and Fe(CO)s partial pressure at 
the gas injection interface. Dateo et al’s main findings in part I of their paper [27] 


43 


determined that increasing the temperature in the HiPco chamber speeds up the 
decomposition of the iron pentacarbonyl, but at the same time increases iron cluster 
evaporation. Increasing the pressure was found to speed up the Boudouard reaction and 
slow down iron pentacarbonyl decomposition. Increasing Fe(CO)s partial pressures 
speeds up both decomposition of the iron pentacarbonyl and the formation of iron 
clusters. When the reduced and full chemistry models were compared in the 0-d model, it 
was found that the two models gave acceptable comparison to each other over the range 
of parameters studied. 


§.5.2 CFD Modeling of the HiPco Process 

In part II of the work [19], Gokgen et al incorporated their reduced three- 
part chemical model for SWNT formation into a Navier-Stokes flowfield solver for 
axisymmetric nonequilibrium flows. The fundamental equations they used were 
expressed in vector form as: 

dU dF 1 d(yG ) II7 H 

at dx y dy y 
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The first 15 equations in the set represent the mass conservation equations for overall 
mass and the 14 species. The latter equations are the conservation of momentum in two 
directions and the conservation of energy equation. The H-vector reflects changes in the 
fundamental equations that are required by the axisymmetric geometry, and the W-vector 
term in the overall equation represents the source vector for chemistry effects. The 
equations were discretized into a finite-volume approach using a Steger-Warming flux- 
vector splitting [ l06 j. The coupled gas dynamics and flowfield dynamics solution method 
utilized an implicit Gauss-Seidel relaxation method [166]. The domain of the HiPco . 
reactor wass pli t in t e r twerbfockspone modeling the- inlet conical “showerbewT 'region,'- 
and the second for the cylindrical downstream portion of the flow. Laminar flow was 
assumed throughout the study. 

Results from the CFD simulation showed that there was an optimum pressure and 
temperature for the growth and development of SWNT that agreed with their previous 
part I of their paper as stated in the previous section. In addition, the variations of 
temperature and pressure showed good agreement with experimental results. An 
interesting conclusion of their study, (not elaborated upon in the paper) is the belief that 
improved mixing in the flowfield domain of the HiPco reactor would not lead to greater 
iron cluster (and therefore SWNT) production. Gokgen, et al., however, allowed that 
work was needed to improve the chemical modeling before they could conclusively 
respond about mixing improvement. 


§.5.3 Other HiPco Modeling Efforts 

Povitsky and Salas [159] also attempted a CFD analysis of the HiPco flowfield 
regime with the goal of improving the mixing of the iron precursors. Their technique was 
to use a combination of an Eulerian-based flowfield solution code named FLUENT [ l67 ] 
for the Navier-Stokes solution of the flowfield without chemical reactions, and a 
Lagrangian solution method to track the particles of the flow and the resulting parameters 
of the flowfield through which they pass. The FLUENT code utilized a k-E turbulence 
model in the solution procedure as implemented by Povitsky et al. In addition, FLUENT 
also provided the Lagrangian methodology for the particle trace simulations. 


The primary purpose of the Povitsky work was to examine multiple geometric 
configurations of the HiPco reactor including the angle of the impinging CO jets and the 
number of jets utilized to determine the optimum configuration of the “showerhead” 
reactor chamber. An interesting result of their work showed that rapid mixing of the flow 
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geometrically does not imply rapid temporal mixing of the flow. The best geometry for 
good heating of the iron precursor flow through the center of the reactor chamber was 
found to be a single central jet containing the precursor gas with three inlet jets for the 
heated CO with an angle of 45 (30?) degrees between the centerline and the impinging 
hot CO jets. 

§.5.3.1 Production of SNWTs/C0 2 

The solution to the chemical kinetics of nanotube production for the HiPco 
process was developed using two methods. In the first method, Dateo, et al [27] and 
Gok 9 en et al. [19] implemented their kinetics model first by applying it parametrically to 
constant temperature and pressure conditions. They then reduced the model and followed 
a similar procedure to show that their reduced model produced about the same amount of 
CC> 2 as their full model. The then included these models in an axisymmetric full Navier- 
Stokes CFD solver and concluded that their reduced model adequately represented the 
results of the full model. (See details in §.2.3.) The second approach was a decoupled 
method in which Povitsky [ l68 ] solved Euler equations with a turbulence model to find 
streamlines. Scott, et al. [31] then used the solutions along streamlines the chemical 
kinetics code of CHEMKIN to find the evolution of CO 2 and other species in the HiPco 
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investigated variations of those models. The time history of the temperature and dilution 
due to mixing were found along streamlines calculated by Povitsky [168] and Povitsky 
and Salas. 169 Sample stream paths are shown in Fig. §.47. The variations were developed 
to investigate the effect of various assumptions about the nucleation rate and the number 
of clusters in the model on carbon creation and cluster size distributions. Fig. 48 shows 
the evolution of carbon nanotubes attached to clusters of iron as a function of iron cluster 
size. One can see that at long times, the largest cluster starts to increase relative to the 
other clusters. This indicates that larger clusters should be included in the model. 
However, the amount of carbon nanotubes is dominated by the small cluster population. 
Therefore, from the standpoint of predicting carbon production, the model is quite 
adequate and does not need augmenting with larger metal clusters. However, the models 
tend to over predict the amount of carbon produced as compared with experiments, as 
seen in Fig. §.49. 

§.5.3.2 Assessment of Nickel as Catalyst for HiPco 

Due to the fact that iron tends to evaporate at fairly low temperatures, compared 
with nickel, it was anticipated that nickel might be a better catalyst for the HiPco reactor. 
However, when nickel tetracarbonyl was injected as the catalyst into the reactor, no 
carbon nanotubes were produced. The reason was not clear, since nickel is very effective 
in the arc and laser ablation processes. Modeling the nucleation and evaporation of metals 
allowed Scott and Smalley [42] to conclude that the higher binding energy and stability 
of nickel, as compared with iron was the principle factor that contributed to no growth of 
nanotubes when nickel-only was included as catalyst. A comparison of calculations of the 
production of SWNTs took into account iron-species rate coefficients and rate 
coefficients adjusted for the binding energies of nickel. It was found that almost no 
nanotubes were produced when the values corresponding to the nickel bond energies 
were used in the rate coefficients for evaporation and for the NiCO bond. It was found 
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that the metal-CO bond energy has a major influence on the production of CO 2 . The 
higher NiCO bond energy almost eliminates production of CO 2 , whereas, the smaller 
bond energy of FeCO leads to significant production. Also, the rate of dimer dissociation 
and small cluster evaporation affects the production of CO 2 by limiting the rate of cluster 
growth, and thus the number of clusters available to catalyze SWNT growth. This 
suggests that if CO reacts with metal clusters and removes atoms from them by forming 
MeCO, which has the effect of enhancing the evaporation rate and reducing SWNT 
production. It appears that the stronger bond between nickel and CO tends to extract 
nickel from clusters. It is like an enhanced evaporation that prevents nickel cluster growth 
in high pressure CO. Therefore, very few clusters are formed that can catalyze SWNT 
formation and growth. 

§.5.3.3 Production of Catalyst Particles Prior to Injection into HiPco Reactor 

An investigation of a potential procedure for injecting premade catalyst particles 
was studied using models of the HiPco process. It was envisioned that a laser could 
dissociate Fe(CO )5 in the injector of the HiPco apparatus. Laser dissociation instead of 
thermal dissociation would free iron atoms for nucleation and growth prior to mixing 
with .hot CO. One might be able to control the size of the iron clusters and thus affect the 
.prft<teaicn>3f^^ iavorahlc . .fabricating tbe deader il^uld-.biu., . 

reasonable just to be able to compute the results and assure that the experiment is worth 
building. A calculation was made assuming that all iron carbonyl is completely 
dissociated prior to injection into the main reactor. This is a truncated model (clusters up 
to n=40). The inlet flow is then allowed to relax to form iron clusters and iron clusters 
with CO attached. From Fig. 50, one can see that it takes about 0.1 to 1 ms for a fairly 
uniform distribution of iron clusters to form at this condition. The model only allows CO 
to bind to Fe n clusters for n>10, therefore, the small clusters have no CO attached, 
whereas, the larger ones all have attached CO. The possibility of the clusters catalyzing 
SWNT formation depends on the residence time in the cool (400 K) zone. The design of 
the inlet can be tailored to meet a wide range of residence times by adjusting the inlet 
tube diameter and length. It is apparent that a CW laser might result in over heating, 
therefore, a pulsed laser might be better. It could be timed such that a slug of gas that is 
dissociated then flows into the main reactor in the time between pulses. Another 
adjustable parameter could be the partial pressure of Fe(CO) 5 , controlled by controlling 
the temperature of its container. If small clusters are needed, then short residence times 
are indicated; but if larger clusters are needed, then longer residence times would be 
preferred. 

§.6 Conclusions 

Modeling of fullerene and carbon nanotube formation involves understanding 
how the species react chemically and evolve in time. Time evolution of conditions 
modeled using various degrees of fidelity in solving the fluid dynamics equations of 
motion through one degree or another using computational fluid dynamics. The kinetics 
involves knowing how they are formed from simpler species, that is, a concept of how 
they grow. This is not an easy question to answer, particularly in the case of single- 
walled carbon nanotubes. In this chapter, we have assumed the fullerene chemical 
kinetics model of Krestinin, et al. that was, to some extent, validated by comparisons with 
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measured production. We have applied it to the arc process and the laser ablation process; 
and we have used it as the basis for carbon cluster growth in models of SWNT 
production. Growth of SWNTs requires the introduction of metal catalysts; and we 
included the nucleation and growth of nickel clusters to represent catalysts in arc and 
laser ablation processes. Unfortunately, the kinetics of nickel cluster growth is not well 
documented in the literature; therefore, we have taken extensive liberties to approximate 
the kinetics of nickel. We have posed the mechanisms of cluster formation and growth 
using the aerosol theory developed by Girshick, which essentially treats cluster growth as 
a kinetically controlled process. Evaporation of these clusters was treated in their theory 
as determined from knowledge of the bulk vapor pressure and surface tension. Several 
variations of the kinetics model of the HiPco process for SWNT production were 
reviewed. The essential results at the current stage of development indicate that these 
models can predict measured trends in production, but the absolute quantity of carbon 
(nanotubes) is over predicted. 

We have reviewed various degrees of approximation in the kinetics models, and 
have shown attempts to reduce the size of HiPco and the carbon vapor models. 
Reasonable comparisons between the full and reduced models were shown. 

It should be emphasized that the kinetics models are approximate and preliminary. 

Tfey-a^iiL.vaJida non. andicfiiyihcnl baagfl ni?_pi)(inparisons with j leasuiements^S- 

pointed out, the nickel evaporation model can be refined by using measured cross 
sections and accounting for internal energy. It is possible that the intermediate clusters 
should be combinations of carbon and a few nickel atoms, rather than forming carbon, 
then nickel, followed by combining the two types of clusters. The present models show 
how we were guided by the observation that metal clusters are seen in the product. 
However, there are molecular dynamics models [' /0 ] that carbon and nickel co-condense. 
That model indicates further coalescence of nickel and carbon into larger cluster that may 
be followed by segregation of carbon from the metal catalysts. There are also some 
electron diffraction measurements that show crystalline metal carbides in the clusters 
seen in the product. 

We have shown computational fluid dynamics models of the arc, laser ablation, 
and HiPco processes. The highest dimension of any of the calculations, thus far, is two- 
dimensional. It is apparent from experiment, that the flow fields are three-dimensional. 
Thus, future work should be extended to 3D for problems when he chemistry is 
decoupled from the flow field. However, it appears that even ID calculations can yield 
important insight and understanding of the flow and its chemical kinetic 
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Table Captions 


Table §.l. Gas phase chemistry for Fullerene Model* 

Table §.2 Thermodynamic Coefficients for Carbon Clusters 
Table §.3 Reduced Fullerene model reaction scheme 

Table §.4 Thermodynamic Property NASA Coefficients Reduced Fullerene Model 
plus some other species used in plasma simulations 

Table 5 Iron Carbonyl and Iron Cluster Model of Krestinin, Smirnov, and Zaslonko 
[29] 

Table §.6 Rate coefficients for iron pentacarbonyl decomposition and 
recombination 

Table §.7 Fits to Vlasov Fig. 2 Recombination Rates 

Table §.8 bull iron pentacarbonyl and iron cluster HiPco model of Dateo, et al. [28] 

Table §.9 Reduced iron pentacarbonvl and iron cluster HiPco model of Dateo, et al. 
[28] 

Table §.10a. Nickel Cluster Evaporation for n=32-2048 based on Girshick’s 
formulas, Eqn. (7)-(13) 

Table §.10b. Nickel Cluster based on Girshick’s formulas, Eqn. (7)-(13) 

Table §11 Reactions to form nickel/carbon nuclei clusters and nanotubes. 

Table §.12 Reduced carbon vapor/nickel catalyst model 

Table §.13 Comparaison of the arc process with other high temperature processes 
used for SWNT synthesis. 

Table §.14 First ionization potential and energy of some element of the plasma. 

Table §.15 Electronic gas phase reactions considered in Hinkov [116] 

Table §.16 Surface reactions considered in Hinkov [116] 

Table §.17 Plasma characteristics calculated by AURORA assuming LTE and 
Tg=T iO n=Te=6500 K. 

Table §.18 Plasma composition calculated for 1=100 A in helium with P=660 mbar. 
[116] 
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Table §.19 Soot nuclei and soot particles reaction rates. Data from Krestinin and 
Moravskii [119] 

Table §.20 Surface reaction considered in Farhat et al. [110] 

Table §.21 Calculated number densities of small carbon clusters and fullerenes at 
three positions in the arc ( - symbol indicates negligible density). Conditions are 
Helium with P= 660 mbar, 1=100 A, dAc=3 mm, and the dilution factor x=20. 

Table §.22 Simplified Carbon Chemistry Model used in laser ablation simulation 

Table 13 Refs. 
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Figure Captions 


Fig. §.l Comparison of iron cluster evaporation rate coefficients at 1380 K versus 
number of atoms in nickel clusters. 

Fig. §.2 Iron cluster recombination rates at 1500 K. Iron atoms attaching to clusters 
of size n-1 forming clusters of n-atoms. 

Fig. §.3 Nanotube bundles emanating from a nickel/cobalt catalyst particle. 
(Courtesy A. P. Moravsky, MER Corp.) 

Fig. §.4 Nanotube bundles emanating from catalyst particles. (Courtesy A. Loiseau) 
J. Gavillet et al, Carbon 40 (2002) 1649-1663. [41] 

Fig. §.5 Carbon cluster and SWNT production for full model (with modified rates) 
compared with reduced model with temperature ramp from 3500 to 1500 K in 1 ps. 

Fig* §*6 Comparison of cluster evolution and production of SWNTs with full model 
(with modified rates) and the reduced model with temperature ramp 
1500 K in 500 ps. 

Fig. §.7 MER’s industrial reactor of capacity -100 g SWNT soot per hour. On the 
left, it shows soot containing up to 20 wt% of nanotubes. 

Fig. §.8. Schematic of ideal structures (a) and imprities (b) in the collected arc 
materials. 

Fig. §.9. Kratschmer and Huffman arc reactor. 

Fig. §.10 Micrograph of the as produced collaret in the arc chamber; a and b are 
SEM micrographs at two magnifications [110] showing bundles of SWNTs and 
impurities obtained with catalyst composition Ni:Y 4.2:1 atomic percent. C) is a 
high resolution HRTEM of the collaret produced in helium showing bundle of about 
60 nanotubes. [171] 

Fig. §.11 Space charge potential a dielectric field distribution in the arc. 

Fig. §.12 Development of the plasma in the interelectrode zone between the cathode 
on the left and the anode of 6 mm diameter on the right. [116] 

Fig. §.13 Measured optical emission spectra showing the predominance of yttrium 
ions. 7 

Fig. §.14. Calculated nickel and yttrium spectra. [50] 


52 


Fig. §.15 Schematic representation of a well mixed plasma reactor [120] and its 
adaptation to predict plasma compositions Y, for given T and T e . [116] 


Fig. §.16 Predicted variations of ions mass fractions w ith the total current intensity. 
[116] 

Fig. §.17. Plasma characteristics calculated by AURORA in the conditions 
Tg=Tj O n=6500 K and variable T e . Hatched zone corresponds to the optimal domain 
of current intensity 80-120 A for nanotube synthesis. [116] 

Fig. §.18 Time evolution of carbon atom mole fraction calculated in helium, at 
P=660 mbar and T=4000 K. 

Fig. §.19 Time evolution of some carbon cluster mole fraction calculated in helium, 
at P=660 mbar and T=2000 K. 


Fig. §.20 Time evolution of fullerene mole fraction calculated in helium, at P=660 
mbar and T=2000 K. 
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Fig. §.21-a. Mole fractions of small carbon species calculated at equilibrium versus 
temperature calculated in helium at 660 mbar. 


Fig. §.21-b. Total fullerene yield versus temperature calculated in helium at 660 
mbar. 


Fig. §.22 Schematic of the flow issuing from the electric arc zone as represented by 
Krestinin and Moravsky. [111,7,119] 

Fig. §.23 Comparison of experimental and calculated yield of fullerenes versus the 
pressure [8], The experimental molar ratio C 70 /C 60 is found to be invariably equal to 
0 . 2 . 

Fig. §.24 Simulation domain and boundary conditions. 

Fig. §.25 African calabash with structure decorated by hexagons and opened on the 
top. 

Fig. §.26 Gas-Phase species reacting at the edge of an open nanotube and top view of 
growing bundles of nanotubes showing surface site density. 

Fig. §.27 Calculated temperature profile in the inter-electrode gap of the arc. 
Conditions are He, P=660 mbar and 1=100 A. 

Fig. §.28 Calculated axial distribution of small carbon clusters and fullerenes in the 
inter-electrode gap of the arc. Conditions are He, P=660 mbar, 1=100 A, dilution 
factor t=20. 
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Fig. §.29 Calculated normalized species abundance at the cathode in nanotube 
growth in helium, P=660 mbar, 1=100 A. In the right of this figure is superposed a 
comparison with Kroto et al. [1] measurements, a) Low helium density over 
graphite target at time of laser vaporization, b) High helium density over graphite 
target at time of laser vaporization, c) Same as b). but with addition of "integration 
cup" to increase time between vaporization and cluster analysis. 

Fig. §.30 Calculated (10,10) nanotube growth rate versus pre-exponential factor A of 
the reactions (SI) to (S3) in Table §.21 for different dilution factors. Conditions are 
He, P=660 mbar, I=100A. 

Fig. §.31 Two calculation domains in Bilodeau et al. [115] (2D) model: (1) arc region 
delimited by the two opposed movable graphite electrodes and (2) all the reactor. 

Fig. §.32 Temperature field in the interelectrode gap (1) and in all the reactor (2) 
calculated with helium, dAc=l rnm, 1=80 A, erosion rate=9.98 mg/s, deposition 
rate=4,71 mg/s and P=13.3 kPa.From Bilodeau et al. [1 151 
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Fig. §.33 Comparison between argon and helium for optimal nanotube conditions, 
Hinkov.[116] 

Fig. §.34 Laser-ablation oven gridding. 

Fig. §.35 “Primary" grid spacing in revised studies by Greendyke et al of laser 
ablation plumes. 

Fig. §.36 C3 mass fraction contours at 200 ps. 

Fig. §.37 Temperature contours at 200 ps for C3 studies of Greendyke, et al.[142] 

Fig. §.38 Carbon plume propagation in Ar at 100 microseconds post ablation onset. 

Fig. §.39 Carbon plume propagation in N2 at 100 microseconds post ablation onset. 

Fig. §.40 Carbon plume propagation in He at 100 microseconds post ablation onset. 

Fig. §.41 Plume expansion factor for varying background gases. 

Fig. §.42 Plume expansion factor for differing inner quartz tube diameters of the 
laser ablation oven. 

Fig. §.43 Streak lines in space for laser ablation axisymmetric calculation 

Fig. §.44 Temperature histories along selected streak lines for laser ablation 
calculations 
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Fig. §.45 Example of evolution of species mass fractions along streak line in laser 
ablation. Arbitrary evaporation rate coefficients were used in calculation, therefore, 
data is not accurate, but only shown as an example of methodology. 

Fig. §.46 Schematic of a generalized “show erhead” HiPco reactor chamber. 

Fig. §.47 X- and Y-coordinates along four typical streamlines versus the axial 
coordinate Z, in mixing zone of HiPco reactor, calculated by the FLUENT code. 

Fig. §.48 Carbon nanotube cluster distributions at various times in trajectory no. 
injO of HiPco reactor. Published as Fig. 11 in JNN Vol. 3 No. 1/2 pp. 63-73 (2003). 

Fig. §.49 Calculated and measured iron fraction in product of HiPco reactor. 
Published as Fig. 12 in JNN Vol. 3 No. 1/2 pp. 63-73 (2003). 

Fig. §.50 Iron cluster distribution at various times in inlet to HiPco reactor, starting 
from fully dissociated iron pentacarbonyl at 35 atmospheres and 400 K. 
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Table §.l. Gas phase chemistry for Fullerene Model* 



Reaction A(cm Vs/moles) 

P 

E/R(K) 

1 . Chemistry of small clusters (Ci 

C + C 0 C 2 

- Cm) 
2 . 00 x 1 0 14 

0 

0 

C + C2 C 3 

2 . 00 xl 0 14 

0 

0 

C 2 4* C 2 C 3 4- C 

2 . 00 xl 0 15 

0 

9040 

C 2 4- C 2 C 4 

2 . 00 xl 0 14 

0 

0 

Ci + C 3 ^ C 4 

2.00X10 14 

0 

0 

C 1 + C 44-4 C 5 

2.00X10 14 

0 

0 

C 2 4- C 3 <— > C 5 

2-OOxlO 14 

0 

0 

C n -m 4" C m <—> C n 

2 . 00 xl 0 14 

0 

0 

n= 6,10 and m=l,n /2 

2. Chemistry of cycles and polycycles (Cn - C31) 

Cn- m + C m ^C n 2.00X10 14 

0 

0 

n=l 1,31 and m=l,15 
3 . Formation of fullerenes 

Cn-m "I" Cm ^ ) C n 

2.00X10 14 

0 

0 

n=32,46 and n-31 < m < 15 

4 . Growth of fullerenes shells 

C n + C C n +i 

2.00X10 14 

0 

0 

n=32,78 excluding n=59,69 
C n + C 2 — > C n +2 

4.00x1 0 08 

0 

0 

C n +2 ^ Cn + C 2 

3.20X10 13 

0 

61900 

n=32,59 excluding n=58 
C n + C 2 — > C n +2 

4.00x1 0 08 

0 

0 

C n +2 ^ C n + C2 

3.20X10 13 

0 

61900 

n=60,77 excluding n =68 
C n + C 3 <-> C n +2+ c 

l.OOxlO 15 

0 

0 

n=32,77 

5 . Formation and decay of fullerenes molecules 

Cgo — ^ C^of 5.00x10 

0 

37745 

C59 + C — > C^OF 

2.00X10 14 

0 

0 

C58 + C2 — > CftOF 

4.00x1 0 08 

0 

-30196 

C6of C58 + C2 

8 . 00 xl 0 12 

0 

61900 

C58 + C3 — > C( 30 F + C 

8 -OOxlO 14 

0 

0 

C6of + C — » C(>\ 

2.00X10 13 

0 

10065 

C60F + C2 — > C62 

2.00X10 13 

0 

10065 

C<jOF + C 3 — > C(, 3 

2.00X10 13 

0 

10065 

C70 — > C70F 

1 . 20 xl 0 13 

0 

37745 

Cfi9 + C — > C70F 

2 . 00 xl 0 14 

0 

0 

C68 + C2 — > C70F 

4.00x1 0 08 

0 

-30600 

C70F — > C 68 + C2 

2.50X10 14 

0 

0 

C (38 + C3 — » C70F + C 

1.40x10“ 

0 

49925 

C70F + C — > C71 

2.00X10 13 

0 

10065 


C 70 F C 2 — ^ C 72 2.00x10 0 10065 

C 70 F + C 3 -»C 72 +C 2.00xl0 13 0 10065 


* Forward rate constants k are calculated assuming Arrhenius temperature dependence 

£ 

k = AT p exp( ) where A is the pre-exponential factor, (3 is the temperature exponent and E is the 

RT 

activation energy. Constants A, (3 and E are from Krestinin, et al. [7]. 


Table §.2 Thermodynamic Coefficients for Carbon Clusters 


H= 1000 to 20000 K 
L= 300 to 1000 K 


i 

An 

A* 

A 31 

Aji 

A S i 

Afa 

A 71 

1H 

2.3205 176E+00 

1. 82683 12E-04 

-3. 1504542E-08 

2.5371 109E-12 

-6.2267852E-17 

8.485 1595E+04 

5.8203836E+00 

1L 

1 .903 1 864E+00 

4.0200238E-03 

-9.039539 IE-06 

8. 1759092E-09 

-2.5854092E-12 

8.4826660E+04 

7.291 76 17E+00 

2H 

4.0753780E+00 

3.3239978E-04 

- 1.2344909E-08 

-7.07 1404 IE- 13 

2.82223 17E- 17 

9.97881 14E+04 

9.6847543 E-01 

2L 

5.4993641 E+00 

I.8891470E-03 

-1. 3226503 E-05 

1 .6061 423 E-08 

-5.8286773E-12 

9.91 40682 E+04 

-7.7 14607 3 E+00 

3H 

4.9032746E+00 

1.11 13409E-03 

-7.533 1984E-08 

-4. 8099305 E- 13 

8.912 1604E- 17 

9.681 2289 E+04 

1. 50528 16E-01 

3L 

3.60420 14E+00 

7.537371 IE-03 

-1. 1558984E-05 

8.8850242E-09 

-2.52872 16E- 12 

9.701 06 17E+04 

6.11 037 10E+00 

4H 

4.9032746E+00 

1.1 i 13409E-03 

-7.533 1984E-08 

-4.8099305E-13 

8.9121 604E- 1 7 

1 . 1492855E+05 

1. 50528 16E-01 

4L 

3.60420 14E+00 

7.537371 IE-03 

- 1 . 1558984E-05 

8.8850242E-09 

-2.52872 16E- 12 

1. 1512688E+05 

6.11 037 10E+00 

5H 

1.15921 18E+01 

8.550293 IE-04 

- 1.281 7463 E-07 

7.7473555E-12 

- 1. 6260603 E- 16 

1.1368121E+05 

-3.51 15388E+01 

5L 

1. 1731935E+01 

-2.568 1572E-02 

7.8642726E-05 

-7.84835 19E-08 

2.61 16987E-1 1 

1. 149521 1E+05 

-2.9295593E+01 

6 

16 

0 

0 

0 

0 

149572 

-61.28986 

7 

19 

0 

0 

0 

0 

157121 

-77.6175 

8 

22 

0 

0 

0 

0 

178812 

-93.94515 

9 

25 

0 

0 

0 

0 

186965 

-110.2728 

10 

28 

0 

0 

0 

0 

178963 

-126.6004 

11 

31 

0 

0 

0 

0 

181413.9 

-142.9281 

12 

34 

0 

0 

0 

0 

183864.9 

-159.2557 

13 

37 

0 

0 

0 

0 

186315.8 

-175.5834 

14 

40 

0 

0 

0 

0 

188766.8 

-191.911 

15 

43 

0 

0 

0 

0 

191217.7 

-208.2387 

16 

46 

0 

0 

0 

0 

193668.6 

-224.5663 

17 

49 

0 

0 

0 

0 

196119.6 

-240.894 

18 

52 

0 

0 

0 

0 

198570.5 

-257.2216 

19 

55 

0 

0 

0 

0 

201021.5 

-273.5493 

20 

58 

0 

0 

0 

0 

203472.4 

-289.8769 

21 

61 

0 

0 

0 

0 

205923.3 

-306.2045 

22 

64 

0 

0 

0 

0 

208374.3 

-322.5322 

23 

67 

0 

0 

0 

0 

210825.2 

-338.8598 

24 

70 

0 

0 

0 

0 

213276.2 

-355.1874 

25 

73 

0 

0 

0 

0 

215727.1 

-371.5151 

26 

76 

0 

0 

0 

0 

218178 

-387.8427 

27 

1 79 

0 

0 

0 

0 

220629 

-404.1704 

28 

82 

0 

0 

0 

0 

223079.9 

-420.498 

29 

85 

0 

0 

0 

0 

225530.9 

-436.8257 

30 

88 

0 

0 

0 

0 

227981.8 

-453.1533 

31 

91 

0 

0 

0 

0 

230432.7 

-469.481 

32 

94 

0 

0 

0 

0 

232883.7 

-485.8086 

33 

97 

0 

0 

0 

0 

265530.6 

-502.1362 

34 

100 

0 

0 

0 

0 

237785.6 

-518.4639 

35 

103 

0 

0 

0 

0 

270432.5 

-534.7916 

36 

106 

0 

0 

0 

0 

242687.4 

-551.1192 

37 

109 

0 

0 

0 

0 

275334.4 

-567.4468 

38 

112 

0 

0 

0 

0 

247589.3 

-583.7745 

39 

115 

0 

0 

0 

0 

280236.3 

-600.1021 

40 

118 

0 

0 

0 

0 

252491.2 

-616.4297 

41 

121 

0 

0 

0 

0 

285138.1 

-632.7574 

42 

124 

0 

0 

0 

0 

257393.1 

-649.085 

43 

127 

0 

0 

0 

0 

290040 

-665.4127 

44 

130 

0 

0 

0 

0 

262295 

-681.7403 

45 

133 

0 

0 

0 

0 

294941.9 

-698.068 


46 

136 

0 

0 

0 

0 

267196.8 

-714.3956 

47 

139 

0 

0 

0 

0 

299843.8 

-730.7233 

48 

142 

0 

0 

0 

0 

272098.7 

-747.051 

49 

145 

0 

0 

0 

0 

304745.7 

-763.3785 

50 

148 

0 

0 

0 

0 

277000.6 

-779.7062 

51 

151 

0 

0 

0 

0 

309647.5 

-796.0338 

52 

154 

0 

0 

0 

0 

281902.5 

-812.3615 

53 

157 

0 

0 

0 

0 

314549.4 

-828.6891 

54 

160 

0 

0 

0 

0 

286804.4 

-845.0168 

55 

163 

0 

0 

0 

0 

319451.3 

-861.3444 

56 

166 

0 

0 

0 

0 

291706.3 

-877.6721 

57 

169 

0 

0 

0 

0 

324353.2 

-893.9997 

58 

172 

0 

0 

0 

0 

296608.1 

-910.3274 

59 

175 

0 

0 

0 

0 

329255.1 

-926.655 

60H 

9.98434 18E+01 

7.8857558E-02 

-3.0608799E-05 

5. 1957690E-09 

-3.2188408E-13 

2.6670488E+05 

-5.4587488E+02 

60L 

-3.3579084E+01 

4.2844440E-01 

-3.171 232 1 E-04 

4.7546257E-08 

2.7677699E- 1 1 

3.0465 122E+05 

1.4832875E+02 

61 

181 

0 

0 

0 

0 

334156.9 

-959.3103 

62 

184 

0 

0 

0 

0 

30641 1.9 

-975.638 

63 

187 

0 

0 

0 

0 

339058.8 

-991.9655 

64 

190 

0 

0 

0 

0 

311313.8 

-1008.293 

65 

193 

0 

0 

0 

0 

343960.7 

-1024.621 

66 

196 

0 

0 

0 

0 

316215.6 

-1040.949 

67 

199 

0 

0 

0 

0 

348862.6 

-1057.276 

68 

202 

0 

0 

0 

0 

321117.5 

-1073.604 

69 

205 

0 

0 

0 

0 

353764.4 

-1089.931 

70H 

1.0677602E+02 

1.0200334E-01 

-3.9562455E-05 

6.7122022E-09 

-4.1568097E-13 

2.9769320E+05 

-5.994 1734E+02 

70L 

-3.1 366983 E+01 

3.34844 10E-01 

9.8587783E-05 

-4.4663983E-07 

2.2008836E-10 

3.4260294E+05 

1.4891 61 4E+02 

71 

211 

0 

0 

0 

0 

358666.3 

-1 122.587 

72 

214 

0 

0 

0 

0 

330921.3 

-1138.914 

73 

217 

0 

0 

0 

0 

363568.2 

-1155.242 

74 

220 

\o~ 

0 


0 

335823.2 

-1171.57 

75 

223 

0 

0 

0 

0 

368470. 1 

-1187.897 

76 

226 

0 

0 

0 

0 

340725 

-1204.225 

77 

229 

0 

0 

0 

0 

373372 

-1220.553 

78 

232 

0 

0 

0 

0 

345626.9 

-1236.88 

79 

235 

0 

0 

0 

0 

378273.9 

-1253.208 

80 

238 

0 

0 

0 

0 

350528.8 

-1269.536 


Table §.3 - Reduced Fullerene model reaction scheme 


A r| E/R 

Cluster Growth 


c 

+ c 

= 

C 2 


2 . OE+14 

0. 

0. 

c 

+ c 2 

= 

C 3 


2 . OE+14 

0. 

0. 

C2 

+ C 2 

= 

C 3 + 

C 

2 . OE+15 

0. 

9040. 

C3 

+ c 

= 

0.100CC 


2 . OE+13 

0. 

0. 

C3 

+ c 2 

= 

0 . 125CC 


2 . OE+13 

0. 

0.0 

C3 

+ C 3 

= 

0 . 150CC 


2 . OE+13 

0. 

0.0 

CC 

+ c 

= 

1.025CC 


2 . 3E+14 

0. 

0.0 

CC 

+ C 2 

=> 

1 . 05CC 


2 . 3E+14 

0 . 

0.0 

CC 


— > 

0.95CC + 

c 2 

3 -2E+13 

0. 

61900. 

Fullerene 

formation 





CC 

+ c 3 

= > 

0 • 7 OC 60p + C 


2 . OE+13 

0. 

0. 

CC 

+ c 2 

=> 

0-70C 60F 


4 . OE+9 

0. 

-30196. 

CC 

+ c 

=> 

0 . 6833333C, np 


2 .OE+13 

0. 

0. 

Fullerene 

disintegration 





C 60F 


= > 

1.45CC + 

C2 

8. OE+13 

0. 

61900. 

Soot 

formation 






CC 

+ CC 

= > 

Z 


4. OE+13 

0. 

0. 

CC 

+ C 2 

= > 

0.525Z 


4 . OE+8 

0. 

-30600. 

CC 

+ C 3 

=> 

0.5375Z 


6.92E+10 

0.5 

0. 

C 60F + 

Z 

=> 

1.75Z 


1.26E+12 

0. 

10098. 

CC 

+ c 3 

= > 

0.525Z + 

C 

2 . OE+11 

0. 

0. 

z + c 3 


=> 

1.0375Z 


4 . OE+12 

0. 

0. 

z + c 2 


= > 

1.025Z 


4 . OE+12 

0. 

0. 

z + c 


=> 

1. 0125Z 


4 .OE+12 

0. 

0. 


Table §.4 Thermodynamic Property NASA Coefficients Reduced Fullerene Model plus some other 
species used in plasma simulations 


AR IVTAN AR 1 G 300.000 20000.000 1000.00 

0 .25199477E+01-0 .228583 82E-04 0 . 74655237E-08-0 . 902 082 51E-12 0 . 360 6373 3E-16 


-0. 75393340 E+ 03 0 . 42492779E+01 0 . 256702 00E+01-0 . 44283132E-03 0 . 98495011E-06 
-0 .88424193E-09 0 . 27879189E-12-0 . 75176371E+03 0 . 40801018E+01 
C GMcB 20K 0C 1 G 300.000 20000.000 1000.00 

0 .23205176E+01 0 . 18268312E-03-0 . 31504542E-07 0 . 25371109E-11-0 . 62267852E-16 
0 .84851595E+05 0 . 582 03 83 6E+01 0 . 19031864E+01 0 . 40200238E-02-0 . 90395391E-05 
0. 81759092 E- 08-0. 25854092E-11 0 . 84826660E+05 0 . 72 917617E+01 
C2 GMcB 20K 0C 2 G 300.000 20000.000 1000.00 

0 . 40753780E+01 0 . 33239978E-03-0 . 12344909E-07-0 . 70714041E-12 0 . 28222317E-16 
0 .99788114E+05 0 . 96847543E+00 0 . 54993 641E+01 0 . 18891470E-02-0 . 13226503E-04 
0 .16061423E-07-0 . 58286773E-11 0 . 99140682E+05-0 . 77 146073E+01 
C3 GMcB 20K 0C 3 G 300.000 20000.000 1000.00 

0 . 49032746E+01 0 . 11113409E-02-0 . 75331984E-07-0 . 48099305E-12 0 . 89121604E-16 
0 . 96812289E+05 0 . 15052816E+00 0 . 3 6042 014E+01 0 . 753737 11E-02 -0 . 11558984E-04 


0 . 88850242E-08- 
C60F 

0 . 1780000E+03 
0 . 3 015100E+06 
O.OOOOOOOE+OO 
CC 

0.1180000E+03 
0 . 2524912E+06 
0.0000000E+00 
Z 

0 . 2380000E+03 
0 . 3505288E+06 
O.OOOOOOOE+OO 


0 . 25287216E-11 
C 60 

0.0000000E+00 
-0 . 9429827E+03 
0.0000000E+00 
C 40 

0.0000000E+00 
-0. 6164297 E+ 03 
0.0000000E+00 
C 80 

0.0000000E+00 
-0 . 12 6953 6E+04 
0.0000000E+00 


0 . 97010617E+05 
0 0 G 

0.0000000E+00 
0 . 1780000E+03 
0 . 3 015100E+06 
G 

0.0000000E+00 
0 . 1180000E+03 
0 . 2524912E+06 
0 0 G 

O.OOOOOOOE+OO 
0 . 2380000E+03 
0 .3505288E+06 


0 . 61103710E+01 
0300.00 5000 

O.OOOOOOOE+OO 
O.OOOOOOOE+OO 
-0 . 9429827E+03 
300.000 5000 

O.OOOOOOOE+OO 
O.OOOOOOOE+OO 
-0 . 6164297E+03 
0300.00 5000 

O.OOOOOOOE+OO 
O.OOOOOOOE+OO 
-0 . 1269536E+04 


00 1000.00 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

000 1000.00 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

00 1000.00 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 


1 

2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 


Plasma Species 

Ni GMcB 2 OK ONi 1 G 300.000 20000.000 1000.00 1 

0. 36971152 E+ 01-0. 82867716E-03 0 . 18628604E-06-0 . 1198693 0E-10 0 . 23752778E-15 2 

0 . 50060051E+05 0 . 72075154E+00 0 . 42140643E+01-0 . 10224755E-01 0 . 25272825E-04 3 

-0. 24064628 E- 07 0 . 78454681E-11 0 . 50323 075E+05 0 . 58841215E-01 4 

Ni+ GMcB 20K ONi IE -1 G 300.000 20000.000 1000.00 1 

0 . 21013259E+01 0 . 71871149E-03 -0 . 11060734E-06 0 . 67741068E-11-0 . 13425799E-15 2 

0 . 14040022 E+ 06 0 . 89762594E+01 0 . 122 60499E+00 0 . 14956310E-01-0 . 30951078E-04 3 

0. 26889803 E- 07-0. 830157 04E-11 0 . 14047983E+06 0 . 16941765E+02 4 

Y IV 25K OY 1 G 300.000 25000.000 1000.00 1 

0 . 79437350 E+ 00 0 . 13579116E-02-0 . 12031792E-06 0 . 31899976E-11-0 . 15582 902E-16 2 

0 . 51410499E+05 0 . 18158732E+02-0 . 66965149E+01 0 . 65805915E-01-0 . 14851900E-03 3 

0 .13342009E-06-0 . 41975353E-10 0 . 51184455E+05 0 . 45676826E+02 4 

Y+ OY IE -1 G 300.000 25000.000 1000.00 1 

0 .30001765E+01 0 . 193 53 870E-04 0 . 10693 682E-07-0 . 93 5773 06E-12 0 . 2 0474488E-16 2 

0 . 12272102E+06 0 . 462 80979E+01-0 . 10156546E+02 0 . 946277 61E-01-0 . 20487420E-03 3 

0. 17857051 E- 06-0. 55138220E-10 0 . 12325328E+06 0 . 57606297E+02 4 


C+ GMcB 2 OK OC IE -1 G 300.000 20000.000 1000.00 1 

0 . 25063181E+01 0 . 163 90834E-05-0 . 41024961E-08 0 . 78583063E-12-0 . 23326192E-16 2 

0.21547022E+06 0 . 43269073E+01 0 . 2 6069164E+01-0 . 48079487E-03 0 . 83960952E-06 3 

-0. 64172956 E- 09 0 . 18061570E-12 0 . 21545410E+06 0 . 38615888E+01 4 

C4 + 121286C 4E ^1 G 0300.00 5000.00 1000.00 1 

0 . 06500180E+02 0 . 04228632E-01-0 . 01790717E-04 0 . 03404812E-08-0 . 02403 978E-12 2 

0. 24662108 E+ 06-0. 11488894E+ 02 0 . 02343 028E+02 0 . 1642 9811E-01-0 . 15279858E-04 3 

0. 07343826 E- 07-0. 15822743E-11 0 . 24662 184E+06 0 . 0982 62 04E+02 4 

C5+ 121286C 5E -1 G 0300.00 5000.00 1000.00 1 

0 . 08078081E+02 0 . 05743464E-01-0 . 02436405E-04 0 . 04638916E-08-0 . 03278909E-12 2 


0. 23886116 E+ 06-0. 01953023 E+ 03 0 . 02115273E+02 0 . 02326331E+00-0 . 02109499E-03 
0.09072734E-07-0.15400926E-11 0 . 23886181E+06 0 . 10976027E+02 


3 

4 


C4 GMcB 2 OK OC 4 G 300.000 20000.000 1000.00 1 

0 .49032746E+01 0 . 11113409E-02-0 . 7533 1984E-07-0 . 48099305E-12 0 . 89121604E-16 2 

0 . 11492855E+06 0 . 15052816E+00 0 . 36042014E+01 0 . 75373711E-02 -0 . 11558984E-04 3 

0. 88850242 E- 08-0. 25287216E-11 0 . 11512688E+06 0 . 61103710E+01 4 

C5 GMcB 20K OC 5 G 300.000 20000.000 1000.00 1 

0 . 11592118E+02 0 . 85502931E-03-0 . 12817463E-06 0 . 77473555E-11-0 . 16260603E-15 2 

0.11368121E+ 06-0. 35115388 E+ 02 0 . 11731935E+02-0 . 25681572E-01 0 . 78642726E-04 3 

-0.78483519E-07 0 . 26116987E-10 0 . 11495211E+ 06-0 . 29295593E+02 4 

HE 120186HE 1 G 0300.00 5000.00 1000.00 1 

0 . 02500000E+02 0 . OOOOOOOOE+OO 0 . OOOOOOOOE+OO 0 . OOOOOOOOE+OO 0 . OOOOOOOOE+OO 2 

-0.074537 50E+ 04 0 . 09153489E+01 0 . 02500000E+02 0 . OOOOOOOOE+OO 0 . OOOOOOOOE+OO 3 

0. OOOOOOOOE+OO 0. 00000000 E+ 00-0. 07453750 E+ 04 0 . 09153488E+01 4 


Table §.5 Iron Carbonyl and Iron Cluster Model of Krestinin, Smirnov, and 
Zaslonko [29] 


1 . Iron Carbonyl Reactions (Krestinin, et 
al.) 

A 

n 

Ea/R 

Fe(CO)5 => FeCO + 4CO 

2.0E+15 

0.0 

20130.9 

FeCO => Fe + CO 

6.0E+14 

0.0 

10317.1 

|Fe + CO + M => FeCO +M 

5.0E+14 

0.0 

0.0 


2. Fe Nucleation 

FeCO + FeCO => Fe2 + CO + CO 

6.0E+14 

0.0 

1006.5 

FeCO + Fe => Fe2 + CO 

6.0E+14 

0.0 

0.0 

1 Fe + Fe + M => Fe2 + M 

6.0E+14 

0.0 

0.0 


3a. Formation small clusters 
Fe2 + FeCO => Fe3 + CO 

6.0E+14 

0.0 

0.0 

Fe3 + FeCO => Fe4 + CO 

6.0E+14 

0.0 

0.0 

Fe4 + FeCO => Fe5 + CO 

6.0E+14 

0.0 

0.0 

Fe2 + Fe + M => Fe3 +M 

5.0E+15 

0.0 

0.0 

Fe3 + Fe => Fe4 

5.0E+11 

0.0 

0.0 

Fe4 + Fe => Fe5 

3.0E+13 

0.0 

0.0 


3b. Evaporation of small clusters 


Fe2 + M => Fe + Fe + M 

1 El 5 ( 1 El 6) 

0.0 

16104.7 

Fe3 + M => Fe2 + Fe + M 

1 E16 ( 1 El 7) 

0.0 

22143.9 

Fe4 => Fe3 + Fe 

1 E13 (2E13) 

0.0 

27679.9 

Fe5 => Fe4 + Fe 

2E14 (7E14) 

0.0 

30699.5 

Fe6 => Fe5 + Fe 

2E1 5 (4E15) 

0.0 

32712.6 

Fe7 => Fe6 + Fe 

5E15 (1 E16) 

0.0 

33719.2 

Fe8 => Fe7 + Fe 

1 E16 (2E16) 

0.0 

34222.4 

Fe9 => Fe8 + Fe 

4E16 (4E16) 

0.0 

35229.0 

3. Fe-atom Exchange 

Fe2 + Fe2 => Fe3 + Fe 

3.0E+14 

0.0 

0.0 

Fe2 + Fe3 => Fe4 + Fe 

3.0E+14 

0.0 

0.0 

4. Coagulation of small clusters 

Fe2 + Fe2 => Fe4 

5.0E+13 

0.0 

0.0 

Fe2 + Fe3 => Fe5 

1.0E+14 

0.0 

0.0 

Fe2 + Fe4 => Fe6 

3.0E+14 

0.0 

0.0 

Fe3 + Fe3 => Fe6 

3.0E+14 

0.0 

0.0 

Fe3 + Fe4 => Fe7 

3.0E+14 

0.0 

0.0 

Fe4 + Fe4 => Fe8 

3.0E+14 

0.0 

0.0 

5. Heterogeneous reaction in the c-phase (n 
Fe + Fen => Fe(n+1) 

>= 5 ) 

2.3E+13 

0.0 

0.0 

FeCO + Fen => Fe(n+1 ) + CO 

2.3E+13 

0.0 

0.0 

Fe2 + Fen => Fe(n+2) 

2.3E+13 

0.0 

0.0 

Fe3 + Fen => Fe(n+3) 

2.3E+13 

0.0 

0.0 

Fe4 + Fen => Fe(n+4) 

2.3E+13 

0.0 

0.0 


Fen + Fem => Fe(n+m) 


0.0 


0.0 


knm = Zn ((n 1,3 + m 1/3 )/ 2 ) 2 ((n+m)/2nm) 1/2 


Table §.6 Rate coefficients for iron pentacarbonyl decomposition and recombination 

High Pressure Ames ;l High Pressure Smirnov b Low Pressure (fall off) Smirnov* 1 


Reaction A(mol,cm 3 , s) t) E/R (K) A(mol,cm'\ s) r| E/R (K) A(mol,cm 3 , s) r| E/R (K) 


Fe(CO) 5 Fe(CO) 4 + CO 

4.62X1 0 2U 

-1.96 

21028.3 

4.26X10 20 

-1.85 

20447 . 

1.22X10 73 

-16.78 

20447 . 

Fe(CO) 4 + CO — » Fe(CO) 5 c 

3.50X10 IU 

0. 

0. 







Fe(CO) 4 — > Fe(CO) 3 + CO 

6.96X10“ 

-2.14 

14346.3 

2.85X10 21 

-1.87 

16598. 

2.89X10 57 

-12.41 

16598 . 

Fe(CO) 3 + CO -» Fe(CO) 4 c 

1.30X10 13 

0. 

.0 







Fe(CO) 3 — > Fe(CO )2 + CO 

8.70X10 19 

-1.58 

14679.8 

2 . 09X10 21 

-1.85 

17921. 

2 . 10X10 47 

-9.11 

17921 . 

Fe(CO) 2 + CO -> Fe(CO) 3 c 

1.80X10 1 ' 

0. 

0. 







Fe(CO) 2 -» FeCO + CO 

3.96X10“ 

-2.29 

18799.4 

4.07X10 19 

-1.46 

13110. 

7 . 92X10 33 

-5.4 

13110. 

FeCO + CO -» Fe(CO) 2 c 

1.50X10 13 

0. 

0. 







FeCO + M Fe + CO + M 

4.87X10 1 '-’ 

-0.90 

2874.5 

1. 073X10 23 

-2.6 

6495. 

1.073X10 23 

-2.6 

6495. 

Fe + CO + M — > FeCO + M d 

1.00X10 15 

0. 

0. 








a Ref. 1 28] 
b Ref. [32] 
‘Ref. [35] 
d Ref. [34] 


Table §.7 Fits to Vlasov Fig. 2 Recombination Rates 


Fe3 

+ Fe => 

Fe4 

1.28E+10 

0.5 

0 

Fe4 

+ Fe => 

Fe5 

1.70E+10 

0.5 

0 

Fe5 

+ Fe => 

Fe6 

2.24E+10 

0.5 

0 

Fe6 

+ Fe => 

Fe7 

2.89E+10 

0.5 

0 

Fe7 

+ Fe => 

Fe8 

1.73E+1 1 

0.5 

0 

Fe8 

+ Fe => 

Fe9 

2.53E+1 1 

0.5 

0 

Fe9 

+ Fe => 

FelO 

3.18E+1 1 

0.5 

0 

FelO 

+ Fe => 

Fell 

3.90E+1 1 

0.5 

0 

Fel 1 

+ Fe => 

Fel 2 

4.71E+1 1 

0.5 

0 

Fel2 

+ Fe => 

Fel 3 

5.60E+1 1 

0.5 

0 

Fel 3 

+ Fe => 

Fel4 

6.57E+1 1 

0.5 

0 

Fel 4 

+ Fe => 

Fel 5 

7.63E+1 1 

0.5 

0 

Fel5 

+ Fe => 

Fel 6 

8.77E+1 1 

0.5 

0 

Fel6 

+ Fe => 

Fel 7 

1.00E+12 

0.5 

0 

Fel7 

+ Fe => 

Fel 8 

1.13E+12 

0.5 

0 

Fel 8 

+ Fe => 

Fel9 

1.27E+12 

0.5 

0 

Fel 9 

+ Fe => 

Fe20 

1.42E+12 

0.5 

0 

Fe20 

+ Fe => 

Fe21 

1.58E+12 

0.5 

0 

Fe21 

+ Fe => 

Fe22 

1.74E+12 

0.5 

0 

Fe22 

+ Fe => 

Fe23 

1.92E+12 

0.5 

0 

Fe23 

+ Fe => 

Fe24 

2.11E+12 

0.5 

0 

Fe24 

+ Fe => 

Fe25 

2.30E+12 

0.5 

0 

Fe25 

+ Fe => 

Fe26 

2.50E+12 

0.5 

0 

Fe26 

+ Fe => 

Fe27 

2.72E+12 

0.5 

0 

Fe27 

+ Fe => 

Fe28 

2.94E+12 

0.5 

0 

Fe28 

+ Fe => 

Fe29 

3.17E+12 

0.5 

0 

Fe29 

+ Fe => 

Fe30 

3.41E+12 

0.5 

0 

Fe30 

+ Fe => 

Fe31 

3.66E+12 

0.5 

0 

Fe31 

+ Fe => 

Fe32 

3.92E+12 

0.5 

0 

Fe32 

+ Fe => 

Fe33 

4.19E+12 

0.5 

0 

Fe33 

+ Fe => 

Fe34 

4.47E+12 

0.5 

0 

Fe34 

+ Fe => 

Fe35 

4.76E+12 

0.5 

0 

Fe35 

+ Fe => 

Fe36 

5.06E+12 

0.5 

0 

Fe36 

+ Fe => 

Fe37 

5.37E+12 

0.5 

0 

Fe37 

+ Fe => 

Fe38 

5.69E+12 

0.5 

0 

Fe38 

+ Fe => 

Fe39 

6.01E+12 

0.5 

0 

Fe39 

+ Fe => 

Fe40 

6.35E+12 

0.5 

0 

Fe40 

+ Fe => 

Fe41 

6.70E+12 

0.5 

0 

Fe41 

+ Fe => 

Fe42 

7.06E+12 

0.5 

0 

Fe42 

+ Fe => 

Fe43 

7.43E+12 

0.5 

0 

Fe43 

+ Fe => 

Fe44 

7.81E+12 

0.5 

0 

Fe44 

+ Fe => 

Fe45 

8.19E+12 

0.5 

0 

Fe45 

+ Fe => 

Fe46 

8.59E+12 

0.5 

0 

Fe46 

+ Fe => 

Fe47 

9.00E+12 

0.5 

0 


Fe47 

+ Fe 

=> 

Fe48 

9.42E+12 

0.5 

0 

Fe48 

+ Fe 

=> 

Fe49 

9.85E+12 

0.5 

0 

Fe49 

+ Fe 

=> 

Fe50 

1.03E+13 

0.5 

0 

Fe50 

+ Fe 

=> 

Fe51 

1.07E+13 

0.5 

0 

! Fits 
Fe4 

to Vlasov Evaporation Rates 
=> Fe3 + Fe 4.255E+13 

0 

24660 

Fe5 

=> 

Fe4 

+ Fe 

4.323E+13 

0 

24720 

Fe6 

=> 

Fe5 

+ Fe 

4.574E+13 

0 

36140 

Fe7 

=> 

Fe6 

+ Fe 

5.287E+13 

0 

34380 

Fe8 

=> 

Fe7 

+ Fe 

5.196E+13 

0 

25020 

Fe9 

=> 

Fe8 

+ Fe 

5.562E+13 

0 

31990 

FelO 

=> 

Fe9 

+ Fe 

5.859E+13 

0 

31250 

Fell 

=> 

FelO 

+ Fe 

7.06E+13 

0 

35970 

Fel2 

=> 

Fell 

+ Fe 

5.095E+13 

0 

38700 

Fel3 

=> 

Fel2 

+ Fe 

5.875E+13 

0 

47940 

Fel4 

=> 

Fel3 

+ Fe 

6.321E+13 

0 

32840 

Fel5 

=> 

Fel4 

+ Fe 

6.45E+13 

0 

41610 

Fel6 

=> 

Fel5 

+ Fe 

6.178E+13 

0 

37050 

Fel7 

=> 

Fel6 

+ Fe 

7.468E+13 

0 

35380 

Fel8 

=> 

Fel7 

+ Fe 

6.98E+13 

0 

34180 

Fel9 

=> 

Fel8 

+ Fe 

7.509E+13 

0 

40220 

Fe20 

=> 

Fel9 

+ Fe 

7.076E+13 

0 

38380 

Fe21 

=> 

Fe20 

+ Fe 

7.40298E+13 

0 

38617.6229 

Fe22 

=> 

Fe21 

+ Fe 

7.592 16E+ 13 

0 

38769.36455 

Fe23 

=> 

Fe22 

+ Fe 

7.78134E+13 

0 

38914.9169 

Fe24 

=> 

Fe23 

+ Fe 

7.97052E+13 

0 

39054.78564 

Fe25 

=> 

Fe24 

+ Fe 

8.1 109E+13 

0 

39206 

Fe26 

=> 

Fe25 

+ Fe 

8.34888E+13 

0 

39319.20281 

Fe27 

=> 

Fe26 

+ Fe 

8.53806E+13 

0 

39444.49644 

Fe28 

=> 

Fe27 

+ Fe 

8.72724E+13 

0 

39565.61064 

Fe29 

=> 

Fe28 

+ Fe 

8.91642E+13 

0 

39682.82685 

Fe30 

=> 

Fe29 

+ Fe 

9. 1372E+13 

0 

39856 

Fe31 

=> 

Fe30 

+ Fe 

9.29478E+13 

0 

39906.55531 

Fe32 

=> 

Fe31 

+ Fe 

9.48396E+13 

0 

40013.50488 

Fe33 

=> 

Fe32 

+ Fe 

9.673 14E+13 

0 

40117.43669 

Fe34 

=> 

Fe33 

+ Fe 

9.86232E+13 

0 

40218.52357 

Fe35 

=> 

Fe34 

+ Fe 

1.0377E+14 

0 

40423 

Fe36 

=> 

Fe35 

+ Fe 

1.02407E+14 

0 

40412.78232 

Fe37 

=> 

Fe36 

+ Fe 

1.04299E+14 

0 

40506.23296 

Fe38 

=> 

Fe37 

+ Fe 

1.0619E+14 

0 

40597.39881 

Fe39 

=> 

Fe38 

+ Fe 

1.08082E+14 

0 

40686.39369 

Fe40 

=> 

Fe39 

+ Fe 

1.09974E+14 

0 

40773.32304 

Fe41 

=> 

Fe40 

+ Fe 

1.1 1866E+14 

0 

40858.28466 

Fe42 

=> 

Fe41 

+ Fe 

1. 13758E+14 

0 

40941.36951 

Fe43 

=> 

Fe42 

+ Fe 

1.15649E+14 

0 

41022.66227 

Fe44 

=> 

Fe43 

+ Fe 

1.17541E+14 

0 

41102.24195 


Fe45 

=> 

Fe44 

+ Fe 

1.19433E+14 

0 

41180.18238 

Fe46 

=> 

Fe45 

+ Fe 

1.21325E+14 

0 

41256.55265 

Fe47 

=> 

Fe46 

+ Fe 

1.23217E+14 

0 

41331.41751 

Fe48 

=> 

Fe47 

+ Fe 

L.25108E+14 

0 

41404.83774 

Fe49 

=> 

Fe48 

+ Fe 

1.2529E+14 

0 

41377 

Fe50 

=> 

Fe49 

+ Fe 

1.28892E+14 

0 

41547.56942 


Table §.8 Full iron pentacarbonyl and iron cluster HiPco model of Dateo, et al. [28] 


Thermal Decomposition and Formation of 
Reaction 


Small Iron Clusters from FeCO Fragment 
A (modern 3 , s) E a /k (K) 


1 . Ames 

Carbonyl Rates 






Fe (CO ) 5 



=> 

Fe (CO ) 4 

+ 

CO 

4 . 620E+20 

-1.96 

21028.30 

Fe (CO ) 4 

+ 

CO 

= > 

Fe (CO ) 5 



3 . 500E+10 

0.00 

0.00 

Fe (CO ) 4 



=> 

Fe (CO ) 3 

+ 

CO 

6 . 960E+22 

-2.14 

14346.30 

Fe (CO ) 3 

+ 

CO 

= > 

Fe (CO) a 



1 . 300E+13 

0.00 

0.00 

Fe (CO ) 3 



=> 

Fe (CO ) 2 

+ 

CO 

8 . 700E+19 

-1.58 

14679.80 

Fe (CO ) 2 

+ 

CO 

=> 

Fe (CO ) 3 



1.800E+13 

0.00 

0.00 

Fe (CO ) 2 



=> 

FeCO 

+ 

CO 

3 . 960E+21 

-2.29 

18799.40 

FeCO 

+ 

CO 

=> 

Fe (CO ) 2 



1 .500E+13 

0.00 

0.00 

FeCO 

+ 

M 

=> 

Fe 

+ 

CO + M 

4 . 870E+19 

-0.90 

2874.50 

Fe + CO 

i + 

M 

=> 

FeCO + 

M 


1 . 000E+15 

0 

0 


Krestinin 


2 Formation and Growth of Small Iron Clusters 
Fe + Fe + M => Fe 2 + M 
Fe 2 + Fe( + M) => Fe 3 (+ M) 

Krestinin Recombination Rates 


4 . 00E+15 0 0 
5.00E+15 0 0 


Fe 3 

+ 

Fe => 

Fe 4 






5.10E+11 0 

0 . 00E+00 

Tl 

(V 

+ 

Fe => 

Fe 5 






3.00E+13 0 

0 . 00E+00 

Fe s 

+ 

Fe => 

Fe 6 






1.10E+12 0 

0.00E+00 

Fe 6 

+ 

Fe => 

Fe 7 






1.50E+12 0 

0.00E+00 

Fe 7 

+ 

Fe => 

Fe 8 






7.90E+12 0 

0.00E+00 

Fe 8 

+ 

Fe => 

F eg 






1.00E+13 0 

0 . 00E+00 

F eg 

+ 

Fe => 

Fe 10 






1.30E+13 0 

0 . 00E+00 

3. 

Thermal Decomposition 

of Small 

Iron Clusters 


Krestinin 









Fe 2 


+ M 

=> 

Fe 

+ 

Fe 

+ M 

1. 00E+15 

0 16102.93 


Fe 3 


+ M 

=> 

Fe 2 

+ 

Fe 

+ M 

1.00E+16 

0 22141.54 


Fe4 



=> 

Fe 3 

+ 

Fe 


1.00E+13 

0 27676.92 


Fe5 



=> 

Fe4 


Fe 


2 . 00E+14 

0 30696.22 


Fe 6 



=> 

Fe5 

+ 

Fe 


2.00E+15 

0 32709.09 


Fe7 



=> 

Fe 6 

+ 

Fe 


5 . 00E+15 

0 33715.52 


Fe 8 



=> 

Fe7 

+ 

Fe 


1 . 00E+16 

0 34218.74 


Fe9 



=> 

Fe 8 

+ 

Fe 


4.00E+16 

0 35225.17 


4. i 

Coagulation of 

Small 

Iron Clusters 



Fe 2 

+ 

Fe 2 => 

Fe4 





5.00E+13 

0 0 


Fe 2 

+ 

Fe 3 => 

Fe5 





1 . 00E+14 

0 0 


Fe 2 

+ 

Fe 2 => 

Fe 3 + 

Fe 




3.00E+14 

0 0 


Fe 2 

+ 

Fe 3 => 

Fe4 + 

Fe 




3.00E+14 

0 0 


Fe 2 

+ 

F e 4 => 

Fe 6 





3.00E+14 

0 0 


Fe 3 

+ 

Fe 3 => 

Fe 6 





3.00E+14 

0 0 


Fe 3 

+ 

F e 4 — > 

Fe 7 





3 . 00E+14 

0 0 


Fe4 

+ 

F e 4 => 

Fe 8 





3 . 00E+14 

0 0 


Fe 2 

+ 

Fe 5 => 

Fe 7 





2.30E+13 

0 0 


Fe 2 

+ 

Fe 6 => 

Fe 8 





2.30E+13 

0 0 



5 . Coagulation of Larger Clusters 


Fe n 

+ 

Fe => 

Fe n+1 

2.30E+13 

0 

0 . 00E+00 

9<n<200 

Fe n 

4 - 

Fe 2 => 

Fe n+2 

2.30E+13 

0 

0.00E+00 

9<n<199 

Fe n 

+ 

Fe 3 => 

Fe n+3 

2 . 30E+13 

0 

0 . 00E+00 

9<n<198 


Fe n + Fe 4 => Fe n+4 


2.30E+13 0 O.OOE+OO 9<n<197 


6. Attachment of CO to Iron Cluters 

Fe n + CO => Fe n CO 1.00E+13 0 0 9<n<201 

7. Desorption of CO from Clusters 

Fe n CO => Fe n + CO 1.00E+15 0 16000 9<n<201 

8. Formation of SWNTs 

Fe n C0 + CO => PFenC OT + (1+3) /2 C0 2 +(l-P)Fe n CO 1.00E+16 0 12500 

9. End of growth of nanotubes (overcoating) 

Fe n CNT => DFe n CNT 1.00E+00 0 0 9<n<201 

10. End of Catalyst Effectiveness 

Fe n => DFe n 1.00E+03 0 0 9<n<201 

P = 1/ ^nNT-l) 0,^=3000 in Dateo, et al.[28] 


9<n<201 


Table §.9 Reduced iron pentacarbonyl and iron cluster HiPco model of Dateo, et al. [28] 


Thermal Decomposition and Formation of Smal 1 Iron Clusters from FeCO Fragment 
Ames Carbonyl Rates 


Fe(CO) 5 => Fe(CO) 4 +CO 

4.620E+20 

-1.96 

21028.30 

Fe(CO) 4 +CO => Fe(CO) 5 

3.500E+10 

0.00 

0.00 

Fe(CO) 4 => Fe(CO)3 +CO 

6.960E+22 

-2.14 

14346.30 

Fe(CO) 3 +CO => Fe(CO) 4 

1.300E+13 

0.00 

0.00 

Fe(CO) 3 => Fe(CO)2 + CO 

8.700E+19 

-1.58 

14679.80 

Fe(CO) 2 +CO => Fe(CO) 3 

1.800E+13 

0.00 

0.00 

Fe(CO) 2 => FeCO + CO 

3.960E+21 

-2.29 

18799.40 

FeCO + CO => Fe(CO) 2 

1.500E+13 

0.00 

0.00 

FeCO + M => Fe + M + CO 

4.870E+19 

-0.90 

2874.50 

Fe + CO + M => FeCO + M 

1.00E+15 

0 


2. Growth and Evaporation of Iron Clusters 

Fe + Fe + M => Fe2 + M 4.00E+ 1 5 

0 

0 

Fe 2 + Fe + M => 3aFec + M 

5.00E+15 

0 

0 

Fe2 + M => Fe + Fe + M 

1.00E+15 

0 

16102.93 

Fe 2 + Fe 2 => 4aFec 

5.00E+13 

0 

0 

Fe 2 + Fe => Fes 1 -00E+ 14 0 0 
Fe 2 + Fee => ( l+2a)Fec 

3.00E+14 

0 

0 

3. Attachment of CO on Iron Clusters 

Fee + CO => FecCO 

1.00E+13 

0 

0 

4. Desorption of CO from Clusters 

FecCO — > Fee + CO 

1.00E+15 

0 

16000 


5. Formation of SWNTs 

Fe c CO +CO => pFecC N T + (l+P)/2C0 2 + (l-p)FecC0 1.00E+16 0 12500 

6. End of growth of nanotubes 

FecCNT + Fee => P DFecCNT + (1- P)FccCnt 3.5E14 0 0 

7. End of Catalyst Effectiveness 

Fee => oc DFec 1000. 0 0 


5 5 

P = l/(2nNT-l), ot = 1/nc, oc — nc/nc 
P " = n c /(n c *-n c ), n NT = 3000, n c = 20, n c = 40 


Table §.10a. Nickel Cluster Evaporation for n=32-2048 based on Girshick’s formulas, Eqns. (§.7)-( 
§.13) 



Reaction 

A (cm 7s/moles) 

e 

E/R (K) 

Ni 32 

— *• Ni|6 + 16Ni 

9.38X10 17 

-0.5 

39127. 

Ni 64 

— * Ni 32 + 32Ni 

1.34X10 18 

-0.5 

40843. 

Nil 28 

— > Nio 4 + 64Ni 

1.96X10 18 

-0.5 

42194. 

Ni 256 

— * Nij28 I28Ni 

2.91X10 18 

-0.5 

43263. 

Nigi 2 

-> Ni 256 + 256Ni 

4.36X10 18 

-0.5 

44109. 

Ni|024 

. — »• Nig, 2 + 5 1 2Ni 

6.61X10 18 

-0.5 

44780. 

Ni 2 Q48 

; ^Ni,o 24 + 1024Ni 

1.01X10 19 

-0.5 

45312. 

Table §.10b. Nickel Cluster 

based on Girshick’s formulas, Eqns. 

(§.7)-( §.13) 


Reaction 

A (cm /s/moles) 



P 

E/R (K) 


Ni 2 
Ni 3 
Ni 4 
Ni 5 
Ni 6 
Ni 7 
Ni 2 
Ni 2 
Ni 2 
Ni 3 
Ni 3 
Ni 6 
Nis 

Nig 
Ni 4 
Ni 8 
Ni,6 
Ni 32 
Niw 
Ni| 2 8 
Ni 256 
Nig i 2 

Ni)024 + 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


Ni 

Ni 

Ni 

Ni 

Ni 

Ni 

Ni 2 

Ni 3 

Ni 4 

Ni 3 

Ni 4 

Ni 2 

Ni 2 

Ni 3 

Ni 4 

Nig 

Ni, 6 

Ni 32 

Ni 64 

Ni i 2 g 

Ni 2 g6 

Nigi 2 

Nil024 


IN 12 

Ni 3 

Ni 4 

Nig 

Ni 6 

Ni 7 

Nig 

Ni 4 

Nig 

Ni 6 

Ni 6 

Ni 7 

Nig 

Ni 7 

Nig 

Nig 

Ni, 6 

Ni 32 

Ni 64 

Ni 128 

Ni 256 

Nig| 2 

Nii024 

Ni2048 


3.70X10 
4.10X10 
4.51X10 
4.90X10 
5.27X10 
5.61X10 
5.94X10 
4.16X10 
4.37X10 
4.60X10 
4.45X10 
4.59X10 
5.06X10 
4.83X10 
4.75X10 
4.67X10 
4.19X10 
9.41X10 
2.11X10 14 
4.74X10 14 
1.06X10 15 
2.39X10 15 
5.36X10 15 
1.20X10 16 


0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Table §.l 1 Reactions to form nickel/carbon nuclei clusters and nanotubes. 


! Nickel/Carbon Cluster Formation 


Z + Nil024 

= > 

ZNiC 

1.E16 

0.5 

0 

C60F + Nil024 

= > 

CFNiC 

1 .E16 

0.5 

0 

C70F + Nil024 

1 

= > 

CFNiC + . 1428571C7 OF 

1.E16 

0.5 

0 

! Nanotube Formation from Z and Ni 




Z + ZNiC 

= > 

0 . 08008008008CNT + ZNiC 

1.E12 

ID 

O 

0 

! Nanotube Formation from C60F, C70F and Ni 




C60F + ZNiC 

= > 

0 . 06006006006CNT + ZNiC 

1.E13 

0.5 

0 

C70F + ZNiC 

= > 

0 . 07007007007CNT + ZNiC 

1 .E13 

0.5 

0 

! Soot inerting 






Z 

= > 

DZ 

1.E2 

0. 0. 


ZNiC 

= > 

DZNiC 

1.E4 

0 . 

0 




Table §.12 Reduced carbon vapor/nickel catalyst model 


c 

+ 

C 

= 

C2 

2.00E+14 

0 

0 

c 

+ 

C2 

= 

C3 

2 . 00E+14 

0 

0 

C2 

+ 

C2 

= 

C3 + C 

2 . 00E+15 

0 

9040 

C3 

+ 

C 

= 

0 . 100CC 

2 . E+14 

0.00 

0.00 

C3 

+ 

C2 

= 

0.125CC 

2 .E+14 

0.00 

0.00 

C3 

+ 

C3 

= 

0.150CC 

2 .E+14 

0.00 

0.00 

1 

Nin Cluster 

Agglomeration 

Rao 

method 


Ni 

+ 

Ni 

= > 

Ni2 

3 . 70479E+12 0.5 

0 

Ni2 

+ 

Ni 

= > 

Ni3 

4 . 09656E+12 0.5 

0 

Ni3 

+ 

Ni 

= > 

0 . 0039062NiC 

4.51E+12 

0.5 0 


Ni3 

+ 

Ni2 

= > 

0. 0048828NiC 

4 . 37E+12 

0.5 0 


Ni3 

+ 

Ni3 

= > 

0. 0058593NiC 

4 . 45E+12 

0.5 0 


Ni 

+ 

NiC 

= > 

1.0009766NiC 

8.04E+13 

0.5 0 


Ni2 

+ 

NiC 

= > 

1 . 0019531NiC 

5 . 96E+13 

0.5 0 


Ni3 

+ 

NiC 

= > 

1. 0029297NiC 

5.03E+13 

0.5 0 



! Rate estimated from WebElements Ni-Ni bond energy 203 kJ/mol 
Ni2 + M => Ni + Ni + M 1 . 00E+15 0 24476. 

! Girshick Evaporation Rates with WebElements energy 


Ni3 

=> Ni2 + Ni 

3 . 37E+17 

-0.5 

28145. 

NiC 

=> 0 . 5NiC + 5 12Ni 

6.61E+18 

-0 . 5 

44780. 

! Carbon cluster agglomeration 




CC + C 

= 1.025CC 

2 . 300E+14 

0. 

0. 

CC + C2 

=> 1.05CC 

2 . 300E+14 

0. 

0. 

CC 

=> 0.95CC + C2 

3 . 2E13 

0. 

61900. 

! Fullerene 

formation 




CC + C3 

=> 0.70C60F + C 

2 . 00E13 

0. 

0 . 

CC + C2 

=> 0.70C60F 

4 . 00E+09 

0. 

-30196. 

CC + C 

=> 0 . 6833333C60F 

2.00E13 

0. 

0 . 

! Fullerene 

disintegration 




C60F 

=> 1.45CC + C2 

8.0E13 

0. 

61900 . 

C60F+C3 

=> 0.7875Z 

2 . 00E+13 

0. 

10065 . 

! Soot formation 




CC + CC 

=> Z 

4.E13 

0. 

0. 

Z + C3 

= 1.0375Z 

4.E12 

0. 

0. 

Z + C2 

= 1.025Z 

4.E12 

0. 

0. 

Z + C 

= 1.0125Z 

4.E12 

0. 

0. 

! Nickel /Carbon Cluster Evaporation 



ZNiC => 0. 

9990234ZNiC + Ni + 0.0009766Z 6.61E+18 0. 

44780. 

! Soot inerting 




Z 

=> DZ 

1 . 0E2 

0 . 

0 . 

ZNiC 

=> DZNiC 

1 . 0E4 

0 . 

0 . 

! Nickel /Carbon Cluster Formation 




Z + NiC 

=> ZNiC 

1.E16 

0.5 

0 . 

! Nano tube 

Formation from Z and Ni 





Z + ZNiC => 0 . 08008008008CNT + ZNiC 1.E12 .5 .0 

! Nano tube Formation from C60F and Ni 

C60F + ZNiC => 0 . 06006006006CNT + ZNiC 1.E13 .5 .0 


Table §.13 Comparison of the are process with other high temperature processes used for SWNT 
synthesis 


Process 

Conditions 

Energy 

density 

(W/m 2 ) 

Vaporization 

Flux 

(g/h) 

Soot 

containing 
SWCNT rate 
(g/h) 

nanotubes 
(% vol.) 

Temperature 

and 

cooling rates 

Arc discharge a 

/ = 100 A; V= 38 V 
r,\— 3 mm; dac = 3mm 
He; P = 660 mbar 

1 .4x 1 0 8 

72 

4.2 

50 % vol. 

r=6500K 
Cooling rate = 7300 
K/ms 

/ = 100 A; V —22 V 
r t \- 3 mm; DAC = 3 mm 
Ar ; P = 1 00 mbar 

7.8x1 0 7 

8 

1.5 

20 % vol. 

T= 1 8000 K 
Cooling rate = 4800 
K/ms 

Pulsed laser b 

A/ = 532 nm 
Epuisr = 250 - 490 mJ 
Vj — 10 — 30 Hz 
</>; spot = 5-6 mm 
tpuhr = 40-50 ns 
A 2 = 1064 nm 
Epu.sc = 300 - 550 mJ 
02 spot = 6-7 mm 
Ar ; P = 500 mbar 

1.3x10 s - 

5.2x10 s 

0.05 

3.3x1 0' 3 - 40x1 0' 3 

60 - 90 % vol. 

Oven at 1200°C 
Target at 3500 K 
Cooling rate = 170 - 
420 K/ms 

Continious 
laser c 

A = 10.6 fim 
0s pot = 1 mm 
Ar;P = 270-530mbar 

4.3x1 0 s 

0.1 

0.13 

80 % vol. 

Oven at 1200°C 
Target at 3500 K 

Solar energy d 

X = sun light 
oven at 2 kW 
power = 1 300 W 
0s pot = 1 5 mm 
Ar ; P = 250 mbar 
Oven at 50 kW 
He ; P = 450 mbar 
Ni/Co = 2:2 % at. 

7.4xl0 6 

8.5x10 s 

1.77x1 O'* 
15 

0.1 

13 

_ 

Target at 3550 K 
Cooling rate = 600 
K/ms 


a: Hinkov [116] 
b: Thess et al., [64] 

c: Munoz et al. [171]; Foutel-Richard [172] 
d: Guillard et al. [173]; Laplaze et al. [174] 


Table §.14 First Ionization potential and energy of elements used in models 



Y 

Ni 

C 

Ar 

He 

1 st ionization potential (eV) 

6.38 

7.63 

11.26 

15.76 

24.87 

Ionization energy (kj/niol) 

600 

737 

1086 

1520 

2372 


Table §.15 Electronic gas phase reactions considered in Hinkov [116] 


number 

Reaction 

A (cm' 3 .s.mor 1 ) 

JL 

E(K) 

Ref. 

Ion 1 

C + e C + + e + e' 

3.3x1 0 20 

0 

214420 

82 

Ion 2 

c + + c 3 -* c 4 + 

9.0x1 0 9 

0 

0 

82 

Ion 3 

c + + c 4 -> c 5 + 

6.0xl0 14 

0 

0 

82 

Ion 4 

C + + e C 

3.6xl0 1(> 

-4.5 

0 

85 

Ion 5 

C60F <-> C60 + + e 

1.33xl0 15 

0 

65900 

84 

Ion 6 

Ni ’+ e‘ Ni 

2.2x1 0 40 

-4.5 

0 

85 

Ion 7 

Y h + e' Y 

2.2x1 0 40 

-4.5 

0 

85 

Ion 8 

He + e’ <-» He 

2.2x1 0 40 

-4.5 

0 

85 


Table §.16 Surface reactions considered in Hinkov [116] 


number 

Reaction 

A (cm '.s.mol 1 ) 

P 

E(K) 

Deposition 

s, 

Ni + WL(S) WLNi(S) (Stick) 

1.0 

0 

0 

s 2 

Y + WL(S) — > WLY(S) (Stick) 

1.0 

0 

0 

S 3 

C + WL(S) ^ WLC(S) (Stick) 

1.0 

0 

0 

S4 

C 2 + WL(S) -4 WLC 2 (S) (Stick) 

1.0 

0 

0 

S 5 

C 3 + WL(S) WLC 3 (S) (Stick) 

1.0 

0 

0 

So 

C 4 + WL(S) WLC 4 (S) (Stick) 

1.0 

0 

0 

S 7 

C 5 + WL(S) WLC 5 (S) (Stick) 

1.0 

0 

0 

etching 

Sg 

C + WLC(S) — >WL(S) +C 2 (Stick) 

0.1 

u 0 l 

0 

recombination of ions 

S 9 

Ni + + e' Ni (BOHM) 

0.4 

0 

0 

S 10 

Y + + e ’ Y (BOHM) 

0.4 

0 

0 


Table §.17 Plasma characteristics calculated by AURORA assuming LTE and T K =T jon =T c =6500 K 


OUTLET CONDITIONS: 

Specified inlet mass flow rate 
Outlet mass flow rate 

(which, based on an reactor density 
and on a reactor volume 
produces a residence time) 
Total ion current density 


= 0.116 g/s 
= 0. 1 16 g/s 
= 0.55 1 xIO 05 g/cm 3 
= 0.6 cm 3 

= <).284xl0 04 s 
= 10.4067 A/cni 2 


Outlet and reactor temperature 

Outlet and reactor pressure 

Outlet and reactor density 

Outlet and reactor mean molecular weight 

Outlet molar flow rate 

Outlet volumetric flow rate 

(based on reactor pressure and temperature) 


6500 K 
0.651 atm 
0.55074x1 O’ 05 g/cm 3 
4.5101 g/mol 
0.25773x10°' mol/s 
21106 cm 3 /s 


= 37833 SCCM 
= 37.833 SLPM 


Outlet and reactor electron temperature Te = 6502.7 K 

Te = 0.56036 eV 


OUTLET CONDITIONS FOR GAS PHASE MOLECULAR SPECIES: 


Species 

Mole fraction 

Number density 
cm' 3 

mol/s 

g/s 

cm '7s 

e' 

8.7x1 O' 04 

6.4xl0 +l4 

2.2x1 O' 05 

1.2x1 O’ 08 

1.8xlO +01 

c 

4.5x1 O' 02 

3.3xl0 +lb 

1.2x1 O' 03 

1.4x1 0 02 

9.5x1 0 +02 

c 2 

2.0x1 a 05 

1.5xl0 +l3 

5.2x1 O' 07 

1.2x1 O’ 05 

4.2xlO’ 01 

C 3 

2.7x1 O’ 08 

2.0xl0 +1 ° 

7.0xl0’ 10 

2.5x1 O’ 08 

5. 8x1 O' 04 

c 4 

6.5x1 0‘ 16 

4.7x1 0 + ° 2 

1.7xl0' 17 

8.0xl0’ 16 

1.4x10'" 

c 5 

2.0xl0' 15 

1.5xlO +03 

5.2xl0' 17 

3.2xl0 15 

4.3x10" 

He + 

5.4xl0’’ 3 

4.0xl0 + ° 5 

1.4xl0' 14 

5.6xl0’ 14 

l.lxlO’ 08 

Ni + 

2.7x1 O’ 04 

2.0xl0 +l4 

7. lxl O' 06 

4.2x1 O’ 04 

5.8xlO +()0 

Y + 

4.7x1 O’ 04 

3.5xl0 +l4 

1.2x1 O' 05 

l.lxlO’ 03 

9.9x1 0 +0 ° 

C + 

1.3x1 O' 04 

9.2xl0 +13 

3.2x1 O’ 06 

3.9x1 O’ 05 

2.7x1 0 + °° 

C 4 + 

1.1x10 ,2 

7.9x1 0 +05 

2.8xl0’ 14 

1.3xl0’ 12 

2.3x1 O’ 08 

C 5 + 

1.7xl0' ls 

1.2xl0 +03 

4.4xl0’ 17 

2.6xl0 15 

3.6x10'" 

C 60 + 

5.3xl0' 41 

3. 9x1 O’ 23 

1.4x1 O’ 42 

9.9x1 O’ 40 

l.lxlO’ 36 

Ni 

1.7x1 O' 03 

1.3xlO +l4 

4.4x1 O’ 05 

2.6x1 O' 03 

3.6xlO +01 

Y 

4.6x1 O' 06 

3.4xl0 +12 

1.2x1 O’ 07 

l.lxlO’ 05 

9. 8x1 O' 02 

He 

9.5x1 O' 01 

7.0xl0 +n 

2.5x1 O’ 02 

9. 8x1 O' 02 

2.0x1 0 +04 


Table §.18 Plasma composition calculated for 1=100 A in helium with p=660 mbar 


Species 

Mole fraction 

Number density 

^ -3 

cm 

mol/s 

g/s 

cm Vs 

e' 

1.9x1 O' 03 

1.4xl0 +l5 

4.9x1 0 05 

2.7xl0'° 8 

4.0xl0 +ui 

c 

4.5x1 0 02 

3.3xl0 +l6 

1.2x1 O' 03 

1.4xl0'° 2 

9.5xl0 + ° 2 

c 2 

2.0x1 0 05 

1.5xlO +13 

5.1xl0 07 

1.2x1 O' 05 

4.2x10'°' 

C 3 

2.7x1 0 08 

2,OxlO +lu 

7.0xlO' 10 

2.5xl0'° 8 

5.7x1 O' 04 

c 4 

6.3x10 16 

4.6x1 0 +O2 

1.6x10 17 

7.9xl0‘ 16 

1.3x10 11 

C 5 

2.0xl0‘ 15 

1.5xlO +03 

5.1xlO' 17 

3.1xl0' 15 

4. 2x10* 1 

He + 

7.0x1 O' 12 

5.1xlO +06 

1.8xl0' 13 

7.2xl0' 13 

1.5x1 O' 07 

Ni + 

1.3x1 O' 03 

9.4x1 0 +l4 

3.3xlO'° 5 

1.9xl0'° 3 

2.7x1 0 + °' 

Y + 

4.7x1 0 04 

3.5xlO +l4 

1.2x1 O' 05 

l.lxlO' 03 

1.0xl0 + °' 

C + 

1.3x1 O' 04 

9.2xl0 +l3 

3.2x10'°° 

3.9xl0'° 5 

2.7x1 0 + °° 

c 4 + 

l.OxlO' 12 

7.7xlO +05 

2.7xl0' 14 

1.3xl0' 12 

2.2x10'^ 

c 5 + 

1.6x10 15 

1.2xlO +03 

4.2xl0' 17 

2.5xl0' 15 

3.5x10'" 

C 60 + 

6.5xl0 41 

4.8x1 0‘ 23 

1.7x1 O' 42 

1.2xlO' 3iJ 

1.4x1 O' 36 

Ni 

7.2x1 O' 04 

5.2xl0 +l4 

1.8x1 O' 05 

l.lxlO' 03 

1.5xlO + °' 

Y 

4.2x1 O' 07 

3.1xlO + “ 

l.lxlO' 08 

9.6x1 O' 07 

8.9x1 O' 03 

He 

9.5xlO‘ ul 

7.0xl0 +l? 

2.5xl0'° 2 

9.8x1 O' 02 

2.0x1 0 +O4 


THERMO AIL 

AR IVTANOAR 1 G 300.000 20000.000 1000.00 

0 .251994 77E+0 1-0 . 22858382E-04 0 . 74655237E-08-0 . 90208251E-12 0 . 36063733E-16 
-0 .75393340E+03 0 . 42492779E+01 0 . 25670200E+01-0 . 44283132E-03 0 . 98495011E-06 
-0.88424193E-09 0 . 27879189E-12-0 . 75176371E+03 0 . 40801018E+01 
C GMcB 20K 0C 1 G 300.000 20000.000 1000.00 

0 . 23205176E+01 0 . 18268312E-03-0 . 31504542E-07 0 . 25371109E-11-0 . 62267852E-16 
0 . 84851595E+05 0 . 58203836E+01 0 . 19031864E+01 0 . 40200238E-02-0 . 90395391E-05 
0 . 81759092E-08-0 . 25854092E-11 0 . 84826660E+05 0 . 72917617E+01 
C2 GMcB 20K 0C 2 G 300.000 20000.000 1000.00 

0 . 40753780E+01 0 . 33239978E-03-0 . 12344909E-07-0 . 70714041E-12 0 . 28222317E-16 
0 . 99738114E+05 0 . 96847543E+00 0 . 5499364 1E+01 0 . 188914 70E-02-0 . 13226503E-04 
0 .16061423E-07-0 . 58286773E-11 0 . 99140682E+05-0 . 77146073E+01 
C3 GMcB 20K 0C 3 G 300.000 20000.000 1000.00 

0 . 4903274 6E+01 0 . 11113409E-02-0 . 75331984E-07-0 . 48099305E-12 0 . 89121604E-16 
0 . 96812289E+05 0 . 15052816E+00 0 . 3604201 4E+01 0 . 75373711E-02-0 . 11553984E-04 
0 . 88850242E-08-0 . 25287216E-11 0 . 97010617E+05 0 . 61103710E+01 
C4 GMcB 20K 0C 4 G 300.000 20000.000 1000.00 

0 . 49032746E+01 0 . 11113409E-02-0 . 75331984E-07-0 . 48099305E-12 0 . 89121604E-16 
0 . 11492855E+06 0 . 15052816E+00 0 . 3604201 4E+01 0 . 75373711E-02-0 . 11558984E-04 
0 . 88850242E-08-0 . 25287216E-11 0 . 11512688E+06 0 . 61103710E+01 
C5 GMcB 20K 0C 5 G 300.000 20000.000 1000.00 

0 .11592118E+02 0 . 85502931E-03-0 . 12817463E-06 0 . 77473555E-11-0 . 16260603E-15 
0 . 11368121E+06-0 . 35115388E+02 0 . 11731935E+02-0 . 25681572E-01 0 . 78642726E-04 
-0 .78483519E-07 0 . 26116987E-10 0 . 11495211E+06-0 . 29295593E+02 
HE 120186HE 1 G 0300.00 20000.00 1000.00 

0 . 02500000E+02 0 . 00000000E+00 0 . 00000000E+00 0 . 00000000E+00 0 . 00000000E+00 
-0 . 07453750E+04 0 . 09153489E+01 0 . 02500000E+02 0 . 00000000E+00 0 . 00000000E+00 
0 .00000000E+00 0 . OOOOOOOOE+OO-O . 07453750E+04 0 . 09153488E+ 01 
Ni GMcB 20K ONi 1 G 300.000 20000.000 1000.00 

0 . 36971 152E+0 1-0 . 82867716E-03 0 . 18628604E-06-0 . 11986930E-10 0 . 23752778E-15 
0 . 50060051E+05 0 . 72075154E+00 0 . 42140643E+01-0 . 10224755E-01 0 . 25272825E-04 
-0 .24064628E-07 0 . 78454681E-11 0 . 50323075E+05 0 . 58841215E-01 
Y IV 25K 0Y 1 G 300.000 25000.000 1000.00 

0 . 79437350E+00 0 . 13579116E-02-0 . 12031792E-06 0 . 31899976E-11-0 . 15582902E-16 
0 51410499E+05 0 . 18158732E+02-0 . 66965149E+01 0 . 65805915E-01-0 . 14851900E-03 
0 . 13342009E-06-0 . 41975353E-10 0 . 51184455E+05 0 . 45676826E+02 


1 

2 

3 

4 
1 
2 
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4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 
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Table §.19 Soot nuclei and soot particle reaction rates, from Krestinin and Moravskii 1 1 19). 


Reaction 

A(cm Vs/moles) 

P 

E 

1. Formation of soot nuclei Z 

c 78 + c 2 -» z 

4.0xl0 08 

0 

-60.8 

C 78 + C 2 — > z + c 

4.0x1 0 08 

0 

0 

n,m = 1-10 
C60-m + C60-n > Z 

4.0xl0 14 

0 

0 

C 70-111 + C 70 - 1 , — > z 

4.0xl0 14 

0 

0 

2. Heterogeneous reactions on soot particles 11 

Ci + soot 44 soot 4.0x10°’ 

0 

0 

C 2 + soot <-> soot 

4.0x1 0 03 

0 

0 

C 3 + soot -4 soot 

4.0x1 0 03 

0 

0 

n=4-79 

C n + soot -4 soot 

4.0x1 0 03 ( 1 /n ) 05 

0 

0 

Coof + soot -4 soot 

l.OxlO 03 

0 

30 

C 70 F + soot — » soot 

l.OxlO 03 

0 

30 


d The rate constants of heterogeneous reactions are given per unit aerosol surface area. 


Table §.20 Surface reactions considered in Farhat, et al. [1 15]. 


Reaction 

A unit (cm-s-moles) 

P 

E(K) 

(Si) CR+C — > CR+CNT 

2.50x10" 

0.5 

0.0 

(S2) CR+C3 CR +CNT+C2 

2.50x10" 

0.5 

0.0 

(S3) CR+C2 -> CR +CNT+C 

2.50x10" 

0.5 

0.0 


Table §.21 Calculated number densities of small carbon clusters and fullerenes at three positions in 
the arc. A symbol « -« indicates negligible density. Conditions are helium with p=660 mbar, 1=100 
amps, dAc =3 mm, and the dilution factor = 20. 



T 

(K) 

C density 
(cm') 

C ? density 
(cm") 

C3 density 
(cm' 3 ) 

C(,o density 
(cm' 3 ) 

C70 density 
(cm' 3 ) 

ANODE 

4000 

7.8xl0 1:> 

l.lxl0 +15 

l.lxl0 +l5 

- 

- 

1 mm from the 
cathode 

6472 

1.2xl0 +l6 

2.9xl0 + ' 2 

2.8xl0 +09 



CATHODE 

1837 

1.7xl0 +l4 

8.5xl0 +13 

6.0xl0 +l3 

2.6xl0 +12 

7.0x1 0 +09 


Table §.22 Simplified Carbon Chemistry Model used in laser ablation simulation 


Reaction Rate Coefficients (where k f =A*T A B*exp ( -T a /T) 



REACTION 


A, cm 3 /mole/s 

B 

T a ,K 

1 

C + e- 

<=> 

C+ +G- + 0~ 

5.46E-04 

0 

214420 

2 

C + C 

<=> 

C 2 

2.00E+14 

0 

0 

3 

C + C 2 

<=> 

C 3 

2.00E+14 

0 

0 

4 

C 2 + C 2 

<=> 

C 3 + c 

2.00E+15 

0 

9040 

5 

O 2 + O 2 

<=> 

C 4 

2.00E+14 

0 

0 

6 

C + C 3 

<=> 

c 4 

2.00E+14 

0 

0 

7 

0 + O 4 

<=> 

C 5 

2.00E+14 

0 

0 

8 

C 2 + C 3 

<=> 

Cs 

2.00E+14 

0 

0 

9 

0 + + C 3 

<=> 

C 4 + 

9.00E+09 

0 

0 

10 

0 + + O 4 

<=> 

C 5 + 

6.00E+14 

0 

0 

11 

0+ + ©- 

<=> 

c 

3.60E+16 

-4.5 

0 

12 

C5 + C 

<=> 

C 6 

2.00E+14 

0 

0 

13 

C 4 + G 2 

<=> 

C 6 

2.00E+14 

0 

0 

14 

C 3 + C 3 

<=> 

C 6 

2.00E+14 

0 

0 


k, 1/s 



n 


Fig. 1 Comparison of Iron Cluster Evaporation Rate Coefficients at 1380 K. 
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Fig. §.2 Iron cluster recombination rates at 1500 K. Iron atoms attaching to clusters of 
size n-1 forming clusters of n-atoms. 
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Fig. §.3 Nanotube bundles emanating from a nickel/cobalt catalyst particle. (Courtesy A. 
P. Moravsky) 
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Fig. §.5 Published in JNN Vol. 4 No. 4 (2004) Fig. 3 
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Figure §.7. 
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Figure §.9. 
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Figure §.10. 
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Figure §.11. 
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Figure §.14. 
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Figure §.17. 
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Figure §.20. 
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Figure §.21 -a. 
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Figure §.21-b. 
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Figure §.22. 
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Figure §.23. 
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Figure §.24. 
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Figure §.25. 
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Figure §.26. 
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Figure §.27. 
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Figure §.28. 
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Figure §.29. 
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Figure §.30. 


31 


15 cm 


◄ 




ANODE 
Axis 3800 K 
Side 3300 K 


1 or 4 mm 


CATHODE 
Axis 3800 K 
Side 3300 K 



t 


(i) 



Arc 
region 

( 2 ) 

Wall — 1 

T=350 K 

i 


▲ 


100 cm 



Figure §.31. 
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Figure §.32. (to be improved) 
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Figure §.34 
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Figure §.36 
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Figure §.37 
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Figure §.40 
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Fig. §.43 Streak lines in space for laser ablation axisymmetric calculation 
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Fig. §.44 Temperature histories along selected streak lines for laser ablation calculations 
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Fig. §.45 Example of evolution of species mass fractions along streak 
line ■in laser ablation. Arbitrary evaporation rate coefficients were 
used in calculation, therefore, data is not accurate, but only shown as 
an example of methodology. 
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